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Abstract 

Theoretical models of seismic sources are developed for 
digital data collected mainly at Edmonton and Mexico city. 
This thesis uses the first digital data ever recorded for a 
major earthquake sequence in Mexico and is the first attempt 
to do a study of source parameters in South West Alberta. 

First, source theory and some empirical relationships 
in seismology are reviewed. Then the source parameters of 
the aftershocks of the Oaxaca, Mexico earthquake of November 
29, 1978, are discussed. The digital data used here was 
recorded at an epicentral distance of approximately 500 km 
and is the first to be recorded in a set of permanent 
telemetering digital seismic stations in central Mexico 
(RESMAC). These source parameters are compared with those 
calculated using data recorded by digital stations located 
above the source region. 

The source parameters of the Rocky Mountain House 
earthquake swarm were studied using data recorded by the 
Edmonton (EDM) digital station. Even though a typical RESMAC 
Station, and Edmonton station have different amplitude 
response curves, and the epicentral distances of the events 
were different, the bandwidth of useful information in both 
systems was between 0.5 and 7 Hz. 

Some events of the Rocky Mountain House earthquake 
Swarm were also studied with a portable digital seismic 
Station. The data obtained from this station was of much 


higher quality than the data from RESMAC or the Edmonton 
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digital station. Finally, seismic data related to the local 
Seismicity of South West Alberta recorded at Edmonton since 
1970 are used to make tentative conclusions about the 


present day dynamics of this area. 
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1. Introduction 
Geophysical data, such as observations made on earthquakes, 
have been described, only partly facetiously, as of type 
III; Inaccurate, Incomplete, and sometimes Incorrect. The 
reason for this lies in the fact that a controlled 
experiment in earthquake seismology is impossible and that 
the configuration of observing equipment is only partly 
designed with basic research in mind. The underlying problem 
in earthquake seismology is then to analyse Type III data 
sets and to extract from them scientifically valid 
conclusions. In this work I investigate 3 such data sets, 
Show that limited conclusions about the nature of the 
mechanism that generates the events are possible, and that 
these conclusions have geological implications. 

With the exception of California and Japan, earthquakes 
do not often happen in areas well instrumented to study the 
response of the earth to this stimulus. Much of our 
knowledge of these phenomena is thus derived from 
observations made at a distance combined with models of 
processes at the source. The verification of the legitimacy 
of these models and the determination of their free 
Parameters is the concern of observational seismology. The 
data on which such investigations are based are the records 
of vibrations induced at various seismic stations around the 
world by the earthquakes. 

Effective analysis of these records is limited by the 


means available to process the data. Until recently the 
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recording medium for the data was photographic paper. 
Although it was relatively easy to determine the time of 
arrival of the wave on such a record a Systematic study of 
the spectral content of the signal was truly painful. Even 
computer aided digitization schemes required numerous 
corrections for unfortunate properties of the data recording 
scheme. Surprisingly, modern versions of the original 
seismograph, a sharp stylus scratching smoked paper, are 
still in fairly common use. 

The appearance of technology which allowed digital 
recording of seismic data (see Aki and Richards Section 
11.1.2 1980) changed the nature of seismic investigations 
dramatically. It now became possible to manage, with the aid 
of computers, very large sets of data. In particular the 
Spectral characteristics and the nature of the radiation 
pattern from earthquakes could be investigated. The work 
done on the LASA network (see for example Engdahl et al 
1970) is one example of the processing that evolved. 
Recently the appearance of cheap digital technology has lead 
to the spread of these techniques from large expensive 
installations to smaller research operations. 

In spite of the developement of technology, the 
distribution of seismic stations in many parts of the world 
is far from optimal for investigation of the spectral 
characteristics of the radiation pattern of seismic events. 
This does not mean data acquired at such installations is 


useless, merely that there are limitations to its value. 
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Since many such installations now exist in the world it is 
useful to ask not only what observations on these 
installations tell us about the spectral characteristics of 
the seismic source, but also to what degree the limited 
scope of the installation limits the value of its data for 
such investigations. 

I discuss here the analysis of 3 sets of such data. The 
first objective of the work was to determine in so far as 
waS possible the nature of the source that gave rise to the 
seismic signal. Clearly understanding of this source must 
have implications for the physical process causing 
earthquakes and a second objective was to connect the 
measured characteristics of the signal to acceptable models 
of the source. The parameters of these models hopefully will 
cast some light on the exceedingly complicated process 
generating seismic events. As a byproduct of the study I 
came to some conclusions about the relationship between the 
capabilities of the observing system and the resolution in 
the source that could be achieved. 

Studies of source parameters, such as those reported 
here, are of great importance in seismology because they can 
show regional differences in the earthquake generation 
process. In other words, regional tectonic processes can 
give rise to variations in the source related parts of the 
spectra of seismic signals. Evidence of such changes is 


contained in the radiated elastic wave energy. 
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One of the more powerful techniques for such studies is 
the analysis of the spectra of body and surface waves in the 
near and far field. Study of the spectra in the far-field 
give uS average properties of the source, such as seismic 
moment (strength of the fault), and average stress drop. On 
the other hand, the study of near-field spectra gives us 
details about the rupture process. Such details include 
complexities of the rupture process, energy focusing, 
maximum and minimum displacements, and velocities and 


acelerations as a function of azimuth. 


1.1 The Data Used 
The data sets I used are 

1. The aftershocks of the Oaxaca (November 29, 1978) 
earthquake as recorded on a digital station 500km from 
the epicentral zone. 

2. The events of a swarm of small magnitude earthquakes 
near the Canadian Rockies as recorded on a 
digital/analogue station 180km from the epicentral area. 

3. And the records of this same swarm recorded on portable 
digital seismic recorders above the swarm. 

This data constituted a substantial, unexploited digital 

seismic data set. It originated from RESMAC in Mexico and 

Edmonton station in Canada. 

At the end of 1978, when the Oaxaca earthquake 


occurred, RESMAC had only three stations in operation, 
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Cerrillo (CRX), Mexico city (MEX), and Acapulco (ACR). In 
addition to the RESMAC stations, the SISMEX stations 
Tonantzint la) i(1tT):,5 Santar Rita (ITC). "and El Pino (1IP), 
provided additional information. Typically with a limited 
network of stations, such as this, epicentral locations near 
the shock cannot be significantly improved. It is possible; 
however, to do spectral analysis of body waves of the 
aftershocks of the Oaxaca earthquake, and compare our 
results with the spectral analysis from the portable digital 
stations deployed by others above the aftershock area. The 
results of this study are described in chapter two. 

Edmonton is the only permanent station in Alberta that 
records in analog and digital format in its Six components 
(short and long period), and the only one in Western Canada 
that records in digital format. It was designed to record 
teleseismic data (McCloughan and Kanasewich 1974), it 
records 18 samples per second in the short period and three 
samples per second in the long period signal. It has an 
amplitude response curve typical of the World-Wide Standard 
Seismograph Network. 

All useful information from this station is archived on 
magnetic tapes at the Department of Physics of the 
University of Alberta. The information recorded at Edmonton 
from the South West Alberta of the Rocky Mountains and 
particularly from the Rocky Mountain House earthquake swarm 
allows 7a Study of this ¢ activity.«Epicentral ;locations are 


generally done by Earth Physics Branch, Department of Energy 
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Mines and Resources in Ottawa, consequently, in order to 
improve those epicentral locations a network of stations 
above the epicentral area is desirable. 

I did the spectral analysis of the body waves of the 
Rocky Mountain House earthquake swarm as observed at 
Edmonton, in order to investigate the earthquake source, and 
the analysis of the refracted shear wave Sn in order to 
calculate focal depths of those events showing clear Sn 
phases in the digital records. 

In chapter four I analyse a fraction of the Rocky 
Mountain House earthquake swarm recorded on a Spregnether 
DR-100 digital instrument. The results obtained are far 
better, and the source parameters are more representative of 
the active area; However, in order to improve this 
information more than one station above the earthquake 
source is necessary. Then it would be possible to obtain 
fault plane solutions and determine principal stress 


orientations in this area. 


1.2 Seismic Source Models 

All these studies explicitly or implicitly require a 
model of the earthquake source. In general the source model 
consists of the intuitively reasonable assumption that over 
a finite surface, the fault plane, within the earth 
(modelled as perfectly elastic) there exists a discontinuity 


in displacement. By analogy with usage in solid state 
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physics this is generally known as a dislocation even though 
the scale of the processes in seismology is much larger and 
Pays no heed to atomic structure. In the more general case 
there can also be a discontinuity in the traction on the 
fault plane. The resulting displacements in the model can be 
calculated using the representation theorem (Aki and 


Richards 1980 equation 2.41) 


Soc t = | at | | | fees EG eto: SPOuav Co) 


" | dt | | 1G, ,(E,t-1;%, 0), (W(E,t) F) 
—©o S 


~ Up (E/T) es 54g E)njG,, 9 (E,t-1 7%, 0) } as (z). 


If the size of the fault plane is small compared with the 
epicentral distance to the detector the displacements can be 
Separated into terms that decay like r~' and terms that 
decay as higher powers of r (r is the epicentral distance). 
Displacements that decay like r~' are called far field 
displacements and terms that decay at higher powers of r are 
called near field (see Aki & Richards 1980 page 73). One 
characteristic of the far field is that it can be separated 
into P and § displacements. In the near field (distances of 
the order of the fault length) we have mixed P and S body 


waves. 
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The problem in source studies is to characterize the 
time evolution of the discontinuities on the fault plane. 
Obviously additional physics is needed for there is nothing 
in the theory that indicates when a discontinuity will 
appear (an earthquake will start) or when it will stop. 
Since the physical properties of faults are largely unknown 
seismologists assume time histories and fit them to spectral 
properties of the radiated elastic energy or to the actual 
Signal shape as observed on seismographs. 

There is a useful simplification of the representation 
theorem. If we introduce the concept of a moment tensor to 
describe the source, the representation theorem acquires the 


form (Aki and Richards 1980 equation 3.22) 


ial a pq np,q 


In general nine couples (called vector dipoles) are 
required to obtain equivalent forces for an arbitrary 
displacement discontinuity in an anisotropic medium. In 
other words, the equivalent forces corresponding to an 
infinitesimal dislocation can be represented as a 
combination of those nine vector dipoles. Those nine dipoles 
together constitute the seismic moment tensor and have units 


of moment per unit area (dynes-centimeter). The components 
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of the seismic moment tensor for a constant slip on an 


infinitesimal fault in cartesian coordinates are given by: 


MP =- M5 (Sis 6 COSA: SiG 26. eS Ss Su) IN ahi Delis 
Oy ar Maa(Siin OtCOSeX ‘cos 26.°+8L Sin 2 SsiniX Sin 26) 2° mM 
y fo) S Ss yx 
, 
Nat =- Ms (cos 6 cosS-A cos 3. + COS) 2067; SkN7A «Sia o,) = Mee, 
Megtsin 6° CoSsoN “Sins 267 Sosin 26 sin ) 2 0) 
yy a C cos Bie 
Me =- mM, (COS G-COS.N Sia %. =—4COS@ 126 SiMe, COS 3.) = tye 
Mes = M Si me’ Sasi md 


Where M, is the seismic moment, , is the strike of the 
fault, » is the rake (rake is the angle between strike 
durection pand.sivp) ped. pusythecdip sand M..) .are,.the 
components of the seismic moment tensor. The constant term 
M, is the static seismic moment and gives us the strength of 
the fault. The dislocation is a function of time, therefore 
the seismic moment tensor is also a function of time since 
it depends directly on the average dislocation along the 
fault. 

The problem of observation interpretation then becomes 
the calculation of properties of the moment tensor and the 
further interpretation requires connection of this moment 
tensor to models of the discontinuities. Steketee (1958), 
Maruyama (1963), Burridge and Knopoff (1964), Brune 
(1970,1971), among others developed a general theory for the 


faulting process in an elastic medium. They assumed that the 
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rebound theory is valid for shallow earthquakes, and 
calculated the displacements using the representation 
theorem. Although their work preceeded the formal statement 
of the moment tensor representation it is easy to derive 
formulas to connect moments with discontinuities. Excellent 
reviews of source theory can be found in Archambeau 
(1968,1975), Ben-Menahem and Singh (1972), Johnson (1979), 
and Sahay (1980). 

Seismologists usually follow either a kinematic or a 
dynamic approach. In the kinematic approach the earthquake 
process is specified by the starting size of the 
discontinuity, the velocity at which slip spreads over the 
fault, and final dimensions of the region over which the 
Slip has occurred. In the dynamic approach to the earthquake 
process, the earthquake is modelled as a shear crack. The 
crack nucleate in a preexisting stress field which causes 
stress concentrations around the tip of the crack which in 
turn cause the crack to grow. Time history of the rupture 
velocity of the crack as well as self-similarity are assumed 
(If a function u(x,t), the solution of the shear dislocation 
for example, is a homogeneous function of x, t of degree 1, 
we say that we have a self-similar solution). 

Given such assuptions the wave equation with appropiate 
boundary conditions outside the source region can be solved 
by developing a Green's function (see Aki & Richards 1980 
sec’ 4.2). This Green's function can be convolved with the 


source term in order to obtain P and S displacements. 
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Clearly such displacements carry information about fault 
Orientation, source strength (seismic moment), and have a 
directivity pattern or radiation pattern. 

Seismologists do not usually have sufficient data to 
calculate all properties of the seismic moment tensor. Often 
they are restricted to determining M,, the seismic moment. 
Several methods can be used to calculate the seismic moment. 
Geodetic measurements estimate the average dislocation at 
infinite periods (D), hence knowing the fault surface, A, 
(for example from the aftershock area) and assuming the 
Shear modulus ( : ), the moment is given by M,=uDA. However, 
a more common method for calculating the seismic moment is 
from the spectra of the far-field displacement of body 
waves, surface waves or free oscillations of the earth. 

The spectrum of the body waves gives the moment at 
periods of the order of seconds. Therefore this moment is 
more affected by anomalies in the spectrum resulting from 
local complexities of the earth structure which give rise to 
Scattering related energy losses. The moment calculated with 
long-period surface waves or free oscillation data is less 
affected by structure complexities, since it is calculated 
at periods of minutes and hours. Values of M, range from 
about 10°° dyne-cm for very large earthquakes to 10'? 
dyne-cm for micro-earthquakes, and 10° dyne-cm for 
microfractures laboratory experiments in rocks (Aki & 


Richards page 49). 
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Observational studies are directed primarily at 
obtaining the data necessary to either check the theoretical 
Studies or investigate the consequences of interpreting 
earthquakes in terms of a specified model. This process can 
be done in two different ways; in the time domain, by 
adjusting the source parameters until a reasonable agreement 
is achieved between synthetic and observed seismograms, and 
in the frequency domain by inferring those characteristics 
of the spectra such as the long period level, high-frequency 
asymptote, and corner frequencies, that can be directly 
bédavea to the source properties. 

Brune's model (1970, 1971) follows this last approach. 
It illustrates both the virtues and difficulties of this 
problem. The Brune model is simple enough to be a reasonable 
and tractable representation of some properties of an 
earthquake. It is in fact remarkably successful. 
Nevertheless it is by no means a complete description. The 
properties of an earthquake derived by assuming Brune's 
model to be correct are best thought of as a more compact 
and transparent representation of the seismic observations. 
Their connection to real stoees se requires further 
analysis. Perhaps the current state of this science is best 
illustrated by the undoubted fact that most observations can 
be fit to such a simple model in spite of its obvious 
physical deficiencies. 

An earthquake source can be modelled (Brune 1970,1971) 


as a tangential stress pulse applied instantaneously to the 
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interior of a dislocation surface (this implies a 
nonphysical infinite rupture velocity). The total shear 
Stress available to accelerate the two sides of the fault is 
the difference between the initial stress and the dynamic 
frictional stress which is of the opposite sense to the 
initial stress and always acts to resist the fault slip. If 
the movement across the fault is such that the shear stress 
acting across the fault has decreased to the value of the 
frictional stress (Orowan 1960 see also Pilant 1979 page 
392), then the stress drop is equal to the total or 
effective stress drop, otherwise we have a fractional stress 


drop. The fractional stress drop is defined by Brune as: 


where o, is the initial stress, c, is the dynamical frictional 


1 


Stress, aiivssthe: tinabtsuress ( 9,204) and is the total or 


effective stress drop. 

The moment, source dimensvonk and fractional stress 
drop for the proposed model can be related to the spectra of 
the far-field displacement of an SH body wave. The long 
period level of the spectrum is proportional to the seismic 
moment, and the intersection of the long period level and 
the asymptote at high frequencies defines the corner 


frequency which is related to the size of the fault. This 


spectrum corresponds to total stress drop. If we do not have 
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a total stress drop (&<))} ‘the long period) level of the 
Spectrum 1s reduced to « times its value for a total stress 
drop. 

Brune obtained the corner frequency by considering a 
balance of energy around the dislocation surface. This 
corner frequency in the spectrum can arise from finite 
energy flow away from the dislocation surface, from finite 
rupture velocities or from finite length of the fault. 
Later, Hanks and Wyss (1972) extended Brune's theory 
empirically to include P-wave spectra for non-circular 
faults. 

The decrease in the amplitude spectra at high 
frequencies proportional to w~? and w™* was first proposed 
by Aki (1967). In the Haskell (1964) model falloff is as 
w'*. In kinematic models the effect of finiteness and 
rupture velocity of the source time function introduces a 
smoothing factor of w ' at "intermediate" frequencies given 
by the envelope of sin(X)/X (see Aki and Richards 1980 sec 
14.1.5). The effect of finite rise time introduces another 
factor proportional to w™'. Thus, the spectra have a flat 
part at low frequencies, a corner frequency, an w™' 
dependence at "intermediate" frequencies, and a w~* decrease 
at high frequencies. Asymptotes in the spectra at high 
frequencies have been observed in the range from 1 to 4. 

The spectrum of SH reflects source properties since SH 
does not suffer conversion of phases at sharp interfaces. 


However, we need to rotate the horizontal components (NS and 
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EW) in order to obtain the horizontal components of motion 
transverse to the azimuth from the epicenter to the station. 
Corrections applied to SH spectra are: free surface 
amplification, (for! SH this is 2 and for SV it) should be 
calculated, see for example Nutlli 1961), radiation pattern, 
instrument response, and attenuation, Q (generally taken as 
independent of frequency). 

Other earthquake source models (Haskell 1964, Savage 
1972, Dahlen 1974, Burridge 1975) are discussed in the next 
section, and compared with the observed results obtained 


with Brune's model. 


1.2.1 Studies related to Brune Model 

Theoretical models which assume a rupture velocity 
predict corner frequencies as a function of azimuth 
(Madariaga 1976, Vargas et al 1980). Theoretical models of 
long narrow faults predict similar corner frequencies for P 
and S wave spectra. Yessie of equidimensional faults predict 
greater corner frequencies for P waves than those for §S 
waves. Molnar, Tucker and Brune (1973) analysing the spectra 
of P and S-waves of the aftershocks of the San Fernando 
earthquake, concluded that their observations were 
consistent with a model which treats earthquakes as 
approximately equidimensional. 

Savage (1972) comparing the spectra of the Haskell 
model (1964) with Brune's model found that they have similar 


relations between the fault dimension and the corner 
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frequency in the displacement spectrum. One of the 
characteristics of the spectra derived from the Haskell 
model is that P and S-wave corner frequencies are similar as 
long as the rupture velocity is subsonic, i. e. with a 
dislocation velocity less than the S-wave velocity. Later 
Savage (1974) explained that the Haskell model he considered 
in 1972 can predict that P-wave corner frequency may exceed 
the S-wave corner frequency if near-sonic or supersonic 
rupture velocities are considered. Laboratory experiments 
along pre-existing shear faults show rupture velocities from 
0.25 to slightly excceding shear velocities (Shamina et al 
1978, Vinogradov 1978). 

Dahlen (1974) also considered a kinematic model in 
order to explain the observed ratio of P to S corner 
frequencies obtained uSing Brune model and later empirically 
extended by Hanks and Wyss. He assumed that the rupture 
nucleates at a single point on the fault surface and spread 
subsonically. He found that the ratio of P to S corner 
frequencies is less than unity at all points of the focal 
sphere. 

Burridge (1975) modified Dahlen's model for shallow 
earthquakes in order to analyse the controversial frequency 
shift of the corner frequency. He considered that the 
rupture nucleated at some point and began to propagate at 
the P speed, slowing down in such a way that the initial 
break governs the high-frequency content of the spectrum. 


This model, prediets that in 70 per cent of the focal’ sphere 
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the P corner frequency is higher than that of S while in the 
remaining 30 per cent, S corner frequencies are higher than 
P corner frequencies. This 1s more or less in agreement with 
the observations, see for example Molnar Tucker and Brune 
ro 3) 

Brune et al (1979) analysed the effects of bilateral 
and unilateral dislocation ruptures on the far-field 
displacements. They found that the pulse shapes are strongly 
dependent on the azimuth. This is equivalent to different 
corner frequencies at different azimuths 

Tucker and Brune (1977) analysed the spectra of the 
aftershocks of the San Fernando earthquake uSing a great 
variety of instruments at different hypocentral distances. 
They found events with one and two corner frequencies. 
Events with a single corner frequency were consistent with 
faulting confined to about 0.1 sec in time or 500 m in 
Space. About half of the larger events had two corner 
frequencies which are consistent with an initial rupture 
that grows during about 2 sec and to about 6 km. They 
explain the events with two corner frequencies as fore-slip 
with a high dislocation velocity followed by afterslip with 
a low dislocation velocity. The calculated stress drops for 
those events lie between 1 to 300 bars. 


Hanks (1981) in a review of 36 source mechanism studies 


comments, 


#03 AL etliiv & 
I 

Cai. 2hcGie 

iy iw Steines w: 


Sivas he 


= a ey HAL 
| : a i 
if ( nek 


Lee 


Veta oa Ty 2 39384 ats ay 2 Lh: ae: he st 
bieit-sat aes coy ee =r as ie no banat ib J oi 
Sada! Met wey site ma huge td ! Leda biabait: rot | at | a ‘ 
vers 1 jatigad up , ‘gi antes nie 
‘ cue be “atts ee Sh Te HES diese! | He: ‘ict : 

orid) ay Ste Be bade ris eres apie Bena 
Seene eo on iei re Shit Gis 7a are Be, phi mil Mes 


‘ER ‘pete ae + Ey SOT YL AE rene It AR ay Ye + ie ft ban ig | 
Lt Sh ee debs =k was bia Age Season 
Wie “Ateals bane. See recat jit aa By nm 

Ad, R08’ 44) Sanbthlyh ued Ripon a 
Ten Sows ser ve tiate Sewer od a ea ated Pest 4 


sadiget Leis se ott i) inedesoring: #96 dvi paid 7 
{403.1 te 6 Pinot ay bs (ie ae, ie duptta on! amoig 
riavedo? sainsz. sa18nips 3 12 Aehses eile rahi es sere utd. ud 

thie E ys og eR “a atu ROU we paeeate nib adoote tit ‘Hea 
3O> Gur Hh Heite He saint ae a ni eae HOR ie gee 
| spies aeipt Fagot ‘ond é 


RN Pep ig.ih ay 53% zeT, 3E ig, nai awk OC Oe 


an Wa 
j P ie om 


Y 
+ 


ie ie 
Ks igi i } fea 
ie Me ey 7 ar, 2 
Dore ~ Pin ra “ws 


"The corner frequency shift is a very common 

condition of the far-field body wave of earthquakes, 

with no discernible dependence on earthquake source 

strength, hypocentral distance, depth, or recording 

device; and the frequency shift is the manifestation 

of an intrinsic property of earthquakes, source 

finiteness." 
His conclusions constrain theoretical earthquakes source 
models. Greater P-wave corner frequencies than S-wave corner 
frequencies are in general observed. However, there have 
been some cases of smaller P corner frequencies. This may be 
due to the analysis of data from a single station which can 
be affected by complicated local effects (see for example 
Bakun et al 1976, Rebollar et al 1981). 

Fault slips for small events obtained using spectral 
analysis have been found to vary substantially with the 
distance at which observations are made. In deep mines they 
can be several orders of magnitude less than those measured 
directly (McGarr et al 1979). However, fault slip calculated 
from spectra of large events usually agrees with geological 
observations. 

Stress drops of the order of kilobars have been found 
in laboratory experiments on brittle fracture of rocks, 
while stress drops of natural earthquakes lie in the range 
of a few bars (1-100). McGarr et al (1979), found that 
tremors in deep mines appear to be associated with stress 
drops of the order of 1 kilobar. Consequently, they 
postulate that a certain amount of stress inhomogeneity is 


required before failure occurs. This may imply large 


concentrations of stress in localized regions. 
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Miyatake (1980a) investigated a three-dimensional 
dynamic crack model in infinite and semi-infinite 
homogeneous elastic model using the finite difference 
method. He considered different strength distributions on 
the fault and investigated the effect of fault strength on 
crack propagation. For low strength the shape of the rupture 
front is elliptic, with a rupture velocity of the order of P 
along the major axis and with a rupture velocity close to § 
along the minor axis. Relatively high strengths produce a 
nearly circular rupture front with a rupture velocity nearly 
equal to S. For a line barrier in the fault (like that 
proposed by Das and Aki 1977), he found that the rupture 
propagated backwards, generating multiple shocks. In the 
case of random strength distributions on the fault, highly 
irregular rupture propagation is generated and in some cases 
unfractured regions remain. 

Miyatake (1980b) analysed the near and far field 
displacement as well as the shape of the spectrum due to 
different distributions of strength on the fault. He found 
that complex seismic waves including multiple shocks are 
generated more by irregular propagation of the rupture than 
by a complex source time function. The analysis of the 
spectrum in the near-field shows a variable slope at 
high-frequencies with almost the same corner frequency for 
uniform strength, line barrier, block barrier, or random 
barrier. Far-field spectra also show a variable slope at 


high-frequencies. 


Ani ‘ne 


‘tobe ra Kats 


HQ ERTS = 
F2i) Ca arty, | 
Sec ye + RY go 35a ee Ly he aay 
a Os Sie ev outa: : iti, ss ee Hi i 
‘aie ide tee ae £10 seine at i a . 
rivseen.,' Joey Gay . sh Ce AeA ee 2a aria gor an enea - 
fal ied 59 ‘sith ea vee hs od Ba 
Mag J ¢ aps 5 ae aH Comey Sat iets wer 7 
Ad ET Ca aedR Mig ic an aiares rane BELEK De é 
i A SGP USN ‘ See ie 
vo eit ies Bt Sones inabaae | iiae | giunge: inane 
‘ * 
Seco ombe in ee a. 7 ) By Ney Aan dat 
iN stn: aes pra e5%. ie’ | 
Pe:2 a6) Sagi ee 4 hae eae baat ania 
2 sD Rete ens i tan hike tit | ae Rie law ee is 7 z 
Che? SH sores) ei bia er cone Ag zits cannes ‘ 7) don a, 
| rval | me 
2 1e ell aete ys Bra lim » pea Si Gath Plates 4 Bie ing : Lae ‘a 


Res Say tet SH 


oa 


foie ncn a ieiguesyi yShaeme ins 
7 


= enc 3 alex aes et) > oem fan fi oo, PI RYOR ke tame 
Je*sgole sidpi ria% cy empe & ears des ody at med 
469 fSristisend rarites Bina sis “peeince | an ee Sey eas 
mekrnss to 73 rit Prose Bice ala detente saath 
8 OSA 2) Seeks ay, % wtsife oe : bs stow #3 relia alt 
i, | Ms ~ elie st ha 
: iT io 


20 


1.3 Empirical Relations in Seismology 

It is an implicit assumption in seismology that 
earthquakes which occur in the same region will have similar 
source properties. This is true in spite of the fact that 
there can be a large range of magnitudes in a particular 
region. This similarity assumption (Aki 1967) implies that 
relationships might exist between various aspects of the 
source in the population of all earthquakes an a region. 
These can be investigated by analysis of statistical 
distributions of source properties. Various authors have 
proposed various such relationships. Since the use of larger 
populations may allow average statements about a source 
region, and since I use this device in the work that 
follows, I now review these relationships. 

Originally (Richter 1935) the magnitude of an 
earthquake was empirically related to the maximum amplitude 
measured in a standard seismograph, the epicentral distance 
and depth of a local earthquake in California. Evidently 
this magnitude was related to the energy radiated from the 
earthquake source. The first attempts to relate local 
magnitude and energy were given by Gutenberg and Richter 
(1942, 1956), they also proposed what they called a unified 
magnitude mb, that is derived from body waves recorded at 
teleseismic distances. 

The usually accepted relationship between local 
magnitude and energy is log(Es)=1.5(Ms)+11.8, where 


according to Richter (1956), Ms (surface wave magnitude) is 
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an approximation to local magnitude for shallow earthquakes. 
(An excellent treatment of the definition of earthquake 
magnitude is given by Pilant 1979 page 385). Another 


independent relation given by Bath (1958) is 
log(Es)=1.44(Ms)+12.24., 


Equations that have been used in energy calculations 
are those obtained by Starr (1928 equation 18) and Knopoff 
(1958 equation 33). Tsuboi (1956), assuming that shallow 
earthquakes have a certain voluminal upper limit to support 
concentrations of stresses, calculated a relation between 
energy and fault area (cm?) given by E=600xA'® (ergs). 
Pilant (1979), using simple physical considerations, 
calculated the energy necessary for a crack fo. return, tots 
original position before rupture, and obtained Es=(M,/ u )x 

ng, where n is the seismic efficiency no is the apparent 
stress, and u is the shear modulus. 

A theoretical relation between energy and source 
dimension was calculated by Gibowics (1975), using Randall 
(1972) expression for the seismic energy Es=(1/2) Ao M./y 
and Keilis Borok (1959) expression for the seismic moment, 


Mo2( 16/7) Nov ct pethat 1S given by 


log(Es)=2log(M,)-3log(r)-27.14 (cgs units). 
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With new data and improvement in the understanding of 
earthquake sources, new empirical relationships between 
source parameters (seismic moment, source dimension, stress 
drop) and local magnitude, surface wave magnitude, body wave 
magnitude and a new earthquake magnitude (Mw) proposed by 
Kanamori (1977), have been developed. 

The first attempt to relate seismic moment to 
earthquake magnitude was made by Brune (1968). He also gave 
a method for calculating the average slip along faults. Wyss 
and Brune, 1968, analysed the surface waves of thirteen 
earthquakes in California with local magnitudes ranging from 
SeeutO.5.5,.0hney found the relationship: logM.=1.4M,7+17. 

Thatcher and Hanks (1973) systematically studied source 
characteristics of Southern California earthquakes with 
local magnitudes in the range between 2 and 7. They found 
empirical and theoretical relationships between local 
magnitude and source parameters. The relationship of seismic 
moment and local magnitude (logM,=1.5M,+16) was similar to 


that found by Wyss and Brune (1968). They also showed: 


Mo=LOG(Mo)—2/ 3 log (r) )=17.8, 
log(M,)=2(M,)+14.2-logAo , 


Tog(2rj=e(2/39M,42.9-(2/3)10g Ao. 5 (cgs-units) 


for earthquakes with magnitudes from 3 to 7. Here, M, 
(seismic moment) is in dynes centimeters, r (source 


dimension) is in kilometers and Ag is stress drop in bars. 
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Source dimensions of these events lie between 0.6 to 25 km 
and stress drops between 0.3 and 200 bars. There is some 
scatter in their data, therefore, these relationships should 
be used with caution in other tectonic settings. They are 
useful in the sense that they can be compared with similar 
empirical relationships from other tectonic settings. 
Randall (1973) combined the spectra of the far-field 
displacement and the response of the Wood-Anderson 
seismograph, in order to calculate a set of relations 
between source parameters and local magnitude. Gibowics 
(1975) following a similar approach, plotted different 
theoretical curves of local magnitude against seismic moment 
as a function of source radius and stress drop. Knowing a 
given event with M, and M, it is possible to estimate the 
source radius and stress drop from those curves. As the 
source radius decreases so does the seismic moment. From the 
graphical relations of Randall (1973) Gibowics calculated 
logM,=M,+17.16, for events with source dimensions less than 
500 meters. This relation is similar to that calculated by 
Rebollar et al (1982a), using Sato-Hirasawa model, for 


Similar source dimensions. Bakun and Bufe (1975) found 


logM,=(1.52+0.05)M,+(16.2+0.1) 
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Kanamori and Anderson (1975) found theoretical 
relations between seismic moment and surface wave magnitude 
using the Haskell model and assuming similarity conditions, 
1. e. the ratio of length to the width of the fault (aspect 
ratio), the ratio of average dislocation to the length of 
the fault (constant strain drop), and the ratio of the 
product of risetime and rupture velocity to fault length 
(constant effective stress or dynamic similarity). Those 
relationships are: Ms=logM, for small earthquakes and short 
risetimes, Ms=(2/3)logM, for earthquakes with magnitudes 
between 6 and 8, and Ms=(1/3)logM, for great earthquakes. 
Further similarity conditions have been investigated by 
Ben-Menahem (1976, 1977), and implicitly by Iida and Aki 
(1972), who included the Mach number, i. e., the ratio of 
rupture velocity to phase velocity. They also obtain a 


relation between seismic moment and source area given by 
oy, 
logM,=(3/2)logSt+log(16A0 /77~) 


This is a linear relationship for a constant stress 
drop. A plot of seismic moment against fault area for 
published data indicates a constant stress drop between 10 
and 100 bars for inter- and intra- plate earthquakes. 
However, inter-plate earthquakes seems to have stress drops 
around 30 bars and intra-plate plate events around 100 bars. 

Geller §(1976))"assuming a Stress drop of 50 bars 
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logM,=Ms+18.89 for 6.762Ms 
logM,=(3/2)Ms+15.51 for 8.12>Ms2 6.76 
logM,=3Ms+3.33 for 8.222Ms28.12 
Ms=8.22 for 1logM,228 


However Sato (1979) argues that such detailed relations 
cannot hold due to the uncertainity of the data. Ohnaka 
(1978) implicitly assumed the constant strain drop 
Similarity and calculated a linear relation between 
log(M,D/L) and local and surface wave magnitude for 
world-wide and California earthquake data (D is average 


dislocation and L fault length); he obtained 


log(M,D/L)=(1.089+0.1)Ms+(8.82+2.05) Ms25.5 
Log (MyD/L)= (158820 2089M,+(9..0721249). for2. OSMis6.8 


In fact M,D/L is the seismic energy. He defines the 
"total force drop" as the product of stress drop times the 


faultareas( AF.) -vanducalculates 
Tog. ME va (1-0740 503 )Met+(120812464) 
Following a different approach Sato (1979), 


theoretically re-examined the empirical relationships 


obtained by other investigators, and suggested 
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logS (km? )=M-4.07 
logD(cm)=0.5M-1.4 


logM, (dyne-cm) =1.5lo0g$+27.5 


Where M is the earthquake magnitude, S fault area, D 
the mean dislocation and M, the seismic moment. These 
relations assume a constant strain drop (dynamic 
Similarity), proportionality of seismic moment to seismic 
energy (constant apparent stress), and the relation 
logT,=0.5M where T, the predominant period of particle 
velocity. He compared results of Utsu and Seki (1954), 
Berckhemer (1962), Bath and Duda (1964), Chinnery (1969), 
Kanamori and Anderson (1975), and Ohnaka (1978). He argued 
that the discrepancy in the coefficients in these relations 
could be due to the use of different data and uncertainities 
in the evaluation of source parameters and magnitudes. 

Kanamori (1977), considering the saturation problem of 
earthquake magnitude for great earthquakes with fault 
lengths greater than 100 km, proposed a new magnitude scale 
as a function of the strain energy drop (difference in 
Strain energy before and after an earthquake). This scale 
does not suffer saturation and gives a minimum estimate of 
the strain energy for a partial stress drop. The strain 
energy is calculated via the seismic moment, i. e. E,=( 4d /2 

H )M,, and the magnitude is calculated using the Gutenberg 
and Richter energy relation given by logE,=1.5Mw+11.8. Singh 


and Havskov (1980) calculated relationships between Mw and 
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seismic moment, according to whether they are intraplate or 


interplate events. They suggest 


Mw=(2/3)10gM,-10.73 for interplate events 
Mw=(2/3)10gM,-10.46 for intraplate events 
Mw= (2/3) 1logM,-10.57 if all events are grouped together 


These magnitudes are close to Ms for earthquakes with 
fault lengths less than 100 km. 

In the next 3 chapters I examine digital data detected, 
mainly from one single station at distances of 515 km, 180 
km, and above the source 6 km in the light of these proposed 
empirical/theoretical relations. Obviously data from a 
Single station can give biased results, but studies from 
ideal networks are generally rare for economic reasons. In 
the next chapters I argue that it is possible to have good 
results from single station data, however, more stations are 
always desirable (see also for example Bakun et al 1976, and 


Pearson 1981 for other studies using a single station). 
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2. Spectra of some Oaxaca Earthquake Aftershocks from RESMAC 
This chapter is based on a published paper (Rebollar and 
Nyland 1980). 

The Oaxaca earthquake of November 29, 1978, (Ms=7.8, 
16.07N, 96.48W) yielded an opportunity for detailed study of 
an earthquake in a subduction zone. This is the first 
earthquake in Mexico recorded with a great variety of 
instruments in the near and far field both before and after 
the main event. This event had been forecast in a gap 
(Ohtake et al 1977), and later trapped (Ponce et al 1978). 

Although the network that trapped the event (Ponce et 
al 1980) consisted of smoked paper recorders, five 
telemetering stations a combination of RESMAC (Lomnnitz and 
Gil 1976) and SISMEX (Prince et al 1973) were recording 
digitally at a distance of about 500 km from the Hee nee 
area. Since November 15, 1978, the recording system for this 
group of telemetering seismic stations in Mexico had been in 
a final testing stage. 

RESMAC stands for "Red Sismologica Mexicana de Apertura 
Continental". The RESMAC system was planned to cover all of 
Mexico. All the events are sent by microwave signals to 
IIMAS at "Universidad Nacional Autonoma de Mexico (UNAM)". 
The main task of the system will be to provide quick access 
to data for earthquakes in Mexico to Mexican and other 
investigators it witl. be fUnctioning in “event detect ron 
mode. The RESMAC system also acquires data from the SISMEX 


network operated by "Instituto de Ingenieria” (II). SISMEX 
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(Prince et al, 1973) supplied data from three of their 
Stations, santa Rita (lIGjee Tonantzintla’ (IIT) 2 and El) Pino 
GULP). 

Two days after the occurrence of the Oaxaca earthquake 
seven portable digital stations were deployed by "Centro de 
Investigacion Cientifica y Educacion Superior de Ensenada" 
(CICESE) and the University of California at San Diego 
(UCSD). This was a joint project between CICESE and UCSD. 
This network was operated from December first of 1978 to 
late April of 1979 (Munguia et al 1979). Several hundred 
events were recorded. Those aftershocks are being analysed 
by L. Munguia at UCSD as part of his PhD thesis at UCSD 
(personal communication). 

The data acquired for the Oaxaca Earthquake of 1978 
illustrates quite well the difficulties inherent in the 
collection of seismological data. Clearly the data has 
deficiencies, but just as clearly it is the best available 
data, particularly for the first 3 days of the aftershock 
Sequence when there was only a very small operation 
monitoring aftershocks in the aftershock zone. An 
understanding of the behaviour of the aftershocks of this 
very large earthquake (in fact the largest world wide of 
1978) is crucial to untangling the physical properties of 
fault zones. 

S-wave spectra of the main event and 21 aftershocks 
have been calculated. Due to the bandwidth of the RESMAC 


system, the corner frequency of the shear wave spectrum of 
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the main event is not observed. I do see two slopes at high 
frequencies (-2.0 and -2.8) intersecting at a frequency of 
1.4 Hz. Comparison of the spectra of aftershocks with those 
obtained from the nearby portable digital stations shows a 
shift of corner frequency to the lower frequencies, 
typically of the order of 0.3 hz and in general, a smaller 
value, typically 4.3 smaller, for the seismic moment. The 
corner frequencies of the P-wave spectra are outside the 
observed frequency band. A magnitude-moment relationship 
given by logM,=1.3M,+16.6 was found. This relationship is 
the same calculated in chapter 3 for events of the Rocky 
Mountain House earthquake swarm and illustrates a similarity 
between inter and intra plate earthquake generating 


processes. 


2.1 Tectonic Background of the Oaxaca Earthquake 

Mexico has one of the most active plate boundaries in 
the world, the Middle America Trench. This plate boundary is 
formed by the subduction of the Cocos Plate below the 
central Mexican part of the American Plate. Subduction of 
the Cocos plate causes shallow and deep (no more than 400 
km) seismic events in central Mexico. Shallow events are 
mainly the result of relative movement between oceanic and 
continental crust or of internal deformation within the 
plates. Intermediate (between 100 and 300 km) and deep 


earthquakes (recorded elsewhere between 300 and 750 km) 
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occur in the interior of plates immersed within the 
asthenosphere and deeper mantle. Abundant literature about 
the tectonics of the Middle America trench can be found 
elsewhere (see for example Cox 1973). 

The state of Oaxaca lies in the morphotectonic province 
of Sierra Madre del Sur (figure 1). Most of this province 
consists of schists and gneisses of uncertain, but probably 
Precambrian or Paleozoic, age. However, the highest mountain 
ranges expose middle Tertiary volcanics, whereas the 
lowlands are underlain by granites (Guzman and Zoltan de 


Cserna 1963). 


2.2 The Oaxaca Data Set 

The spectra and source parameters of the main shock and 
21 aftershocks were obtained from the data recorded at the 
Cerrillo station (CRX). The distance from Cerrillo to the 
epicentre of the main shock is approximately 515 km. For 
studies in Oaxaca , the Tonantzintla (IIT) SISMEX station is 
one of the best stations of the system. Unfortunately 
calibration information for the SISMEX stations for 1978 and 
1979 is not available to us. Three stations of SISMEX (IIC, 
IIT, IIP) and two stations of RESMAC (CRX, MEX) recorded the 
main shock and all subsequent aftershocks (with local 
magnitudes greater than 3) of the Oaxaca event. The main 
shock saturated must of the stations, however, Mexico and 


Cerrillo were exceptions. Mexico was moderately saturated 
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and is extremely noisy. Cerrillo is quiet and suffered no 
Saturation. Acapulco station was installed 30 hours after 
the main shock. 

Due to the instrument response of the RESMAC system 
(figure 2), Source parameters for larger events are obscure. 
However, it was possible to get the slope at high 
frequencies and source parameters for events of magnitude 
between 3.6 and 5.0. The source parameters at the CRX 
station were compared with the spectra of some events from 
the nearby CICESE-UCSD portable stations (Munguia et al, 


1979). 


2.3 Instrumentation 

The response curve of a typical RESMAC station is shown 
in figure 2 with the magnification curve for the CICESE-UCSD 
stations. A RESMAC station consists of a Mark seismometer 
(with period of one second and damping 0.7), an amplifier, 
two filters in cascade (period of 0.1 sec and damping 0.7), 
and an analog-to-digital converter (A-D). A PDP11 
mini-computer used for data recording. The maximum output of 
the analog-to-digital comverter (A-D) is #131072 “counts”, 
which corresponds to an input of +10 volts. The maximum 
velocity magnification is 162 counts/(micron/sec). The 
output of the A-D is recorded 36 times/sec giving a Nyquist 


frequency of 18 Hz. 
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The SISMEX system transmits data by radio link to 
Instituto de Ingenieria, where the recorded events can be 
digitized. While SISMEX records continuosly on analog 
magnetic tape, RESMAC does real time event detection on its 
input chanels. 

The configuration of SISMEX and RESMAC stations in 
November 1978 made it difficult to get reliable epicentre 
locations in the Middle America trench. The seismograms 
(figure 3) are very complicated making it difficult to 
recognize different phases. 

During fifty two hours after the main event the RESMAC 
system was triggered 151 times. During this interval of time 
fifty five events (appendix 1) were recorded but only nine 
(table 1) were clear enough at Cerrillo station for spectral 


analysis (table 2). 


2.4 Data Analysis 

From the travel time curves for a continental structure 
at normal depth (33 km), the phases that would appear in the 
seismogram at a distance of 510 km are: the refracted wave 
Pn, the direct wave Pg, the reflection PgPg (very close to 
Pg), the conversion SgPg, the refracted Sn wave, the direct 
Sg wave and finally the reflection SgSg. For this distance 
the surface waves do not appear very soon after the arrival 
of thesS-wave Consequently |; lt as unlikely that the energy 


spectrum is contaminated by surface waves. Instead the 
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Table 1.... Source parameters derived from RESMAC data 
recorded at Cerrillo. 
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Table 2.... Events recorded at RESMAC and SISMEX stations 


mie as bye 


Sve 
t 
G2 
Teneo Sb 
BMAD | ' ¥ } / 
Fy Pret LAN, 
Ces 4 
ay 
ies 
bet 
coe 
re ieee 
ud we vm ye a 
Rakes Peo hea 
hy 'ya ae 2a | 
a a oe Mala ‘ 
hie Ch eee o 
2) Beno ree. eas Sat ae ae 
CA pall e ae 10 Sy hhhp ; 
Ae TAT eee pect 
fiwgae' Byes! ee tenn Cras 
emay @  oiemas, we 
WL BIE Atta a Sot ee 


‘Sasi pars ne a ene 
D Hike hdl ied . 
eB Sct tw Hece aor 
) | eo shea| 

natn fl aes 
arte oe 


Bsteesats > +: ait 
‘cea, lygaeel 
“ Agee Me Mi ; 


‘sae i 


hats 124 One stabil ft ee A ee 
‘ aera Thay a ih, ba 


; i w x VAY iP an 


ik 


ge 


5 anne. SS) ate 


' a 
ene wg 


Ae ae 

be iy 
a) } 
ore eae 


as, eee “Ate 


$Y shi ef ev i0 
Se.f5 7S SUG 
a6. 4€ 68, ee ‘05 
ge.60) va a Hutte. 


se BE. Ad, gt UF 
el. Na @6, BT, ree v' aK 


bb ae tt at. tdaven . 
0, NG ua 04 
@h.3%,, oe Owen x 


oA ry 
ah 
ry 


Naty 


ei) 


fifty two hours after the Oaxaca event. 


Spectra will have energy from the phase SgSg and the 
Scattering due to reflection near the station. 

Spectra of P and S waves recorded at CRX were corrected 
for instrumental response and seismic attenuation. The 
epicentral distance was estimated from the average Sg-Pg 
times assuming a Vp velocity of 6.5 km/sec (based on 
unpublished refraction results at Pinotepa Nacional). This 
distance has an error of +35 km. The uncertainty in the 
distance yields an error of one order of magnitude in the 
moment. Several values of Q were tested to gain more insight 
into the value of this parameter in the region. Since very 
little is known about Q in this area, anelastic attenuation 
waS approximated with a Q of 500 for S-waves and a Q of 1242 
for P-waves (as suggested by Anderson et al 1978). From the 
analysis of the spectra, the reliable bandwidth of 
frequencies was determined to be between 0.2 to 7.0 Hz. 
Consequently, the earthquake signal was band-pass filtered 
outside this range of frequencies. 

The different time windows that were used in 
calculating the shear wave spectra (figure 4) show some 
increase in the spectral amplitude with an increase of the 
sample length. All spectral calculations were made with 2000 
Gata points, (55.6 seconds) of data. This implies the 


possible presence of scattered P and S energy. Reyes et al 
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Figure) 4. ...eVvariation of the spectra with increasing 
sample length. 
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(1979) used hypocentral determinations of aftershocks to 
establish 2 trends in the distribution of aftershocks. One 
part of the aftershock zone dips at an angle of 20° while 
another part suggests a reverse fault zone dipping at 70° 
toward the trench. There are a large number of aftershocks 
on the first portion, on the second activity is scarce but 
is supported by a composite focal mechanism. 

In order to compare RESMAC spectra with CICESE-UCSD 
data (figure 5), the spectra must be corrected for radiation 
pattern and free surface reflection of SV waves. Two of the 
events used in our comparison to CICESE-UCSD results are 
associated with the zone that dips at 20°. This admittedly 
weak result is the basis for assuming that all events 
analyzed here slipped on fault planes parallel the main 
thrust zone in a down dip direction. This corresponds to a 
dip of 20 degrees and an azimuth of 60 degrees. Nuttli, 
(1961), calculated the amplitude of S-waves at the free 
surface for different angles of incidence. Using Nuttli's 
formulas the ratio of the displacement amplitude at the free 
surface to the amplitude of the incident wave is 0.82. The 
previously mentioned unpublished refraction results (Singh, 
personal communication) were used to arrive at this number. 
An average depth of 35 km was assumed. At 60 degrees azimuth 
the average radiation pattern correction for Sv waves 
(assuming a rupture velocity of 0.9 Vs) is 0.84 (Madariaga 
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Figure 5.... Two examples of spectra for events recorded 
both by CICESE-UCSD and RESMAC stations. 
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Spectral amplitudes, corner frequency and slope at high 
frequencies (table 1) decrease for different values of Q 
(table 3). The corner frequency also tends to decrease but 
in an erratic way. 

The response characteristics of the RESMAC system were 
such that it was only possible to get the slopes at high 
frequencies of the main shock, (-2.0 and -2.8 (figure 6)). 
The intersection for these slopes is at 1.4 Hz. The first 
aftershock of magnitude 5.6 was recorded nearly 13 minutes 
after the main event, a corner frequency at 0.7 Hz was 
found, although this means a small source dimension for a 
relatively large event. The second large aftershock of 
magnitude 4.8 was recorded approximately 74 minutes later, 
this spectrum does not show any corner frequency at all in 
contrast with the first aftershock, instead it is possible 


to see a slope of -1.3 at high frequencies (figure 7). 


2.5 Source Parameters 

The constants used in this analysis of source 
Parameters (Brune 1970, 1971) are: S wave velocity=3.75 
im secpmdensity p=2.35.9r/cm>. Table 2 shows the results (of 
these calculations. 

Two (figure 7) of the 21 aftershocks spectra show 
poorly defined corner frequencies. I compared (table 4) 
source parameters of the spectra of some events registered 


at both RESMAC and CICESE-UCSD stations. In general there is 
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Table 3.... RESMAC spectral parameters for diferent values 
of Q compared with the results from the CICESE-UCSD nearby 


stations. 
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52.49 
0.40 
83.63 
sla zfal 


9 = slope at high frequencies 


Table 4.... RESMAC source parameters for Q=500 compared 
CICESE-UCSD results. 
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Figure 6.... Spectra of the main shock and the first two 
aftershocks recorded at RESMAC 
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not good agreement between corner frequencies (the largest 
difference is 0.8 Hz). CICESE-UCSD corner frequencies are 
larger than the RESMAC corner frequencies. This difference 
results in a discrepancy of 1.7 km in the source radius and 
a large discrepancy in the stress drop. This bias could be 
due to scattering in the crust (Dahlen, 1974), but it seems 
more likely to be evidence of variations in the focal 
mechanisms. 

Unfortunately only CRK was sufficiently quiet for 
Spectral analysis. It is evident from the seismograms and 
the time windows which have been used, that the RESMAC 
spectra contain scattered P and S waves. Consequently, 
scattering could increase the estimated fault size. 

RESMAC moments are smaller, sometimes by a factor of 
close to 5, than CICESE-UCSD moments. The disagreement could 
be due to: 

1) Differences in bandwidth between CRX and the 
CICESE-UCSD stations. 

2) The bandwidth of CRX. The station was not originally 
planned to provide spectral information at teleseismic 
distances. 

3) Uncertainties in the radiation pattern. I have only 
one station. 

4) Energy losses due to conversion of the SV signal and 


due to scattering in the largely unknown structure between 


Oaxaca ‘and Mexicc. 
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Consequently calculations of the moment give a possible 
lower limit. The spectra of P-waves were calculated and do 
not show any corner frequency. This could be due to the 


lower amplitude of the P-waves. 


2.6 Seismic Moment and Local Magnitude Relationship of the 
Oaxaca Aftershocks 

Using twenty aftershocks of the Oaxaca event I 
calculated an empirical relationship between seismic moment 
(M,) and local magnitude (M,). This relationship was 
compared with a similar relationship for events from the 
Rocky Mountain House earthquake swarm (See chapter three). 
Even though these seismic sequences were generated in 
different tectonic enviroments (the Oaxaca aftershocks are 
interplate events, and the Rocky Mountain House earthquake 
Swarm is an intraplate sequence). Both sets of data seems to 
fit the relationship given by logM,=1.3M,+16.6 for 
Magnitudes between 1.5 and 5.5 (figure 8). A more expanded 


discussion of this relationship is given in chapter three. 


2.7 Conclusions 

Source parameters of 21 aftershocks of a large 
earthquake on a subduction zone were obtained. Corner 
frequencies lie near one Hz, showing smaller values than 
those obtained with the nearby stations. This could be due 


to scattering along the path length and near the station. 
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—-~— SATO-HIRASAWA (1973) logM.=M,+17.6 
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Figure 8.... Comparison of the Oaxaca seismic moments with 
the seismic moments of the Rocky Mountain House earthquake 
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Seismic moments for the observed events with magnitudes from 
5. Ono. oranges from s.4 5x02 2tdyne-cmito 343x104 *Mdyne-cm, 
with source dimensions from 1.2 to 2.7 km and a maximum 
Stres drop of 4.7 bars. Consequently source dimension and 
Stress drop are lower limits. P-wave spectra do not show any 
corner frequency, being very sensitive to scattering and 
noise level. Seismic moment is the more reliable parameter. 
The magnitude seismic moment relationship for the 21 events 
analysed fits the relationship given by logM,=1.3M,+16.6. 

In this chapter I calculated the spectra of the Oaxaca 
earthquake aftershocks at about 515 km from the source. I 
can consider this a “large” epicentral distance. In the next 
chapter I analyse the refracted phase Sn and the spectra of 
the events from the Rocky Mountain House earthquake swarm 
detected at Edmonton, those events are approximately 180 km 
from Edmonton. Therefore, I can consider that spectral 
analysis at an "intermediate” distance from the earthquake 
source. The data in the next chapter will also permit 
comparison of processes on a subduction zone with the poorly 
understood processes that give rise to intraplate 


seismicity. 


‘at 


a mI 
7 age f- 
ate =: 


i 
bal: 23' 35 art oad sich | eed Ate ;. oe 


ney eat oath atnek Rented neal {e . deaengory io , waa 


hy san 


Bo Aue = vt 


he, - a bE ‘ 
a er dl i att 
: sesmeit: | 

i A ity. 

‘i Re. 


Rr: Ost Meats id # ry ‘% Pera sk ite a 
ay. Cena ate ioe oe ipedipd te 


im yell ks bets 


ORFLE OR 1 Cie en oP 444 oie i vi halts 


i, ee 


¥ 
td 
yr 


5 tem Seti yey alate, ve 
Me Ac alee | 
Aine Nea ate ee gyi Pacey 


t 
ye q 


PIT MSS: See uebt Het 


Sy) MOS MAT at a6. ai | ee 


ry uA 


yj patties onal et Ew 1 BAS oe Bas, eek: 
2 
eee ens: on paint reeoimey seabuoee, Roe : ' 


3. Focal Depths and source Parameters of the Rocky Mountain 
House Earthquake Swarm from the digital data at Edmonton 
It is often assumed in earthquake seismology that similarity 
exists between large and small earthquakes (Aki 1967). This 

is the motivation for comparing a sequence of small 
earthquakes in an intraplate environment with the 
aftershocks of a large earthquake on a plate boundary, such 
as the Middle American Trench near Oaxaca. 

In this chapter I make an analysis of the Rocky 
Mountain House earthquake Swarm, such an intraplate 
Sequence, from records taken at Edmonton (EDM) and Suffield 
(SES). From January 1976 to February 1980 one or both of 
these stations detected 220 events with magnitudes less than 
4, Some of these events show well defined Sn, Sg and Pg 
phases and small variation in the difference of Sg-Sn and 
Sg-Pg at EDM. Analysis of the theoretical travel times using 
a structure determined for central Alberta yields an average 
focal depth of 20+5 km and an average epicentral distance of 
175+5 km SW of Edmonton for 40 of these events. Since Sn was 
not clear on the remainder, it was not possible to get focal 
depths for all the events. 

Seismic moments of 80 events with local magnitudes from 
1.6 to 3.5 were found to be in the range of 6.0x10'* to 
7.9x10?° dyne-cm. A relationship between local magnitude and 
seismic moment was logM,=1.3M,+16.6; Similar to one 
determined in California. A theoretical realtionship using 


Sato-Hirasawa (1973) model was calculated with different 
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rupture velocities and source dimensions: It was found that 
it is strongly dependent on these parameters. Using a 
rupture velocity of .9 times the shear velocity and a source 
dimension of 1 km the model give logM,=M,+17.6 The 
Sato-Hirasawa model fit well our data. The Gibowics (1975) 
theoretical relationship given by logM,=M,+17.1 also fits 
Our data. Source radii, where they could be determined, were 
between 500 to 60 meters and stress drops were from 0.01 to 
Peo oars. 

The depth of focus, the low stress drops, and the 
Statistical similarity to other natural earthquakes 
sequences suggest that at least part of this swarm is of 
Natural origin. 

. This chapter is based on a paper submitted to Canadian 


Journal of Earth Sciences (Rebollar et al 1982) 


3.1 Introduction 

Destructive earthquakes are not common in Canada, 
although large events occur off the coast of British 
Columbia. Of all the seismic activity detected in Canada 
approximately two-thirds occurs iG Western Canada (figure 
9), mainly in British Columbia, Yukon, and the North West 
Territories (Whitham et al 1970). 

The Rocky Mountains, one part of the North American 
Cordillera, and the Alberta plains show some seismic 


activity. Seismic events occurring on the Alberta plains are 
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Figure 9.... Seismicity of Western Canada. 
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usually small (Local magnitude M,<4) and thus are not often 
detected on more than a few seismic stations. Some may be 
triggered by secondary and tertiary recovery methods during 
hydrocarbon extraction. A somewhat higher seismicity level 
exists in the foothills of the Rocky mountains. This 
Seismicity can be observed at Edmonton (EDM), Suffield (SES) 
and Penticton (PNT). The Rocky Mountain House earthquake 
swarm (figure 10) is a good example. 

Since January 1976 EDM has been recording signals from 
a seismic swarm near Rocky Mountain House, appoximately 185 
km SW of Edmonton. Due to the small magnitude of these 
events they have been seen mainly at EDM. Larger events (M, 
near 3) have also been seen at Suffield station (SES) which 
is about 350 km from Rocky Mountain House. Therefore, this 
study is mainly concerned with the analysis of the Rocky 
Mountain House earthquake swarm detected at Edmonton and 


Suffield. 


3.2 Analysis of the Data 

The University of Alberta seismic observatory (EDM) 
records 3 components each of short and long period data in 
digital and analog form. The analog system is a part of a 
Canadian and Worldwide Network of Standard Seismograph 
stations. In the digital system the seismometer is coupled 
With ansamplifier and a Butterworth filter. The data are 


recorded at 18 samples per second (short period) and 3 
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Braure. )0..aGeologic provinces of ‘they sSouthernuRocky 
Mountains and locations of some seismic events recorded at 
EDM, .SES;,) and PNT. 
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samples per second (long period) and is written in a 2'S 
complement binary as a 14 bit bipolar word. The maximum 
signal is 13 bits plus sign (+8192 counts) and this is 
equivalent to +10 volts level at the imput to the analog to 
digital converter. There is a WWVB bit inserted at bit 
position 15 (McCloughan and Kanasewich 1974). The response 
curve of the curve of the system (see Aki and Richards 
chapter 10 for details on magnifications curves) is given in 
Cf£iqures 11). 

About 5 events per month apparently originating in 
South West Alberta have been detected on analog records at 
EDM (figure 12). There were apparently random bursts of up 
to 25 such events. This is unusually high activity compared 
to the local activity detected at EDM since 1970 (figure 
AS ite 

Some of the events of the swarm are simple, with well 
defined Pg, Sn, and Sg phases. However, some of the large 
events are complicated showing Similar amplitudes of P and § 
waves. This could indicate variations in the source 
mechanism (figure 14). Pn is usually barely above the noise 
level. Neither the complex nor the simple events appear on 
the long period records. This suggests relatively deep foci. 

From January 1976 to February 1980, it was possible to 
read Sg-Sn ne Sg-Pg in 41 events out of 220 (table 5). 
Reading errors of Pg and Pn are of the order of 0.05 sec 
when the onset is clear, and 0.1 sec or greater for Sg and 


Sn, when a clear phase was recorded. The average time 
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with no apparent frequency magnitude relation. 
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Figure 13.... Local events (S-P times of less than 60 sec) 
recorded at Edmonton station since 1970 
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Figure 14.... Events recorded at Edmonton from the Rocky 
Mountain House earthquake swarm, left - complex event, right 
- simple event. 
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DATE Pg-Frn Sg-Sa Sg-Pg 
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Table 5.... Pg-Pn, Sg-Sn and Sg-Pg phase differences 
recorded at Edmonton station (N North South and E East West 
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difference for Sq-Sn was 1.9+0.2 sec and for Sgq-Pg was 
21.7+0.7 sec. The consistency of this difference suggests 
that the activity is confined to a small area. The average 
Pg-Pn difference of 0.6+0.1 sec was based on 6 of the 41 
events. This group of events shows remarkably consistent 
differential travel times and as a consequence I feel 
justified in treating the data as a group. The average of 
the differential times should be representative of the group 


and it is this representative value I analyse. 


3.3 Analysis of Refracted Phases 

These average differential travel times have 
implications for the average depth of seismicity at Rocky 
Mountain House as observed at EDM. Obviously it would be 
better to have data from nearby stations to determine depths 
but equipment and technician limitations prevent this except 
for short intervals. In-order to estimate source depth I 
assume a structure derived from seismic refraction profiles, 
calculate travel time curves for Pg, Sg, Pn, and Sn as a 
function of source depth and use the average Sg-Sn, Pq-Pn, 


and Sq-Pg times to constrain the depths. Since Sg-Sn and 
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Sg-Pg were tightly constrained we can place limits on the 
depths from which the seismicity originates if our structure 
is correct. I have tested the entire range of possible 
plains structures and conclude that no reasonable structure 
permits shallow source depths for the portion of the Rocky 
Mountain House seismicity which has observable Sn at 
Edmonton. 

This conclusion is important. It seems highly unlikely 
that seismicity at depths exceeding 10 km is triggered by 
the relatively small changes associated with hydrocarbon 
recovery operations on the Strachan gas field which is near 
the swarm. The events studied at EDM are most likely not 
induced. This does not preclude the possibility that other 
activity exist and is related to hydrocarbon recovery. The 
possibility that the basement structures near the foothills 
are seismically active, albeit at a low level, also has 
interesting implications for the geodynamics of the North 
American plate. 

Several observational constraints can be placed on the 
problem. The first, and most important, is the requirement 
that Sn must arrive before Sg. This must hold no matter what 
velocity model is used and no matter what source depth is 
used. The second derives from the spread in Sg-Sn. Although 
it was .2 seconds in my data I accepted solutions in which 
deviations of up to 0.35 seconds appeared. A third 
constraint, which is also crucial, is that Pn arrives before 


Pg. Here again I accept the pessimistic view that 0.35 
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second deviation from an observed average is acceptable. 

Given these constraints I must solve an inverse problem 
which has no unique solution. I determine the range of 
velocity models, source depths and epicentral distances 
which fit our observations. In order to explore the problem 
I first derived what appeared to be a reasonable velocity 
model for the problem and then searched for a source depth 
and a distance to the activity that satisfied our average of 
Sg-Sn and Sgq-Pg. 

The crust is nearly 45 km thick in much of southern 
Alberta Clowes and Kanasewich (1970), Chandra and Cumming 
(1972), Cumming and Chandra (1975). Ganley and Cumming 
(1974) found a possible Moho from reflection record at 35 km 
near Edmonton. However, those profiles are too far from the 
Rocky Mountain House activity. The nearest reversed seismic 
refraction profile which gives reliable deep crust and 
mantle velocities is that determined by Richards and Walker 
(1959) from an approximately NS refraction profile near 
113.5°W and between 50.8°N and 51.9°N (table 6). Obviously 
this profile serves only as a guide to structure but it can 
be used as a starting point for average velocity depth 
models on the wave path from Rocky Mountain House to 
Edmonton. The seismic work that was done closer to this area 
does not provide reliable lower crust or mantle velocities. 
In order to interpret the data I have at Edmonton such 


velocities are far more crucial than the shallow velocities. 
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'Layer Thickness Vp (km/sec) Vs (km/sec) 
(km) 

2.0 31.0 Veoa 

fe} Glen Sn 

22.5 6.2 Soe! 

3. 0 Wess 4.15 

Half Space la & 4.73 


FEE ES Ss mn The velocity structure (Richards and Walker 
1959 


At a distance corresponding to that between Rocky 
Mountain House and Edmonton the Sg-Pg time does not depend 
on the depth of focus. It 1s thus possible to adjust this 
depth to satisfy observations of refracted S (Sn) which I 
obtained from digital data archived at Edmonton. It was 
however not always possible to match refracted P (Pn) as 
well. This is not surprising, for Pn is not well observed, 
even on the records I used. I found that the seismicity was 
associated with depths in excess of 10 km. Best estimates of 
the sedimentary section in this area suggest that the 
maximum depth of sediments in this area is approximately 
16000 feet or 4.8 km (Bally et al 1966, keating 1966). 

Obviously such depth calculations are model dependent. 
One way to explore the degree of this dependency is to 


determine the effect of small changes in the model. Another 
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plausible structure (table 7) also yields the results that 
the average depth of the activity is greater than 10 km. The 
predicted values of the time differences for various phases 
(table 8) for this model are such that the average of the 
observed values indicates depths of the order of 20 km and 
epicentral distances on the range 170 to 180 km. Again Pg-Pn 
does not fit well. On the travel time plots for events at 
various depths with a dipping Moho and a flat Riel 
discontinuity we note that only once the source depth 
exceeds 15 km does Sn arrive before Sg. I cannot distinguish 
the effect of Moho dip on the travel times hence I compared 
my data (figure 15) with the travel times predicted for a 
flat Moho. 

The data. at EDM does not contradict the existence of 
Shallow activity at Rocky Mountain House. Only 41 of the 220 
observed events in the digital system showed positive Sn-Sq 
times. Small Sn-Sg or negative values of this phase 
difference would lead to masking of Sn in the Sg wave train. 
Those events for which Sn cannot be picked can be explained 
as occurring at depths of less than 10 km. The sharpness of 
an observed distribution of Sn-Sg suggests a concentration 
Of activity at 20+5 km depth. 

In this kind of calculation a compromise between 
computational effort and results is required. At this stage 
I have a strong suggestion, but certainly no proof, that the 
average depth of a significant fraction of the seismicity is 


in excess of 20 km. In order to explore this further I 
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3 2.68 135 
fied 4.59 PN 
29.4 62:20 Shs Sl) 
13.0 Te 20 4.15 
Half Space 3.20 4075 


Table 7.... A modified version of the Richards and Walker 
velocity structure 


Depth Er. Ds Pg-Pn Sg-Sn Sg-Pg 
5 km 170 km -1.06 -1.71 20.90 
5 km 185 km -0.47 -0.68 22.68 
10 km 170 km -0.51 -0.75 20.92 
10 km 185 km 0.07 0M 27 22.70 
15 km 170 km 0.06 0.23 20.95 
15 km 185 km 0.64 125 22873 
20 km 170 km 0.65 27 200 
20 km 185 km 1:23 2.28 22478 


Table 8.... A selection of predicted differential travel 
times for the modified Richard and Walker model 
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Depth = 5km 
Sg — Pg = 20.90 sec 
Sg-Sn = -1.71 sec 
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Figure 15.... Theoretical travel time curves for the Alberta 
model at different depths.The best match iofpea,s. SoA ara sq 
with the theoretical travel time curves was for a source 


depth of 20 km. 
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simplified the characteristics of the velocity model and 
defined a priori reasonable bounds beyond which the solution 
of our inverse problem could not lie. I admit that 7 
Parameters are unknown to some degree. These are the average 
epicentral distance of the activity, the average depth of 
the activity, the P wave velocities of two layers and a half 
Space and the layers thicknesses. I assumed that Poisson's 
ratio was 0.28 in order to relate S wave velocities to P 
wave velocities. | 

It is possible with a computer search algorithm (ZXSRCH 
from the IMSL library) to determine N possible heptads of 
numbers constrained within a 7 dimensional rectangle such 
that the N heptads are evenly distributed in the possible 
range (table 9) of the solutions. For no particularly good 
reason, other than the desire to limit computation, I set N 
to 1000 and let the computer find those heptads which define 
solutions to the inverse problem, Only two models fit (table 
10). They both require focal depths in excess of 25 km. 

Obviously a computer search of this type does not 
rigorously prove that dramatically different solutions do 
not exist between the points tested. I do however have now a 
reasonable basis for the statement that at least some of the 
seismicity at Rocky Mountain House originates at depths 
which suggest natural rather than induced causes. I also 
state that this depth can be determined if Sn-Sg can be 


observed. 
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Table 9.... The bounds on the search for 
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the velocity structure 
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Table 10.... The only satisfactory models in a suite of 


reasonable ones 


There is also shallow seismicity at Rocky Mountain 


House. One week of observations with Sprengnether DR100 


1000 


recording system (Rebollar et al. 1981) found mainly shallow 
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activity but one event with an S-P time of 2.1 seconds. 
Unfortunately this event cannot be located, for the Earth 
Physics Branch portable analog Stations, and the Edmonton 
digital station were down. A program of recording at Rocky 
Mountain House is continuing in order to acquire data for | 
* several swarm events recorded both at Rocky Mountain House 
and EDM. This would provide final confirmation of deep 
seismicity at Rocky Mountain House and raise the very 
interesting problem of its cause. 

Even with good observations of Sn the calculation of 
depth is problematical. A major change in the structure 
between EDM and Rocky Mountain House would destroy the 
conclusions. Then the identification of the phase I call Sn 
can be doubted. Procf will consist of simultaneous | 
observation at EDM and at Rocky Mountain House of a deep 
event. I suggest this may be difficult. In 5 years there 
were 48 events on which Sn could be picked at EDM, so it may 
take some time to establish the case for seismic activity in 


the basement structures of the Rocky Mountain foothills. 


3.4 Spectral Analysis 

In order to extract good source information, high 
sampling rate high-dynamic range and a broadband system are 
desirable. However, there are some cases were digital data 
of lower quality will still produce some spectral 


information. Unfortunately the digital data for the Rocky 
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Mountain House Earthquake swarm is of the latter type. EDM 
was designed to detect teleseismic activity. It has a 
Sampling rate of 18 samples per second and a reliable band 
of frequencies between 0.5 and 7 hz. 

Our basic data are the 3 digitized components of the § 
phases of the signals identified as originating in the Rocky 
Mountain House area. These were always band pass filtered in 
this range of frequencies. By suitable rotations I isolated 
SH on the transverse component and combined the vertical and 
radial components of SV to isolate pure SV wave. 

Clowes and Kanasewich (1970), studying the attenuation 
in Southern Alberta, found a Q of 300 in the sedimentary 
layers and a Q of 1500 for the basement. Therefore I 
approximate the attenuation with a Q of 1000 in order to 
take into consideration the sediments. The signal was 
corrected for attenuation, instrument response and distance. 
I use a sample length of 20 seconds in which is included Sn, 
Sg and scattered waves nese the.station. The seismic spectra 
was calculated and smoothed with a Daniel window of 0.2 Hz 
(figure 16). Corner frequencies were difficult to recognize 
in many cases due to the narrow bandwidth. 

Source parameters can be derived from these spectra 
(Brune, 1970, 1971). Since only one station was used in this 
study, the radiation pattern was approximated as the rms 
aversae over the focal sphere for S waves (0.63). I 
calculate the correction for amplification at the free 


surface for SV waves using Nutlli's formulation (1961). The 
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Figure 16.... Spectra of typical events. The amplitude of 
radial and tangential components has been increased in order 
to plot them together. The spectra were smoothed with a 
Daniel window -DW- of .2 Hz 
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correction for free surface amplification is 1.04 for SV 
vertical component and 1.75 for SV radial component 
considering an angle of incidence of 24.5 degrees. 

The source parameters, seismic moment M,, source radius 
rand stress drop are related (Brune 1970, 1971) to the 


Spectral characteristics by 


gs 
M = 4m7pRB 2(0)/RA dyne-cm 
tee O37 B/ft | km 
Ao = Q 3719 
106PR (O)£ /10 bars 


where B=3.37 km/sec is the S wave velocity, p=2.9 g/cc is 
the density,Q = spectral amplitude at zero frequency, 
fe=corner frequency, R=epicentral distance, R = radiation 
pattern, A free surface amplification. 

Computations were made on the spectra of 78 events. 
Usually the seismic moment calculated with the SH wave 
components was higher than that calculated with SV (table 
1?) SeThis could’ be an effect of the tault orientation. ‘The 
average seismic moment of each event was used. The corner 
frequencies was identified (between 1.8 and 2.9 sec) for 34 
events. Those corner frequencies gave source radii from 500 
to 600 meters and stress drops from 0.01 to 1.5 bars (figure 


17). Obviously at this epicentral distance and with this 
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Table 11.... Seismic moment, source dimension, stress drop, 
and local magnitude of some events of the Rocky Mountain 
House earthquake swarm 
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Figure 17.... Plot of local magnitude, source radii and 
Stress drop. 
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Narrow band system I sample only a restricted number of 


events. 


3.5 Seismic Moment and Local Magnitude 

Local magnitude is a useful measure of small 
earthquakes. Small earthquakes do not generate large puceaee 
wave amplitudes at periods of 20 seconds and at teleseismic 
distances the amplitude of P waves is small, so Mb and Ms 
are hard to determine. 

Since I deal here with a relatively large collection of 
data which is individually rather uncertain, average 
behavior is probably the best representation of the spectral 
information. In addition to seismic moment local magnitude 
is relatively easy to determine and I have calculated local 
magnitudes for all events for which I have obtained seismic 
moments. Seismic moments range from 6.38x10'* (M,=1.6) to 
2.21%107,% (Mi =3.5) dyne-cm. A plot of M. against M, ‘shows 
more scatter (figure 18) than that on similar plots (Wyss 
and Brune 1968, and Thatcher and Hanks 1973). My linear 


trend is similar to theirs. The best fit to our data is 
LOGiMo= lsM, +1676 
This equation is a complicated function of source 


parameters, nevertheless, it can be interpreted in terms of 


a constant stress drop fault model (Kanamori and Anderson 
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Figure 18.... Comparison of theoretical relationships 
between M, and M, for a circular fault with our experimental 
relationship, squares are the moments calculated using 
Thatcher and Hanks (1973) relation 
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1975, Sato 1979). I used the Sato and Hirasawa (1973) model 
Pom a circular fault in order to find a theoretical 
relationship between seismic moment and local magnitude. In 
general we can represent the spectral amplitude of the 


far-field displacement in the form 


A(w) = MAB(v,u,t,$,L,or r)/47pRB° 


Where 

op = density 

R = epicentral distance 
v = P or S wave velocity 
u = rupture velocity 

T = rise time 

r = source radius 

Diz fault lenoth for a rectangular) Eault 
R = radiation pattern 
@ = azimuth 

M,= seismic moment 

w = angular frequency 

B = source function 


Local magnitude is defined as 
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8 1 
M, (0.8Hz)=logA-logA,-1.35 


where A is the recorded trace amplitude for a given 
earthquake at a given distance, and A, is that for a 
Particular earthquake selected as standard (Richter 1958), 
and -1.35 is the correction for Edmonton station. Hence, 
Substituting the spectral amplitude evaluated at a frequency 


of 0.8 Hz we have 


5 


BB 
= ————_3 Se ‘ Bs 
ML logM | Oe Bl 25.05 


If we assume that all the events of the swarm have the 
Same rupture velocity of 0.9 times the shear wave velocity, 
a shear velocity of 3.52 km/sec, a density of 2.9 gr/c?*, and 
rms average radiation. pattern of 0.63 and a source dimension 
2r of 1 km (this number comes from the spectral analysis) we 


get 
logM,=M,+17.6 


for the Sato-Hirasawa model. This model seems to fit our 
data. The intercept (17.6) of the relationship, using 
Sato-Hirasawa model, was modeled for a wide variety of 
rupture velocities, shear wave velocities, densities, source 
dimensions, radiation patterns, and azimuths, using again 


the subroutine ZSRCH (IMSL, 1979). The best fit was found to 
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82 
logM, =M,+17.1 


assuming a shear and rupture velocity of the same magnitude 
(3.5 km/sec), a density of 2.8 gr/c*, a source dimension of 
100 meters a radiation pattern of 0.5 and a azimuth of 71°. 
Therefore, this search for a best fit tends to give smaller 
source dimensions, than those observed in our spectral 
analysis. 
Gibowics (1975) calculated a similar theoretical 

relation using Randall (1973) graphical relations, assuming 
a circular fault and the asymptotic behaviour of the source 


function, given by 
logM,=M,+17.16 


for source dimensions of less than 500 meters. This is 
essentialy the same relationship that I found. 

It should be noted that the stress drops reported here 
are just that, fractions of an unknown absolute stress. The 
stress drops reported here are at best lower bounds on the 


actual values of the absolute stress. 


3.6 Conclusions 
Some of the events detected at Edmonton from the Rocky 
Mountain House Earthquake swarm appear to have a source 


depth of nearly 20 km. This conclusion is based on the match 
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of refracted phases with theoretical travel time curves for 
a model of Central Alberta. 

A portable digital seismic station, operated during 
August of 1980 at 52.26°N and 115.13°W, approximately 160 km 
SW of EDM, indicates mainly shallow activity (less than 6 
Km). Nevertheless, there is one event with S-P time greater 
than 2 seconds, indicating a distant or deep event. 
Therefore, two kinds of activity, deep activity possibly 
associated with thrust faults of the basement of South West 
Alberta and shallow activity possibly induced by secondary 
and tertiary recovery methods in the Strachan gas field may 
exist. 

Seismic moments calculated for events with local 
maqnucude from 1.6) tOus.5, Vary rom 6.s6xi0. itor 2. 2x10 2° 
dyne-cm. A relation between seismic moment and local 
magnitude was found to be logM,=1.3M,+16.6 for magnitudes 
from 1.6"to 325. This relationship 1s similar to that found 
in California by Wyss and Brune (1968) and Thatcher and 
Hanks (1972) for events with magnitudes in the range from 
3.0 to 6.0. A theoretical relationship was calculated using 
Sato-Hirasawa (1973) model. A source dimension of one km and 
a rupture velocity of 0.9 the shear velocity was used. This 
relationship seems to fit our data best, and is essentialy 
the same calculated by Gibowics (1973). Source dimension and 
stress drop do not show any linear relation to local 


magnitude; this could be due to the poor resolution of our 


system. 
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In the next chapter I analyse the Rocky Mountain House 
earthquake swarm activity recorded at "short" epicentral 
distances (less than 6 km) by a Sprengnether DR-100 digital 
intrument. This provides further information on the 
reliability of source parameter estimates and a better 


insight into shallow activity. 
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4. Source parameters from shallow events in the Rocky 
Mountain House earthquake swarm 

This chapter is based on a paper in press (Rebollar et al 
1982). I report here source parameters of the Rocky Mountain 
earthquake swarm derived from three component digital data 
recorded: on a nearby station. During 6 days in October 1980 
21 events were recorded. Focal depths for these events are 
in the range of 2+2 km. Eleven events with local magnitudes 
from 2.1 to 2.8 yielded source parameters. Corner 
frequencies of the S-wave spectra were found in the range 
6.2+0.5 HZ giving source dimensions of of 160+10 meters. The 
corresponding P-wave corner frequencies are in the range 
8.6+3 Hz. The ratio of P to S corner frequencies varies from 
0.9 to 2.1. There is a path effect between 13 and 16 Hz, 
that could have affected these ratios. The average fall-off 
over 3 components at high frequencies varies from -1.8 to. 
-2.3. High stress drops, ranging from 47 to 263 bars and 
apparent stresses from 2.5 to 23 bars, were calculated. Five 
events have remarkably similar characteristics in the 
frequency and time domain. For these events the ratio of 
minimum strain energy W., according to Kanamori (1977), and 
the energy calculated using the integration scheme of Hanks 
and Thatcher (1972) was 3.7+0.5. A theoretical value gives 
3.1. The seismic efficiency ranges from 0.2+0.04 to 
0.17+0.8. Large seismic moments for a relatively small 
magnitudes were found. Some of these Spectral 


characteristics are best explained as the result of 


85 


ts, Jo: NACL Og) aye ot 3 

LG DRM | ok oats anit i dhs remer9g 99 i ae Hid aisll 
esa) fad) pas sad ai ead > soa na ; 
938 soAeVs 
 Behosringam, Deo id iw  egeews 
oe raers ia adermnies, ‘ 
Ses. Bt.) a AD A dae B12 ge: . . 8 59mg 

; RY 

ott sevaven Ot ecbit! Bp de pha taamia feuuts enti aH Bs 


ei 
Sree ait Hi ete ap toneupett soma vena | Tia 


mort) a9 at getsnaupers tampon 2 of & iq nite se ee 
Mee ed 1, aa 


i 


ni a Ung e4 Neonved: sausage Aeeq. ‘ at Ts \ wa 
tionfien e;.vevs ant papic cw “wees bpsoet in: eva f 
oy BAT?) meat ae hmee > oil ah a6 i 


Ag 


She e2ac 1a a? Ts mox> ani “ 4 
avira bate tinbes  s2am a hee BS eke SL eed ropeandy 
etd) ta a) vat tase ted e ‘tents ) Pidesageen eat: 
Io prdes aaa sites oobi a .atwoob mts ee utp 
bas eee jromenel or eninaeas 7 28 eee amt 
# exit be seid “tan 


esiish 30 "alpine ine aa idempeons ae 


pevip wi isy sot setsey 4 bs Er ss at +8 oer 
pt, 80. CeS8 teas 2 44 ity 9 ‘she: | ps8 Se 
[leme : feviseteady & ; ns 2 “ 
iangoeha ‘oe ath | itis? wiew penne 
to tluedy ody tee Bae Peat wie esisabieroed, 


a ‘ thie ies ; i Si 
oy , f f Gel we 


displacement along a smooth fault. 
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During late September and early October of 1980, a 


local array of five analog vertical portable seismic 


Stations was operated by R. J. Wetmiller of the Division 


Seismology and Geothermal Studies of Earth Physics Branch, 


Dept. of Energy Mines and Resources. One University of 


Alberta digital 3 component portable station was deployed 


near Rocky Mountain House Alberta in order to 
seismicity of the area. This area is the most 
seismicity in the South West Alberta. Some of 
has been detected at Edmonton (Rebollar et al 


the first time that this kind of analysis has 


study the 

active zone of 
this activity 
(Swe arThis .1s 


been done in 


South West Alberta at short epicentral distances (less than 


5 km). 


Events recorded by the analog and digital instruments 


were confined to an area of approximately 8 square km near 


S202 Nand 1157) 14.W at a depth of 422 km) (figure 19) 


(Wetmiller 1981). Those events may be confined ina 


hydrocarbon bearing sedimentary layer from the upper 


Paleozoic. 


This study is mainly concerned with the evaluation of 


focal depths, and source parameters with a generally 


accepted source model (Brune 1970,1971). These source 


parameters are compared with already published data, 


theoretical models and results from laboratory experiments 


of shear crack displacements in an effort to explain their 


particular features. 
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Figure 19.... The location of the Rocky Mountain earthquake 
Swarm. The dash square indicates the seismic locations 


reported by Wetmiller 1981 an the locations of the digital 
Station. Almost all the locations were near the gas wells 
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4.1 Digital Recorder 

A three channel Sprengnether digital recording system 
model DR-100 was used in this experiment. This recording 
System was connected to a five-pole antialiasing Butterworth 
filters and to Mark product L-4C-1.0 Hz seismometers, to 
provide a single 3 component portable station. The sampling 
rate used in this experiment was 100 samples per second. 
Absolute time was synchronized with the WWVB signal. A 
response curve for the combined recorder-seismometer for 
maximum gain of 120 and 60 db is given in (figure 20) The 
Station was located at 52.23°N and 115.27°W, and operated 
from October 3 to October 9 1980 figure 19. The station was 
set up in sediments from the Cretaceous, and operated with a 
gain of 66 db due to local background noise (mainly wind 


generated and gas field activity). 


4,2 Event Locations 

In order to analyse the events we need a structural 
model. Richards and Walker (1959), from an approximately NS 
refraction profile near 113.5°W and between 50.8°N and 
51.9°N, reported two sedimentary layers above the basement. 
The uppper one is of Mesozoic age with a thickness of 2 km 
and a P velocity of 3.6 km/sec. The lower one is from the 
Paleozoic with a thickness of 1.5 km and a P velocity of 6.1 


km/sec (figure 21). 
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Figure 20.... Magnification curves of the digital station 
at 60-ands:120#db.s Themstation was? operatede at-'66r dbs 
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Pigure 21.... Travel time curves for different depths using 


Richards and Walker model 
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One of the events shows no discernible P or S waves. 
Five events are remarkable similar in all three components 
Bottomvot (figure 22) and’ (figure 23). This jis) not, the first 
time that this kind of earthquakes has been observed. Geller 
and Mueller (1980), reported four similar earthquakes in 
central California at epicentral distances from 20 to 51 km. 
This also has been observed in the Victoria earthquake swarm 
of 1978 in the Valley of Mexicali, Mexico, and in the 
aftershocks of the Oaxaca deeecunike of November 28 of 1978 
(L. Munguia personal comunication). This similarity implies 
that they are clustered in space so that they sample the 
Same path to the recording station. The maximum amplitude 
was recorded in the EW component, where there is a dominant 
Signal with a period of about 0.2 sec, that correspond to a 
wave length of 0.6 km assuming a S phase velocity of 3.0 
km/sec. | 

The S-P times of these events are in the range 1.0+0.4 
seconds. The P arrival was read in the vertical component 
and the S arrival in the horizontal. component. I calculated 
theoretical travel time curves uSing the Richards & Walker 
model (1959). A travel time curve for a source depth of two 
km gives an epicentral distance of 3 km, for a source depth 
of 4 km this correspond to an epicentral distance of 2.65 
km. Deeper source depths do not predict the observed S-P 
time, for example a travel time curve for a source depth of 
6 km predict a minimum S-P time of 1.2 seconds above the 


epicenter. Therefore these events have a source depth of 
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Figure 22.... Similar events and the spectra of the _ 
transverse T and radial component R of the S wave. PZ is the 
P-wave spectrum calculated from the vertical component, 
arrows show corner frequencies. 
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Figure 23.... Similar events and the spectra of the 
transverse T and radial component R of the S wave. Pz is the 
P-wave spectrum calculated from the vertical component, 
arrows show corner frequencies. 
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approximately 2+2 km figure 21. 

These events do not have a large content of high 
frequencies, even though the station was very near the 
epicenter. It is well known that irregular rupture 
dislocations with variable slip and complex strength on the 
fault enhance high-frequency waves (Miyatake 1980a, 1980b). 
Perhaps these events correspond to a process in which a 
barrier is being repeatedly broken and healed. In this case 
the events share the same dynamic and kinematic 
characteristics. Another possibility is a collinear shear 
crack (Rudnicki and Kanamori 1981). Geller and Mueller 
Suggest that the clustering and similarity indicates that 
these earthquakes can represent stress release at the same 
asperity. 

One event with a relatively large S-P (2.1 sec) time 
was recorded, therefore is not possible to determine the 
approximate depth or epicentral distance. I filtered the 
Signal for this event with different bandwidths and still 
observed clear P and S waves. Perhaps it is one of the deep 
ere analysed using Sn phases at Edmonton. Unfortunately 


this event was not recorded at Edmonton (figure 21-b). 


4.3 Calculation of the Spectra 


Source parameters were calculated using Brune's (1970, 
1971) model. This theory has been found to be valid in 


general both at short epicentral distances of a few 
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kilometers and at teleseismic distances. See, for example, 
Hanks and Thatcher (1972), Hanks and Wyss (1972), and Hanks 
(1981) among others in an extensive literature. 

Corrections applied before the calculation of seismic 
moment are: free surface reflection of SH or SV waves, 
radiation pattern, seismic attenuation, geometrical 
Spreading, and instrument response. In this study eleven 
reliable events were used in order to get source parameters. 
Epicentral locations, azimuths and angle of incidence 
calculated by Wetmiller (personal comunication) were used. 
These events were rotated in order to have a "pure" SH wave 
in the transverse component, and the horizontal component of 
the SV wave in in the radial term. Corrections for free 
surface reflection of the radial and vertical component of 
SV waves were calculated using Nuttli's (1961) formulas. 
Taking into consideration short epicentral distances and the 
sedimentary layers of the South West Alberta, a Q of 150 was 
used to correct for attenuation in both S and P spectra. 
However, errors in the choice of Q are not critical because 
of short epicentral distances. 

The sedimentary layers of the South West Alberta 
consist of shales, siltstones, sandstones, and carbonates 
with accumulations of hydrocarbons ranging in age from the 
Paleocene Period to the Cambrian period in the Paleozoic Era 
(Bokman 1963). Therefore, I assumed a P wave velocity Gi) .5 | 
mmiyeec at censity of 2.5 19m cc 1anGea (Porsson caciono fe0 23a, 


in order to take into consideration the low velocity body 
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waves of sediments reported by Richards and Walker (1959). 
Sample lengths of 1.0 and 1.5 seconds were used for S wave 
spectra and 0.5 seconds for P wave spectra. All the spectra 
were filtered between 1 and 45 Hz and smoothed with a Daniel 


window of 1 Hz. 


4.4 Source Parameters and Discussion of the Spectra 
Eighteen events with S-P times of the order of one 
second were analysed; of those eleven yielded source 
parameters. The spectra of these events show a significant 
path effect between 13 and 16 Hz in the transverse or 
North-South component and vertical component of the spectra. 
This effect is less prominent in the radial or East-West 
component. This path effect makes it difficult to recognize 
corner frequencies, and sometimes it suggests a spurious 
corner frequency at high frequencies, see figures 25 and 26. 
In the spectra of P waves this effect is even more 
noticeable, moving the corner frequency of four events up to 
11 Hertz (figure 24). Spectra of some of the events are 
shown in (figure 25) and (figure 26). Similar events show a 
great resemblance in the spectra (figure 22 and 23). 
Laboratory experiments have shown, that the spectra of 
waves radiated from the displacement of a homogeneous fault 
have a simple spectrum i.e. a low level amplitude an the 
fall-off at high frequencies with a few peaks decreasing in 


amplitude. However, displacement spectra from an 
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Figure 24.... P-wave spectra of the events of figures 25 
and 26. These events show a possible path effect in the 
corner frequencies; P-wave corner frequencies were usually 
higher than S-wave corner frequencies 


| ) 
‘ 


fe ee 


: 3S bas 
ene § 

: , é DAs, actin 

at ee Tits. ah ee 


q M y : ts) 7 


98 


*satztauabowoyurt yyzed atqtssod soys 

e1}z098eds aulog *‘Satouanba1rjy 19u109 moys sMOoriy *stsA[Teue styy 
UT pasn SjzUaAa 9Yy JO BAeM-S AY} JO JUaUOdWOD JY asriaasueIly 
pue ‘y Tetpea ‘Z TeOTAaA |BYQ Jo e1}0Dadg *°°"°GZz aaNnbtIg 


29°01 80 1000861 1t°90°~O LOO 086i 24S'S0'€0' 190 0861 
0001 ZH 1=M0 OSI=O 344 007 on 0°00! ZH |=MO OSI=0 944 901 on oool ZH 1=M0 OSI=O0 34¥43 901 ot 
0000°2 0008" ¢ 0002" 0000°6 000°'o 0000°0 0000°2 0009°% o00z"t 00080 00080 0000°0 o0n0'2 0008"! oooz-' 0008-0 000°-6 0000°8 
a Bat a gg 4 4 4 nt gee (ag gg lg ost (- 4 at bt ak kn ggg 
voc eee ty cag 
o o 
ee zz e- w ore 
~~ AY 
vee eo-gm Ire 
9 om 
wt-o 62 oO ere 
m m 
oo-z- ~ oer ~ 90°e 
Oo ie) 
or-z- = 99-2 >t 20-8 
u Wp vi 
oe-z-m z9-z- 1m Guz 
o2-e- \ err. eerz- 
: aee- < ; vz 1 c9-2- 
00°2- ~ ole- ov-z 
eco ‘Ape : ere 
8 3 | 
oe-2- g ' TGs are 
v 
mem oer 06-2 
99-2- ‘2- vez 
sto wee 
ov-2- ~ oez- = ou? 
ial is) 
ve-z- % zee F ose 
“ wn 
so-z-™ vem tore 
ze-4- 0-2 ez-2 
ore ; ae » Ue 
Fe ws ¥ Pate a oer Seal v vice 
00: 1- ae out as 00° 2- 
oe ws pS 
a ng: 29°8- 
o6-2- 9e°8 2 
ew ou ec'g- 
vez eereaes 
ors ze'e- ve-e- 
Reto m 
o9-2- 7 ose. = ore. 
oO Oo 
ov-2- 7 ove Ft oe s- 
S a 
oc-z- im see. ze-2- 
e2-2- eee: VA eo-2- 
alee age oa-2- 
2 gre 2 ze z 
Ash o0:8s oe 


336-42 N30 3348 907 


" 338-3 320 3348 967 


338-3 N20 3268 SCT 


18 


: ere 


ze 


tz" goles 2Sectrs 
ae J 


s 


o.e- Gf Spt caspre Geeg 


im. 


> 55 


Senet 
a a 


= 
£ 


APET GST 


tje ss 


= ios 26 


ME oT 


Ske Oy $36 


= 
gue 


To" VEG 


| 


: 
é 
% 
a 
z 
: 


‘ee - 3 


ns ae a 


= oa 
ae 
a 


i cid 5 
eas an) 
‘ot a 


stun cad ba ie, ie ; 


99 


*satqzrauabowoyut yyed atqtssod smoys 

eiqzoads auwosg *satouanbaay aaui109 moys sModiy “*stsATeue styy 
UT pasn SzuaAa ay} JO aAeM-S AY JO JUusUOdWOD | aSi1aASUeT} 
pue ‘y [Tetpea ‘Z TeOTIAAA ay JO e1}39AEdS °°°°9z7 BINbIg 


9i'ct'80' 190 0861 8¢'00'S0' 1900 0861 Se€81'€0 190 0861 
0001 ZW 1=MO OSI=0 344 901 on) 0001 Z4 1=M0 OSI=0 344 901 ot ; 
! : , : : : ) 4 4 0001 ZH Ue MU OSI=0 aNd 9 t 
eee rte Ue nln ban On Gee agree 2 pesthpll on reer n, ACSAROR CAGED SCA sores eopee een woge0 | 80090 00008 


ehh ge 


8 
PGB 
234$ 907 


° 
J3S-W3I N3O J36S 207 


* 93S-wa N30 


<« 
a |e SS SS SE RES pS gp 
e723 2-8 $93 3 3 3 sss 28 $8 2 38 
SS) Se Gt Ropes oe ROW ON GN Gh Gg 
J3S-WI N30 J3dS 907 23S-W3I N30 J3dS 907 
ae 
Cn 
——F 4 
¢ 
\ 
\ 
/ 
i) ¢€ 
$338 28 58 8 SS) AeSr se) ap i2e 5) S 
i ed - 4 SB HH RR 
J3S-wW3 N30 33dS 907 23S-WI N3O J3dS 901 
2 ped 
ye a 
\ 
\ 
€ 


oe-¢- 


338-43 N30 23dS 907 


34S 9017 


e 
* 339-49 N30 


2 
338-W9 N30 334s 907 


Sux 


e 


coe 2boc rte 
raed Fu 


Secere 


RIGS: 


Giese tee 
FPpe 


_ 
= 
a 


Feb Et se 
P-AGLG OFT 


LEE ;OE- 
us 


ape 
oe ges 


Se tae “ASLETEs TT 5" 


ARES, 
i Or 
a 


100 


inhomogeneous fault (with asperities or obstacles) have a 
different form; the low level amplitude and the peaks in the 
fall-off at high frequencies are comparable in magnitude 
with the low level trend, suggesting two corner frequencies 
(Vinogradov 1978). Therefore, these results could explain 
high peaks at "high" frequencies in our spectra. However, 
the same experiments show that shear displacement along a 
homogeneous fault (without asperities) gives a high seismic 
moment for a relatively small magnitude. This is in 
agreement with our results. Since seismic moment is a more 
reliable parameter we conclude that those effects present in 
the observed spectra are more likely due to local 
inhomogeneities. 

Corner frequencies were found in the range of 5.8 to 
6.3 Hertz. This is equivalent to source dimensions from 161 
to 148 meters. Some corner frequencies in the P wave spectra 
dbeestrong!ysatlected by pathvettects: asacan be seen wing ehe 
Spectra of figure 24. Hanks (1981) pointed out that source 
Parameters calculated at a single station are more likely to 
be affected by path effects. 

Stress drops for those 11 events were consistently high 
ranging from 47 to 263 bars. Previous work has found stress 
drops in the range of a few tenth of a bar to 100 bars for 
events with local magnitudes from -1.3 to 3.4, Wyss and 
Brune (1968), Douglas and Ryall (1972), Thatcher and Hanks 
(1972), Bakun et al (1976), Jonhson an McEvilly (1974), 


Fletcher (1980), Marion and Long (1980), and Rebollar et al 
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UrSe 1). 

Bakun et al (1976) reported stress drops of the order 
of 245 bars calculated from the spectra of P waves. Hartzell 
and Brune (1977) found stress drops in the range from 1 to 
636 bars for 61 events of the Brawley earthquake swarm with - 
magnitudes from 1 to 4.7. 

Stress drops of possible deep events of the Rocky 
Mountain House earthquake swarm are of the order of a few 
tenths of a bar (Rebollar et al 1982). Since the shallow 
events have high stress drops this seems to contradict the 
observations that the stress drop increases with depth 
(Fletcher 1980, and Hartzell and Brune 1977). However, this 
observation could indicate that shallow events might be an 
indirect consequence of the extraction of oil and gas 
whereas deep events are more likely of tectonic origin. 
Figure 19 shows the epicenter area reported by Wetmiller 
(1981). Those epicenters tend to lie near the gas wells of 
the Strachan fields. 

Source parameters of the five similar events are more 
homogeneous, reflecting almost the perfect match in the 
seismograms. A plot of local magnitude versus seismic moment 
(figure 27) shows consistently large seismic moments for 
this range of magnitudes. However, they follow the same 
slope as the relation logM.=1.3M,+16.6 (Rebollar et al 
1981). Large seismic moments for relatively small magnitudes 
have been observed in laboratory experiments of shear 


sliding along a smooth fault (Vinogradov 1978). This could 
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explain those high moments, even though laboratory 
experiments are an over-simplification of a real earthquake. 

The radiated energy E was calculated using Hanks and 
Thatcher (1972) integration scheme. They found an analytic 
expression by integrating the far-field shear displacement 
Spectra proposed by Brune (1970), assuming a complete stress 
drop (€=1) and a fall-of at high frequencies of -2. Thatcher 
and Hanks (1972), however pointed out that uncertainities in 
the fall-of at high frequencies can give a misleading 
measure of the energy. In our case the average slope in the 
three components of the shear spectrum at high frequencies 
varies from 1.8 to 2.3 (table 12). Therefore, errors in the 
calculation of the radiated energy due to the slope at high 
frequencies are small. Most of the energy of small 
earthquakes is radiated in short period waves, hence, 
because of the small epicentral distances litle energy can 
be lost. Therefore, the integration scheme can give a better 
estimate of the radiated energy. 

Energies range from 2.4x10'* ergs (M,=2.1) to 2.1x10'’ 
ergs (M,=2.8). Radiated seismic energies calculated using 
Gutenberg and Richter, and Thatcher and Hanks empirical 
relations derived for California give smaller values, as we 
Can See from (figure 28) This could mean that the 
integration scheme gives a better estimate of the radiated 
energy at short epicentral distances (less than 4 km). 

Corner frequencies of P wave spectra were consistently 


larger than those of the S wave spectra. Only one corner 
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frequency of P waves was found to be less than the 
corresponding corner frequency of the S wave spectra table 
12. P wave frequencies seem to be more affected by local 
path effects. Those path effects are easily recognized 
because they appear in the same range of frequencies 
independent of the obvious variation in the radiation 
pattern that can be recognized by inspection of the 
seismograms. Ratios of the P and S frequencies (corner 
frequency shift) were found in the range of 0.9 to 2.1. 
However, the average ratio of fp to fs in the similar events 
is 1.2+0.2, more in agrement with the general observation 
that fp/fs is greater than one, and with theoretical models 
that treat earthquakes as equidimensional faults (Brune 
1970, Madariaga 1976, Burridge 1975), or long and narrow 
faults with near-sonic or transonic rupture velocities 
(Savage 1974). Hanks (1981) suggests that the frequency 
Snacteiscan intrinsic characteristic. of ithe farstireid 
spectra of body waves depended: of source strength 
(Seismic moment), hypocenter, epicentral distance or 
recording device. 


The minimum strain energy drop W, (Kanamori 1977) is 
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Which assumes a complete stress drop (Opa Oe ors vi the 


Orowan (1960) condition is met (oc =0,). Substituting the 


2 
moment and stress drop according to Brune (1970,1971), we 


get 
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Where p is the density, 8 is the shear wave velocity, R is 
the epicentral distance, fc is the corner frequency, ©, is 
the low level amplitude of the SH spectra, and 1.26 is the 
product of the average radiation pattern and the free 
surface reflection of SH waves. The radiated seismic energy 
Es according to the integration scheme of Hanks and Thatcher 


(1972) is given by 
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Taking the ratio of W, and Bs we get Wes Uae 


If we assume that 1/3 of the seismic energy 1S contained in 


the P wave we have 


This means that even though W, is a minimum estimate of the 
energy, the integration scheme gives even lower values of 
the radiated energy. Values of W., Es, and the ratio of 


W,/ES are given in (table 13). This ratio for the similar 
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(Kanamori, 1977) (Hanks and Thaccher 1972) 
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ini ding 
Table 13.... Comparison of minimum strain energy accor 
to Kanamori 1977 and Radiated seismic energy by Hanks and 


Thatcher 
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events is 3.7+0.5, in good agreement with the theoretical 
result (W,=3.1Es) which does not include the contribution of 
the P-wave energy. However, events 1 to 6 give an average 
ratio of 6.8 with a large standard deviation of 3. 


The apparent stress is defined (Wyss 1970) as 


Where n is the seismic efficiency, o is the average shear 
Stress, vp is the shear modulus, E,is the radiated seismic 
energy, and M, is the seismic moment. Uncertainties in the 
evaluation of the apparent stress are a direct consequence 
of the uncertainties in tne evaluation of the radiated 
energy. Therefore, the seismic efficiency is one of the most 
uncertain parameters in seismology. However, Wyss (1970) 
calculated a seismic efficiency of 0.1 for deep and 
intermediate earthquakes in South America. Apparent stress 
calculated by the Wyss formula range from 2.5 to 23.2 bars 
fablew!2-sAsplotvor stress Gropaversuspapparentmseness can 
be compared with that calculated by Hartzell and Brune 
(1977). Even though there is some scatter it follow a 


Similar trend (figure 29). Those differences indicate 
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regional variations in the apparent stress and consequently 
the state of stress in the crust of the Rocky Mountain House 
area. 

The ratio of stress drop to apparent stress for similar 
events is almost constant at 7.4 with a standard deviation 
of 1.1. This ratio suggests a possible way to calculate 
absolute stress, however, in order to do this it is 
necessary to know the seismic energy and the frictional 
energy. 

Following Wyss (1970), it is possible to calculate an 


approximate upper bound for the seismic efficiency given by: 


Seismic efficiency calculated in this way for the similar 
events gives a constant value of 0.2+0.04. This reflects the 
need to pick the correct cornet frequencies in the spectra. 
Events 1 to 6 give 0.17+0.8 (table 12). These high values 
apparently suggest a high conversion of potential energy to 


elastic seismic energy for these shallow earthquakes. 


4.4.1 Further Field Work 
On September 19, 1981 I deployed, in the Strachan gas 


field, a Sprengnether DR-100 digital seismic station in 
order to detect deep seismic activity from the Rocky 


Mountain House earthquake swarm. This station is located at 
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52.23°N and 115.27°W, at the same place where it was located 
during the study of 1980. This station has been recording 
Shallow activity (S-P times of less than one second), like 
that reported at the begining of this chapter, and possible 
deep activity. 

During three months of: observation two events were 
recorded both at the nearby station and at the analog and 
digital station of Edmonton. The first event occurred on 
September 23, 1981, and was recorded at 11 hours 37 minutes 
and 50.5 seconds at Edmonton. It has an S-P time of 20.5 
seconds and does not show any refracted phase in the analog 
record. This could be due to its small magnitude (M,=1.8). 
This event recorded in the Strachan gas field station shows 
a clear P and S wave. Due to drifting problems with the 
internal clock of the Sprengnether intrument in very cold 
weather, it is not possible to read with accuracy the 
arrival time of the signal. Even though the S wave saturated 
the instrument, it was possible tormread aneS=P time, of% 1/542 
seconds, that is equivalent to a maximum epicentral distance 
of 11.5 km assuming a half space with a P wave velocity of 6 
km per second. This also could be a possible deep event. 

The second event was recorded at Edmonton on November 
25 O81 ear 14 hours, 25 Minutes, 12.9. Seconds. s inicrevent 
with a local magnitude of M,=2.8 show clear refracted phases 
in the analog system; it has a Sg-Pg of 20.8 seconds and a 
Sg-Sn of 1.2 seconds (A typical Sg-Sn in the digital records 


is 1.9 seconds). The Sprengnether station show clear S-P 
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times of 1.54 seconds, that is equivalent to an epicentral 
distance of 12.6 km. Therefore, those events with S-P times 
of 1.5 seconds recorded in the portable digital station, if 
not definitely proving the existence of deep seismicity at 
least show the posibility of events with depths of more than 
10 km, in agreement with our conclusion of chapter three. A 
definite solution of this problem requires the deployment of 
a network of three digital seismic stations for at least six 


months of continuous operation. 


4.5 Conclusions 

Events with S-P times of about one second, and source 
depths of 2+2 km, show different amplitudes of P and § 
waves, indicating dislocations with distinct fault 
orientations. Similar events have small variations in source 
Parameters, perhaps indicating a process in which a barrier 
or asperity is being repeatedly broken and healed. Stress 
drops and apparent stress were consistently high for all the 
events. Hydrocarbon recovery operations can cause 
concentration of stresses in surrounding areas. This could 
result in brittle fracture in the sedimentary rocks or the 
granitic basement by activating ancient faults, joints, or 
by rupturing along pre-existing cracks. Laboratory 
experiments reveal a lesser accumulation of strain energy, 
and consequently stresses, in a closed fracture as compared 


with an open fracture (Shamina et al 1978). Therefore, high 
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stresses could indicate a highly fracture zone, since stress 
drops represent minimum tectonic stresses related to the 
seismic events. 

Radiated seismic energy calculated using W,=(Ao/2u)M, 
(Kanamori 1977) and Hanks and Thatcher (1972) give a ratio 
of 3.1; in good agreement with the observed ratio of 3.7+0.5 
for the similar events. Larger energies were obtained for a 
given magnitude using Hanks and Thatcher (1972) method than 
those calculated using the empirical relationships of 
Gutenberg and Richter and Thatcher and Hanks for events in 
California. This could mean that the integration scheme 
gives a better estimate of the radiated energy at short 
epicentral distances (less than 4 km). Values of 0.2+0.04 
and 0.17+.8 for the seismic efficiency were found, 
suggesting a high conversion of potential energy to elastic 
energy for these shallow earthquakes. 

Corner frequencies of the S-wave spectra were found 
between 5.8 and 6.7 Hz, given source dimensions from 150 to 
169 meters. The ratio of corner frequencies of P-wave and 
S-wave spectra give a value of 1.2+0.2 for the similar 
events, as is usually observed (Hanks 1981). However, corner 
frequencies for events 1 to 6 give ratios ranging from 0.9 
to 2.1. These ratios could have been affected by local 
inhomogeneities. Evidence for such anomalies is found in our 
Spectra between 13 and 16 Hz. 

A plot of local magnitude versus seismic moment gives 


systematically large seismic moments for relatively small 
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magnitudes. However, they follow the same trend as the 
relationship found for deep events by Rebollar et al (1981). 
Laboratory experiments show that displacement along a smooth 
fault is a possible mechanism of earthquakes with large 
seismic moments and a relatively small magnitude (Vinogradov 


1978). This could explain those high moments. 
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5. Other seismicity of South West Alberta and Conclusions 
Seismic records from the Edmonton seismological observatory 
Since 1970 provide evidence for active seismic zones in the 
South West Alberta and in the Rocky Mountains. There usd 
seismic zone near Rocky Mountain House, and another near 
McNaughton Lake. An earthquake occurred in the Willmore 
Wilderness Provincial Park (The Willmore event) on October 
9, 1977. The seismic moment of this event calculated at 
Edmonton from the SV spectrum was 6.7+2x10?? dyne-cm. The 
seismic moment of the McNaughton lake earthquake of May 14, 
1978 calculated from SV waves at Edmonton, was 7.6+2x10?3 
dyne-cm. This compares reasonably with the value calculated 
from surface wave spectra (Rogers et al, 1980) at several 
stations. 

During the last 10 years there has been a denser 
distribution of seismic stations in British Columbia than in 
Alberta. Alberta had only (table 14) two seismic stations in 
continuous operation (EDM and SES), making epicentral 
locations of micro-earthquake activity in South West Alberta 
difficult. Earth Physics Branch of Canada usually reports in 
its monthly bulletin, unlocated” events detected’ in a single 
station. Those events recorded at Edmonton usually have 
local magnitudes between 1.7 and 2.5 and typically S-P times 
between 19 an 22 sec. 

Seismic zoning maps of Canada, (Whitham and Milne 1972, 
Whitham et al 1970, Whitham 1975}, suggest that the South 


West Alberta and the Alberta Plains lies in the zero zone. 
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In this zone the annual probability that accelerations of 
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This probably 


reflects the low level of seismicity in South West Alberta, 


but it could also reflect the detection threshold of the 


permanent stations (M, near 3.5). 
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5.1 Previous Studies 

A seismicity map of Canada (figure 9) shows a great 
scatter of events along the Rocky Mountains. Historical and 
recent earthquakes have been located close to McNaughton 
lake in the Rocky Mountains (Rogers and Ellis 1979, Rogers 
et al 1980). This is the most active zone in the Rocky 
Mountain Trench. Some of this activity (events with 
magnitudes greater than 3.4) has been detected at Suffield 
(SES), Penticton (PNT), and Fort St. James (FSJ), and 
located by the Seismological Service of Canada, Department 
of Energy Mines and Resources, Ottawa. 

Previous studies of seismicity in Western Canada have 
been concerned with the offshore seismicity of British 
Columbia, its tectonic relation with major plate boundaries 
and the detection threshold of the Canadian Network in 
Western Canada (Milne et al 1978). This study is mainly 
concerned with the description and interpretation of records 
of micro-earthquake activity in South West Alberta recorded 
at the Edmonton seismological observatory (EDM). 

Of all the events recorded at Edmonton from the Rocky 
Mountain House earthquake swarm (table 15), fifteen events 
with clear Sn refracted phases were located by Ottawa. 

A recent example of activity in this area, the Willmore 
earthquake (frgure 30), occurred on October 9, 1977. It was 
located at 53.63°N and 118.29°W had a focal depth of 18 km, 
a reported magnitude of 3.3 (Ms) and 4.4 (mb). The 


McNaughton Lake earthquake of May 14, 1978, (Rogers et al 
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Figure 30.... Location of Snipe, Willmore and McNaughton 
events and some events located by EPB from the Rocky 
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1980) showed .8 right-lateral strike slip faulting combined 
with .5 of thrust movement. The preferred fault orientation 
is N10°W. The maximum principal stress axis is nearly 
horizontal oriented in a northeast direction. 

Earthquake swarms have been observed before and after 
the impounding of the McNaughton lake reservoir (Ellis et al 
1976). Ellis and Chandra (1981) analysed the seismicity of 
the McNaugthon lake reservoir. They found an average of one 
earthquake with magnitude M,23 every 2 years in the period 
of 1963 to 1972, and four earthquake swarms sequences in the 
period of 1973 to 1978. They suggest that the activity is 
mostly of tectonic origin. Earthquake swarms have also been 
observed near Bella Coola British Columbia and near Mould 
Bay (Milne et al 1969). Rogers (1981), based on the 
decreasing age of the Anahim volcanic belt from the coast to 
the Rocky mountains, the high level of seismic activity and 
the relocation of recent and historical earthquakes, 
Suggests the possible existence of an active hot spot 
beneath the Rocky Mountains, near McNaughton lake. 

In the Plains an earthquake occurred in 1909 in 
Southern Saskatchewan with a probable magnitude of 5.5 
(Stevens 1977). More recently an event occurred in 1968, 50 
km north-east of Bengough close to Amulet with a magnitude 
of mb=2.9 (Horner et al 1973). In 1972 the Bengough event 
occurred at approximately 49.35°N and 104.93°W. This event, 
in spite its small magnitude (mb=3.7), was throughly 


analysed by Horner et al (1973). They studied intensity 
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isoseismals, magnitude, epicenter location, focal depth 
(approximately 10 km), and focal mechanism. In order to 
calculate the focal mechanism es: used the radiation 
pattern of Love waves for a shallow event. The Bengough 
event was probably strike slip, striking approximately N 30° 
Be 

No historical earthquakes have been reported in the 
South West Alberta or the Alberta Plains. However, the Snipe 
Lake earthquake, of magnitude mb=5.1, occurred on March 8, 
1970. This event was felt over approximately 100,000 km? in 
west central Alberta. It was located near 55° N and 116.5° W 
with a probably focal depth of 9 km (Milne 1970). Even 
though the magnitude of this event was relatively high for 
this part of Canada it was not recorded clearly in many 
Stations. Therefore, no fault plane solution is available, 
and the probably cause of this event is unknown. It occurred 
in an area where oil and gas is produced. 

Milne and White (1958), using three seismic stations in 
the Crowsnest Pass area of Alberta and British Columbia 
during a period of three, years found a correlation between 
Goaliminingvactivityrand seismic events. sAneaverage, of five 
bumps per day was recorded in active mines. Using a small 
array (figure 30) at Fernie, Hasmer, and Anderson showed 
that the locations were confined near the coal mines. Local 
events detected at stations located at Fernie, Coleman, 
Turner Valley and Banff showed that there are 


micro-earthquakes (with magnitudes less than 3) extending 
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from the International Border to Banff. The distribution 
seemed random in the South West Alberta and in the Rocky 
Mountain Range. However, this could be due to a poor time 
control. No relationship between mine tremors and local 


earthquakes was found. 


5.1.1 Seismicity in South West Alberta as Seen by Edmonton 

Seismograms recorded at the Edmonton seismological 
observatory have been read on a routine basis since 1963. 
Only first arrivals were read from 1963 to 1969 but S-P 
times have been available since 1970 (table 16). 

A histogram of local events (S-P times less than 60 
seconds) shows all possible events from the South West 
Alberta (figure 31). Care was taken to avoid mine blasts, 
which usually occur during working hours and at fixed 
intervals of time. The largest events were Perera with 
S-P times at SES, PNT, and FSJ; generally located by Ottawa 
table 15. A relatively high level of seismicity was observed 
from 1970 to 1972 (figure 31-a). Then the seismicity 
decreased to an average of 2.2 events per year. Figure 31-a 
only includes what I call the "Normal Seismicity of the 
South West Alberta". It does not include the Rocky Mountain 
House earthquake swarm (RMHES), that apparently started in 
1976. Figure 31-b includes the Rocky Mountain house 
earthquake swarm. Events from the RMHES are easily 
recognised since they generally have an average S-P time of 


21.7+0.7 seconds at Edmonton. 


sp Lot chlis TG "Re Meth ai vies 
geil TO elae dans banat” aoe : Pia ial 


tH Seeyese Wor va - Mt 


a wt! ry iefure a et 


seek Eth alia ane eau Tost? satpedei 


ca RY eT i, | 
ro eed SH- afensin See dtuee hal gets naei mi 
aa Tt: rohit 7 oe waters minal gn a 


(ee: greets? Bee kaoslchse tend 
tat a 
faet wey’ boon Boke rata 
opal! BOY geeke sn ate via! 


we R- A) ssaeiee ao i nso 


fi, 


an! 


Pm stn or age ute oe shot rhs) ac 
. ify ie ia HNO ‘esa 
aa ened vied ah sab ehtes 
Levey wteyeebnan Mamma de Semel wii 
wt Wri] oe me ‘ae ie ea 358 
2 stale the vad ¢ " 


joe re menue? ra 


, ad i } 


igter a 


a tie? 


not, 9 wed 


Ley | 


eet AY, thu te hi see 
. an » 


pottnwe we 


MONTH DAY 


YEAR 

JAN1676 
JAN2476 
JAN2576 
JAN2676 
JAN2776 
JAN2776 
JAN2776 
JAN2776 
JAN2776 
JAN2776 
JAN3076 
FEB0876 
FEB0876 
MARO0676 
MAR0676 
MAR1076 
MAY0176 
MAY 1376 
MAY 1376 
JUN 1876 
JUN 1876 
JUN2176 
JUN2276 
JUN2776 
JUN2976 
JUL0876 
JUL0876 
JUL3076 
AUG0876 
SEP0176 
SEP 1376 
SEP 1476 
OGT15/7'6 
OGTA 7.6 
OCT157 6 
OCTZ376 
OCT 2776 
NOV0376 
NOV2976 
DEC0176 
DEC0276 
DEC0976 
DEC 1876 
DEC2376 
DEC3076 
DEC3076 
JAN0577 
JAN0577 
JAN0577 
JAN0577 


ARRIVAL TIME 
HOUR MINUTE 


SECOND 

oO: .6an 60S 
¥6.20:9 222 
VO eo/i0 
Rasen 
01.48.07 
04.13.00 
04.21.39 
OS caretct 
0513552 
OS 2S: colsl 
Wi Oi7s720 
18.49.02 
23.04.14 
G9: 43'0:15.0 
20.44.49 
03.49.54 
06.44.09 
04.48.04 
06.26.43 
Ne lekZiOi 5 by 
14.28.04 
Zileri e2o 
03338 ..014 
Li Sennes sre 8) S) 
00.40.43 
05.24.08 
05,226:.:0.0 
Pheer eho 
WGiet6...02 
PRES us eel 
10.24.48 
02.44.43 
LO soe. 05 
09.48.21 
Use Se) 
OZ ierSile 40 
GG.20.05 
11.42.54 
02.28.47 
Wis S011 
OStes01-n 10 
Oi722.075.0'3 
i OeonleeO © 
Z2e iO cel) 
(Pee. Si Paley 
Tearorie0n 
03.10.44 
03 5424 
10:548°535 
Zoic Ojenle 


N 
N 


NNN — 
ONOW 
7 e@«© e@© e # e@ 


N 
—- OO 


NN NY 
- OO 
e e e 


e e e e e e e e e e e e 
910190 MWWO—-WHOWO-—- A INITINMNMODOINNVRONDIOUWINODOWODOO0O 


124 


. e 


e e @ e e 


a e e 2 
WOONAVNDAIHNAOANHNDWOON KLPIWN -—UTPHYNWH D-H IW P= NHWOW 


WW NH]NN HH NWH NN HAH NWNHNNNWNHYDNYNNMNNN NH WH ND = fh = 
Cy at ee ee eo e@ e e e @ @ e 


+ 
. . 
we 
ay oe 
On 
eA Ts 
« at 
~y 
You 
. ‘ 
‘ 
a 
6 , 
tages 
4 
‘1,4 
r Pa) 
. 


= 7 
it 


ral | fe ae <> 7G 


i Sitar Th ¢:: 


oA - ee 
es = 


= 
i 
= tp i 


Sis Heeet cs coe 


4 


ew — “y me A = 


—— én a} Gu Sb 2 


o> +> 
~~. 
al 


ce we 


sa? 


i= — i Bin 
23534 


y 


= 
> 
= 


= © 
ee 
or 


=< % 
> 


4 
— . 


La 
om 


- ig 
13. #2 <4 =~ 


a = 
ay eo? Or 


~ 
2» 


5. aoe 


.. 


pee Le Moe 


ao. 


—?n 


a ge ae 
yi 


vis 


* 


as 
~~ ‘uf 
te ceed 


aa = 


"id: oe 


Bets i : 


Geta” 


JAN0677 
JAN0877 
JAN0877 
JAN0877 
JAN 1477 
JAN1977 
JAN1977 
JAN 1977 
JAN3077 
JAN3177 
FEB0Q377 
FEB0477 
FEB0477 
FEB1077 
FEB1177 
FEB1577 
FEB2477 
FEB2577 
MARO 177 
MARO577 
MAR2377 
MAR3077 
APR1077 
APR1277 
MAY2277 
JUN0477 
JUN 1877 
JUL0277 
JUL0377 
JUL0577 
JUL2477 
JUL2477 
JUL2477 
JUL2477 
JUL2477 
JUL2477 
JUL2477 
JUL2477 
JUL2577 
JUL2577 
JUL2677 
JUL2677 
JUL2677 
JUL2777 
JUL2777 
JUL2777 
JUL2777 
JUL2977 
JUL3077 
JULSd77 
AUG0277 
AUG0477 
AUG1177 


02.44.30 
04.40.13 
07% 058.36 
N2esaie ae 
h6%2.9e oe 
(Elba s hac 
144'24438 
Zia io 
OSe39525 
18.40.00 
17.30.48 
14.34.25 
168 26607, 
205.252 50 
05220632 
03.52.06 
18.42.28 
04.00.19 
07.49.39 
We abies 12 
VS) A) Shewets, 
RET ene he Oy 
04.26.11 
Welasos 
TGS Sine 
O6559% 37, 
04.06.46 
20240219 
10.44.37 
ZUR 4 0n 02 
O05. 56 
105 0S..06 
11.24.48 
IM332:4 
Tite 23 0 
6p. 06. 13 
Wir Ue 2c 
288 Oo 4 
Os. 3.6% 53 
08% 59.34 
OSe ttn Z5 
Ofer, Sere 13 
Oo lao: ia 
O2e 1c 09 
04.27.06 
07.54.01 
h2e 468 57 
12 Oe to 
1520, Gre 
022.587 56 
06.54.54 
04.06.03 
OBiioile 25 


oo 28©« © e@© e@ @ oe (See eee “se 61h (eo) sex ele “eo Fe. et 686. oe. Oo 6. oe 
NDNAIWODAINVNOMNWOWAAYIKHWDOMNMDODOWAO~T=—ON=-WODN 


NNNNNNMNNNNMN ANN WH ND NH W 


=’ IND) GOs ok ae 


W—= = = N= WNHNMN WN WN NN NH 
° e ° ® e e e @ ° 


e 
NNN -OWAIKPUONNWODWDAWO~) 


e ° e ® e ® @ @ e 


e e e e e 
mM NM WD OE WW JP = = = = ON 


® e e e ® e 
“S—-NOBOON VI 


e e 
Iw —- 


rooms 


VWZo 


& 
a 

J 

Sie oe 
ee a 
a, ) 
—_ 


. : a oy ° J 
; : a i ; roa 
an e yg 
ef Wie. } ) & miss j 
Me oe 
i aa dof E 
AD) 
e's i oe 
: at 
, > * 
, 4 
ae a lee aw 
5 , ¥ ri 
» v 
als 
tie 7 G5, 
} hs nN : 
ove 
# aN u% 
i ats 
& ee 
| aby 
oe S 
t ee aor 
; oo? i 
tak tht 
. Fis tht 
& ane a 
os pes 
\ yi 
Tee hat, 


eo aS eee oO ; 
; a 7 . 


e 
~~ 
—s 


- 


2555 & “J 
a Ph ia¥ Peat HS 


> 

e 
——~ 
; —— 


A] 
. 
a 


—— Le 


te igh ee itt 


pe: 
ek 
vet 
sey 
og $0 
mo. 4 eon as 
i a Pee, SE;50 10 9 
2 fos he re ranea: 
sae , i mG ‘i as oe. 82 
SUS eed. SF), Le 
a i : ft aed 82490 . - : 
a> Nd ean ack a7 'l 
ie Bk es Boy BH 

; : ab, es! G 


i 
—. -— 


—F 


= ya 
cai cm: 


wn 


, te 
A ‘ pA) 
' P ' 


- 
ad 
eS es 7 

= 


AUGT 377 
AUG 1477 
AUG 1477 
AUG 1477 
AUG 1477 
AUG1777 
AUG2377 
AUG2777 
SEP0877 
SEEN 7 
SEP h27 7 
OGTO2 77 
OCT0277 
OGT2Z6'7'7 
OCS] 7 
NOV1177 
NOV 1977 
DEC0877 
DEC1077 
DEC22/7 
DE@Z37) 7 
DEC2477 
DEC2877 
DEC3077 
DEC3 177 
JAN1378 
JAN1478 
JAN1778 
JAN1778 
JAN2178 
BEB S78 
FEB2678 
MAR1078 
MAR1478 
APR1478 
MAY2078 
JUN2678 
JUN2878 
JUL0378 
JUL0678 
JUL1878 
JUL1978 
JUL1978 
JUL2778 
JUL3178 
AUG0378 
AUG0878 
AUG0978 
AUG0978 
AUG1278 
AUG 1478 
AUG1578 
AUG1578 


OSs 2.2 
O09 sO 
Orincalenres 
O6n: 598819 


250 


NO 

At 
e e e e eo e e @ e e@* © ‘ene ©€ 6 ee, 6 “ese 6 e e@ e e e@ ee e Cy 
WOW PN SPW WrIYW = DOYODAYAHAD$AMHNHOLSFWOO NNN DOAN NM LO VIVO U1O OO 


- hb 


Nh 


NONMNNNWN-NNNNWNHNN —N W — —- NN W DY 
e e e ° ® ° e e e ° e e e e e e ° ° e e 


eo 
\O OV 


WWNHH]-NNN WN —]= H—NWNHNNN W = = — — 
e e e ° e e ° ° ° e e ° e ° © e e e e ° © 
Hm OVOV WO YIM SH ONDAIAYNPWWW POO ~AIN O10 


ONON ON PH HWOAWHA—-woOWd— & OO OI ~) NM 


—" 


126 


my rj 


eee & Sia 


~ 
ai, 


or 


* bs 
> ee 
2 ety 


c= eas 


orn 


~ 2 — 


we! 2 3 


oo 
— 


ce 


mn 


o 


err 1a 


AUG 1678 
AUG1778 
AUG2478 
AUG2778 
AUG2878 
AUG2978 
AUG2978 
AUG2978 
SEP 1278 
SEP2078 
SEP2678 
SEP2778 
OCT0278 
SCTad7s 
OCT2Z578 
DEC 1178 
DEC3178 
JANQ179 
JAN0279 
JAN0279 
JAN0579 
JAN1079 
JAN2879 
FEBO179 
FEBO979 
FEB0N979 
Pes iy 9 
FEB1179 
MAR2879 
MAY0479 
MAY2179 
MAY2579 
MAY 2879 
JUN0679 
JUL2479 
AUG0O179 
AUG0679 
AUG1279 
SEP0879 
SEP2579 
OCT0679 
OCP NY 
OCT IGA 9 
OCP, OF.9 
OCT 21779 
OG 2Hi/,9 
O€T2279 
NOV1679 
DEC 1379 
DEC2479 
JAN0380 
JAN0880 
JAN0880 


~42 
46 
«20 
ys) 
24 
03 


2g 
30 
21 


00 
“oo 


24 


04 
-09 
so 


32 
2 


44 


01 


Oo 


14 


239 
x02 


"3 


ees 
OF 
> OF. 
aD. 
46. 
47. 
AS Oya 
rey 
46. 
10. 
ees 
Bas 
DIS 
43 


2A 


seo 
OE 
49 
49 
oN 
LS 
ONE 
46 
4. 
ABeohe 
45 
oul 
my 
ROX Ses 


206 
52) 
50g 
20S 
sai) 
«20 
A | 
. 33 
A RS) 
47 
De) | 
o2i0 
Oe 
600 
07 
- 00 
ou6 
Pais 
06 
ooU 
<u) 
34 
<o8 
pei 
41 
Py 
43 
02 
17 
26 
«316 
05 
0 
Lo 
Sis) 
07 
01 
nets 
247 
54 
oft 
oie 
201 
30 
2a 
04 
wZ0 
34 
58 
A EAP 
45 
44 
30 


NM NM NN NM 
— — = CE © 
° 


NM NM MND 
-_-O-—-0O0 


N 
sk 


Nh MN NO 
SY eMOTs, 


MM NM NH NW NHN 
(e@ypicy 5 (ay V5 
se e es 


N BN 
Qos 


N NM NN LN — Nh 

(e) o-O-O oo 
e o @ e eo 8© e e @ e (8 e ° o—"@ Ver _*e %e7 te Ne ee e e e© @© @ e@ eo e © @ @ e @ e Cr ee ee | 
OWOODDWOWOUWOUONWWMOWUOWON UO UTUMONOWUFPN FLOW WAWYOOWYWROOO—o 


N 
a1 
e ° 
oOo 


NNNNM —NMNNNNN NM NN ND = — = YN W = HN ND lO 
e ce ° °o © e e eee ee e 


MOS] NNNNNNMN WH — ND —= NH 
eo ° Ce) 


NN — 


ON — NY = NN ND PP 


e e 
S=WOPYTO == =- ON —DOWWOWOTIO-—-WHoONWD 


® ° e ° e e @ e ° e 
ONorP-—-aNoO~ IWWOWOfHLWHhD 


e e e ° ° ° 
OWWO PWN ND — © 


COow 


a] + 
a 
4 
' 
ec 
=’ 
Sie 
Gq A 
“| 
¢ 
ie @ Ws 
ae 
é 
— 


me | 
a: = 
ee 
i 


= 
oot 
woe 


Sew 


ky 
Po 


Lhe 
we ee 


ee 7 
_——— - 


ere fo: 
4 


~ 


an or ee 


: 
7 Pam) poe = 
Mie WHF Te 


© Ou 


+ 
ey => 


hy. 
4 3 


- 
Ba ot 
oS 
ae i 
Oy iw f > 
$5 
(i « 
fF, 
@ 
¥ 
er a 
4 we 
an 


40 


tae 
/ bt eat 
i a8. ES 


a5 


JAN0980 
JAN 1080 
JAN1280 
JAN 1880 
JAN2380 
Pep 730 
FEB 1980 
FEB1780 
FEB1980 
MARO 180 
MAR0580 
MARO0580 
MAR0680 
MAR1980 
MAR2680 
MAY0580 
JUN0680 
JUN 1280 
JUL2580 
AUG0380 
AUG0680 
AUG2380 
SEP0180 
SEP0280 
SEP0280 
SEP0480 
SEP0480 
SEP0480 
SEP0480 
SEPOS80 
SEP0580 
SEP0580 
SEP0580 
SEP0680 
SEP0780 
SEP0880 
SEP0880 
SEP0880 
SEP0880 
SEP0980 
SEP0980 
SEP1080 
SEP 1080 
SEP1180 
SEP1580 
SEP1680 
SEP W/ 8.0 
SEP 1780 
SEP1780 
SEP 1880 
OETO 250 
OCT0980 
OCT2980 


~43.47 
~40.45 
~14.44 
a hf A of 


08225 


eC Lele 
~39.04 


WA is 


7o0. 04 
4939 
720.06 
54.50 
Fie) ota 
Pic AS 
eOS e053 


Tees 
05.34 


034.12 
~44.,43 
SOO. S| 
bt Pa iad 
See as 
eee ed 
Aol Gn) 
Agree 
-49.54 
“2 ONE RO) "2 


e 


34.22 


EIS seeps) 
rye) ee: 
-47.34 
-20.48 
-41.49 
200.06 
fe She wea Aes 
Seoke ey A 
oe Ocemic 
Aap PAY) 
047.34 
.07.48 
1oo8US 
Oia a 
-46.31 
.49.08 
~34.54 
Boe) Shea) 
[SSNs 
Bes) eK 
ee oales Ce 
e000 
niet s9) 2, 
2053.0 
Sree ies 


e 


. e ® ° ° e °e e e e e e e e e 2 e e e e e e e 
ONDODONODOLFPODWOWODMOWOBDWWWHANDOWO=—NM=DONWOND 


WWNHWWNHWWNHWWNYANNN = HENAN NANNNWWNHWANNND @2 NN =D = ooo oy oo ND 
e ° o ¢@ ° ©; '@¢" ¢ J¢. “se Cn) ry ee e e ee @ ° e oe e e® @ e 


N G& © 
. e 


e e 


e @ 


° e ) ° ° ° e ° ry r) e e e ° e 
WIPO OKDWWOSPOATPOW =A WOW DINDIIINNDOUNWOOWONKKHWAMNWWOWUDMDWDOUWWWN 


e 
— WO) 


128 


ag Oe ee —_ ie 
- incr a ee eee 
OF ie Se ie eee 


er 


fo 


gabe, 


- 


> --+-- 


rib 


Sy Phi 


- = 


i F 


OCT3 180 
OCT 316.0 
OCT S160 
0CT3180 
DEC 1380 
JAN1981 
MARO98 1 
APR1581 
APR1681 
JUN 1681 
JUN248 1 
JUL 1481 
JUL278 1 
JUL288 1 
JUL308 1 
JUL308 1 
AUG0 181 
AUG038 1 
AUG208 1 
AUG278 1 
AUG298 1 
SEP0181 
SEP048 1 
SEP238 1 
OCT158 1 
OGT98 71 
0CT2081 
OGT2 1651 
OCT228 1 
NOV058 1 
NOV1281 
NOV 1381 
NOV 1481 
NOV 1481 
NOV 1481 
NOV 1781 
NOV258 1 


TADL Cl 1 Ors oars 


00.40.37 
00:550%:23 
O5.1623 1 
16.237, 45 
09.31.48 
Wer oem 
12..00239 
eS 2S) 
12.0675 24 


2.0 
Zale) 
2075 
2020 
PA Ae) 
2050 
20S 
22.0 
ro 
20.0 
220 
20S 
EAGT 
Zo 
202.9 
20.8 
Zilte0 
ZO. 6 
Z20).6 
20a 
19.4 
7218) 218) 
24.0 
Z0n0 
21.0 
Zits 
Ziiew) 
ZO0mo 
Ss, 
Zils 
Zaliets 
2S 
ZOE 
20.8 
PAS) ARS) 
PSRs: 
20.8 


NINN W=| WN NNNNN HNN HH NWNNNWNHNDN —H— NNN — NNN NNN NH 
e e e e e ° e ° ° ° ° e e e e ° ° ° e ° e ° e e ° e e ° © e e e ° e ° e 


O- —-$ OWO FF WMONAVNAWNWO OO WO HA OMOWOKHTD AWAIT LHS OS 


Table of all the events recorded at Edmonton 
from the Rocky Mt. House earthquake swarm 
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Figure 31.... a)Histogram of local seismicity recorded at 
EDM including RMH earthquake swarm. b) Without RMH 
earthquake swarm. c) Histogram of some well defined S-P 
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5.1.1.1 Relation to the Mica Array 

Aside from the station at Edmonton the nearest seismic 
facility to the Rocky Mountain House Earthquake swarm is the 
Mica Array. I have examined an extensive sample of the 
records from this analog array of vertical seismometers and 
conclude that the S phase is too emergent as a rule to aid 
materially in refining locations by means of S-P time. In 
any case, exact locations are of secondary importance to 
this study. Its thrust is towards source mechanisms. The 
Gata at Mica are all recorded on analog magnetic tape, 
Spectral studies such as I report here would require 
uncertain corrections for non-linearity in this recording 
medium and after the fact digitization. I have examined 
analog records very carefully; convincing evidence of the 
presence of the Moho refracted phases cannot be found from 
analog records alone. Finally, and most important, the ray 
path from Rocky Mountain House to Mica passes through 
considerably more complex geology than that to Edmonton. 
Even if the data could be reliably reduced its 
interpretation would be substantially more uncertain. 
Nevertheless I report here those results I have from Mica 
data. 

Some events recorded at CUM station from the RMHES have 


an average S-P time of 24.5+0.8 seconds (table 17). 
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Events from the Rocky Mountain earthquake swarm 
detected at CUM 


5.1.2 Summary 


From the histogram of number of events against S-P 


times from the figure 31-c, and associated locations (Table 


15) we can suggest. 


1) Events with relatively small S-P times must lie in 


the Alberta plains (figure 32). Whether these ‘events are of 
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Figure. 32.1... WocCaleeVvents recorded at EDM. Right, " an 
event with S-P=12 sec. Left, . - an event with Sg-Pg=27.3 
and Pg-Pn=2.4 and was located by EPB at 51.95N and 115./76W. 
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tectonic origin or due to the oil or mining operations in 
the Alberta plains cannot be answered by this general 
discussion of the seismicity. 

2) Events with S-P times between 20 and 30 seconds are 
probably close to Rocky Mountain House, approximately 180 km 
SW of Edmonton (figure 33). This is the most active zone in 
South West Alberta. I showed in chapter 3 that some of those 
events were mainly deep (nearly 20 km). However, most events 
detected with a temporary array of analog and digital 
seismic stations were not deeper than 4 km (Wetmiller 1981, 
Rebollar et al 1981 b). 

There may be two kinds of activity associated with the 
Rocky Mountain House earthquake swarm. Deep seismicity (near 
20 km) more likely of tectonic origin, and shallow 
seismicity (not deeper than 4 km) possibly associated with 
secondary recovery activity in the Strachan gas field. 

3) Events with S-P times around 35 seconds 
(approximately 300 km from EDM), could come from any place 
between Calgary and the Willmore park area. Therefore, more 
evidence from a close network or other stations is needed 
(table 18). 

4) S-P times between 40 and 60 seconds probably 
correspond to earthquakes that come from the area of 


McNaughton Lake (figure 34). 
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Figure 293... .cnxampler oLmthes Rocky Mountain House events 
detected in the analog stations at EDM, SES, and PNT. Some 
of those seismograms do not show clear phases. 
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Figure 34..,. Bvents with) S-Potimes greater than’ 40) sec 
probably coming from the McNaughton lake area. Right event 
fas So-Pq=44,.S0-Sh-=i 1] 5. vancduedgren=7. 7 ssec helt event has 
Sg-Pg=46 and Pg-Pn=11 seconds. 
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5.1.3 Energy release and b value of the Rocky Mountain 


Earthquake swarm 


Lacking an empirical relationship for the evaluation of 


radiated seismic energy of small earthquakes in Canada. I 
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used Gutenberg and Richter's (1942, 1956) energy-magnitude 
empirical relationship in order to calculate the energy 
release of the Rocky Mountain House earthquake swarm. 

Gutenberg and Richter considered the radiated energy 
contained in a spherical shell centered on the earthquake 
source. They assumed that this energy was radtrated ina 
Sinusoidal wave train and recorded on a standard strong 
motion instrument. Therefore, knowing epicentral distance, 
hypocenter, the physical constants of the standard 
intrument, acceleration, period and travel times of the 
Signals, they found an empirical relationship given by 
logE=11.3+1.8M,. Later, Gutenberg and Richter (1956), ina 
review of their original paper found the relationship given 
by logE=9.4+2.14M,-0.024M’,. 

This empirical relationship depends on the theoretical 
study of seismic radiation at short epicentral distances, 
therefore it can give a reasonable estimate of the energy 
for small earthquakes from the Rocky Mountain House 
earthquake Swarm. 

An average of 10'* ergs per month was released during 
the 58 months since 1976. The contribution of small 
earthquakes to the total energy release is negligible 
(figure 35). The total energy released during this period 
was 5.6x10'’ ergs equivalent to a single earthquake of 
magnitude (Ms) 3.9. The total energy release of possibly 
deep events, i.e. earthquakes that show clear Sn refracted 


phases, was 5.49x10'’ ergs. This mean that 98 % of the total 
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Figure 35.... Energy release, cumulative energy release, 
and strain release from the Rocky Mt. earthquake swarm 
calculated at EDM since 1976 to 1980 
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energy release was released by possibly deep events. This 
mean that largest earthquakes (with magnitudes between 2.8 
and 3.5) are mainly deep. The cumulative seismic strain 
release (figure 36) yields a maximum strain release of 
Peoxt0* ergs®*, 

The strain release as a function of time (figure 35-a) 
shows a maximum release of strain energy during the first 
year of the swarm (1976). After that the strain is 
accumulated and represented by the flat part of the plot. 
However, whether the strain energy is released in small 
earthquakes or it is accumulated, cannot be answered from 
this short period of observation. 

Using 242 events I calculated the cumulative number of 
events versus local magnitude (figure 36) or the Gutenberg 
and Richter frequency-magnitude relation given by 
logN=a-bM,, where N is the number of earthquakes for unit 
time, a and b constants and M, local magnitude. I used some 
of the magnitudes Ss enlenes by EPB and reported in their 
bulletin. For events not reported by EPB, I calculated the 
local magnitude according to Richter 1958 (see for example 
chapter 2). 

There are three main factors that limit the accuracy of 
the evaluation of the b value. The magnitudes are uncertain, 
small events are undetected, and there are few events due to 
the short time of observation (Milne et al 1978). 
Uncertainties in the evaluation of magnitudes are difficult 


to estimate if the magnitudes are calculated at a single 
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station. 

Edmonton station detects events from the Rocky Mountain 
House area, with magnitudes greater than 1.2, aS can be seen 
from the frequency magnitude relationship, however, the 
system loses its sensitivity for magnitude 2 from this area. 
The b value for the Rocky Mountain House earthquake swarm is 
0.8 and lies in the range of values (between 0.6 and 1.5) 
found by Everden (1970) using world wide seismic data, and 
is Similar to the values found by Milne et al (1978) in 
Western Canada (those values range from 0.65 to 0.82). 

The frequency-magnitude relationship has been found to 
vary in different tectonic regimes. Higher values of b have 
been found in spreading centers (Sykes 1970, Reichle 1975) 
and in volcanic areas (Suzuki 1959). Small values, near 1.0, 
are observed for intra-plate earthquakes, like this 
sequence, and in transform faults. This parameter usually is 
associated with the state of the stress (Mogi 1962, Scholz 
1968, Berg 1968, Wyss 1973). However, Knopoff (personal 
communication) suggests that it is more likely to be 
associated with the roughness of the fault. 

An analysis of the bulletin edited by the Earth Physics 
Branch of Ottawa, and the archives at Edmonton reveals that 
the Rocky Mountain House earthquake swarm is not a common 
geophysical event. The common pattern in the South West 
Alberta is scattered events without an aftershock pattern, 
much like the Willmore event (Ms=3.0), and the event (local 


magnitude 3.4) that occurred on December 22, 1981, with an 
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onset at Edmonton at 11 hours 36 minutes and 51.2 seconds. 
This micro-earthquake was located by E.P.B. at 51.95°N and 
115.76°W, approximately 220 km from EDM (see figure 32). 
However, this pattern of events without aftershocks could 
indicate that the aftershocks are too small to be recorded 
in the permanent stations. Therefore, the Rocky Mountain 
House earthquake swarm is an interesting phenomenon that 


needs to be studied more closely. 


5.1.4 The Willmore Earthquake 

Earth Physics Branch (EPB) located the Willmore event 
at 53.63°N and 118.29°W with a shallow depth and assigned an 
Ms=3 an mb=4.4 (see Nuttli, 1973), The’ U;. SS. National 
Earthquake Information Service (NEIS) located the event at 
53.68°N and 118.03°W and assigned a mb=4.4 (figure 30). 

This event (October 9, 1977) was not recorded clearly 
in many stations, nevertheless I plotted ten first motions 
(table 19), assuming a focal depth of 18 km for this event. 
Obviously there is no unique solution for the focal 
mechanism. Among those solutions it is possible to fit a 
Strike slip fault, with either of two possible fault plane 
Orientations, one of dip 60> andea dipi direction of M207 and 
the other with a dip of 60° and dip direction of 24°. None 
of those agree with the strike of the Rocky Mountains. A 
normal fault with either a dip of 20° and dip direction of 
Sue tor’a dip 70. andr dipedirection Of (2745 16Walso possible. 


(figure 37) 
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Seismic Station 
Mia Dainard (DAT) 
Mt. Thompson (THO) 
Mt. Cummins (CUM) 
Mica Creek Village 
(MCV) 

Edmonton (EDM) 
Penticton (PNT) 
Victoria (VIC) 
Yellowknife (YKC) 
Inuvik (INK) 
Mould Bay (MBC) 
Albuquerque (ALQ) 


Latitude 
52.199°N 
52.689°N 
52.087°N 
52.008°N 


5a. 222 N 
49.317°N 
48.518°N 
62.478°N 
68.307°N 
76.241°N 
34.942°N 


Longitude 
118.384°W 
119.120°W 
118.212°W 
118.561°W 


113.140°W 
113.140°W 
123%9417°W 
114.473°W 
133.520°W 
119.360°W 
106.457°W 
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Table 19.... Table of the seismic stations used in the plot 


of the equal-area projection of the Willmore earthquake 
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DD = 274° 


® COMPRESSION 


(S) RAREFACTION 


Figure 37.... A plot of ten first motions in an equal-area 
projection of the focal sphere showing probable fault 
Orientations (D=Dip and DD=Dip Direction) 
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Edmonton station recorded this event in digital and 
analog format. I was able to calculate the seismic moment. 
Only the vertical component was not saturated at EDM (figure 
38). Therefore, I used the spectrum of the vertical SV 
component of the shear waves in the analysis. The Sg-Pg time 
recorded at Edmonton was approximately 37 seconds. This 
gives an epicentral distance of 325 km assuming a half space 
with a P-wave velocity of 6.4 km/sec. The EPB location is 
330410 km from Edmonton. The spectrum was calculated 
assuming a Q of 1000, independent of frequency, a density of 
2.9 gr/cm*, an average radiation pattern of 0.63, a shear 
wave velocity of 3.7 km/sec, an epicentral distance of 
325+10 km/sec, and an angle of incidence of 49°. Figure 38 
shows the spectrum calculated at EDM. The logarithmic low 
level of the spectrum for the Willmore events is -3.1, and 
that gives a seismic moment of 6.7+2x10?? dyne-cm. 

In order to investigate the difference between the 
seismic moment calculated with body waves at a single 
station and the seismic moment calculated from surface waves 
at several stations, I used the McNaughton Lake earthquake 
‘of May 14, 1978. The moment calculated using the amplitude 
equalization method (see Aki and Richards 1980, chapter 7) 
was 4+2x102° dyne-cm (Rogers et al 1980). At Edmonton 
(figure 38) this event saturated both horizontal components. 
Hence, assuming the same physical constants as those used 
for the Willmore event, and the seismic moment calculated by 


Rogers et al (1980), I get a logarithmic low level amplitude 
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Figure 38.... Seismograms and spectra of Willmore (Ms=3) 
and McNaughton (M,=4.8) events. Upper (Willmore) event show 
EW and NS saturated components. Bottom (McNaughton) shows 
the vertical component at EDM 
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of the spectrum of -2.4. The value observed at EDM is -2.2, 
which yields a seismic moment of 7.6x10?° dyne-cm (figure 
38). However, the low level amplitude of the spectrum could 
have an error due to the bandwidth of the digital system of 
Edmonton. 
If we consider the orientation of the fault (strike N 

10 W, dip 38 SW, .8 Strike slip, .5 thrtist), the correction 
for the radiation pattern is 0.45 (see Aki and Richards 1980 
page 115). Considering this correction I get a seismic 
moment of 107* dyne-cm, if no correction 1s made for the 
radiation pattern we have 4.8x10?° dyne-cm, therefore, the 
seismic moment calculated at EDM gives 7.6+2x107° dyne-cm, 
two times higher than that calculated by Rogers et al 
(1980). Kanamori and Anderson (1975) comment, 

"The seismic moment is one of the most reliably 

determined instrumental earthquake source 

Parameters. For large earthquakes, in particular, 

the value of M, determined by different 

investigators seldom differ by a factor of more than 

two. For small earthquakes, the uncertainity is 

usually somewhat larger". 
In view of this comment I consider the agreement 
Satisfactory. However, probable errors in the evaluatiuon of 
these seismic moments are: errors in the epicentral 
distance, radiation pattern, attenuation, path 


inhomogeneities, and in our case a probably narrow 


bandwidth. 
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5.2 On the relation of Source Parameters to Tectonics 

In this work I have derived a number of measures of 
earthquake behaviour in various areas. Their connection to 
geology, which was discussed above is not easy to establish 
unambiguously. The physical properties of the system are 
Simply to complicated and the quality of the data is not 
hagh, 

First an obvious question arises: How realistic are 
those source parameters? I cannot answer this question 
precisely, since there are many unknowns in the source 
region, and in the earth's crust in which the signals 
travel. Corrections for attenuation seem to be not critical, 
Since Q has been observed experimentaly to be between 100 
and 1000 for most earth materials, in a range of frequencies 
from 10-5 to 10° seconds (Pilant 1979). A major problem in 
our observations is the loss of energy due to scattering and 
conversion of phases that could give rise to changes in the 
estimates of source parameters with a given model. Chapter 
two shows that source parameters calculated at distances 
greater than the fault length are usually smaller than those 
calculated at distances of the order of the fault length. 

I studied two different earthquake sequences, at 
relatively large epicentral distances (180 and 500 km). Even 
though the amplitude response curves are different, for a 
typical RESMAC station and the Edmonton station, the useful 
information in both systems was between 0.2 and 7 Hz. 


Outside this range of frequencies the ratio of signal to 
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noise was small. Therefore, in order to get more 
representative information about the earthquake source, it 
will be necessary to complement spectral analysis of body 
waves with spectral analysis of surface waves when the 
epicentral distance exceeds about 150 km. A typical example 
of this kind of study is that done by Hanks and Wyss (1972). 
The seismic moment of the Willmore and McNaughton lake 
earthquakes were calculated from the spectral analysis of SV 
waves recorded at the Edmonton digital station. Those 
moments are 6.7+2x10?? dyne-cm for the Willmore event and 
io r2xwO07e dyne-cm, for, the McNaughton lake event. This last 
seismic moment compares moderately well with that calculated 


with surface waves (Rogers et al 1980). 


5.2.1 Speculations on Tectonic Implications 

Results from individual earthquakes or earthquake 
Sequences, are important to an eventual understanding of the 
earthquake source and plate tectonics models. Orientation of 
principal stresses (calculated from fault plane solutions), 
magnitude of stresses (calculated from in-situ measurments), 
seismic moment and source dimension (calculated from 
spectral analysis), are important parameters of an 
earthquake source and the values of these parameters are in 
the final analysis controlled by the driving mechanism of 
plates. Seismic moments and source dimensions have been used 
to calculate the rate of slip of major plate boundaries 


(Brune 1968), and time recurrence of large earthquakes 
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(Singh et al 1981). 

Intraplate earthquakes have been observed to lie along 
preexisting zones of weakess within areas of youngest 
orogenesis (Sykes 1978), or to be apparently due to effects 
of local topographic features in the plates (Stein 1978). By 
close analysis, Sykes (1978), found that seismic activity 
tends to be concentrated near ends of major oceanic 
transform faults, along preexisting zones of deformation, or 
along faults in old fold belts within the thicker 
lithosphere of the continents, hence, the study of 
earthquake sequences in intraplate enviroments may reflect 
varioations in strength within plates. 

From spectral analysis of seismograms we can derive 
stress drops, seismic moment and source dimension. 
Earthquake stress drops of large intra- and inter-plate 
earthquakes are probably not greater than 100 bars. This is 
based on a large number of observations independent of 
source strength. McKenzie and Jarvis (1980) estimated a mean 
stress on plate boundaries of 150450 bars with a probable 
upper bound of 500 bars, in agreement with seismically 
determined stress drops. Higher stress drops (630 bars) of 
micro-earthquakes in the sedimentary basin of the Imperial 
Valley of California have been observed (Hartzell and Brune 
WOT e7L ee 

In this work I compared an interplate sequence of 
events (Oaxaca aftershocks) and an intraplate sequence (the 


Rocky Mountain earthquake swarm). They differ in some 


once wit oo RRQ Rneeed ‘vad, whiiig ‘x: | 
Jeapatioy Oe nae eats be gy 

aye ele mes ee on) on 3a a & 

A a YG? See bh: Wa 4s. ee ahh 


opiate cade Sao? aay) 


— 


— 


My pam aceon ee ee aed on wee cananiaes: por 
icgemtiokes To asnder booed gaa! papal (uatve® ce 
rei nhae sto eee ee aaled thot bap nd ett | 
1? eon |} sae a ae x0 
Seal Pere ee wi nh PT EN a bd dun hy fey, Hestiavnas | 
Seeteiy my agent at } | 
(R250 ey ee ele < clbgrd Re oe  % og ag serevee | - iy 
extents |S uate ha ety via tae Sab 
gherst rs ok ea ee mga wogiae ste 
ATES NOT MES “sea kercthe pidedore wie set " 
i¢ 348 coe, BHR Eh poe aye lave apaed re na es 
md here le 1a ark ea AMR deere 
Bekah A OR AROS TES, t Siteth eT ‘ogg sd 20 #15 iQ) tio 
Yi eas OBESE sib OC oa ba a yee Gn We esasead 
to Cezan Bagh Se gi ade sate weet Ree TIS neared 
Ys! es ténn | wagon ayo) ina aE 
equié Bin 2-2 4-0) Gate see ee ei det La Bee (etiev 


ed 


ue 


| Arsen. 
i-i$wwnes wir eee ite Egeaied 2 21c8 gift wt 


atts} "alien yanz Tey (hee es (Pagede ssi is gsaxso) esnana. ? 
paren 30h (om ee Bhmipcs ie. nietavens £20, | 
r + _ 


TS2 


respects. The Oaxaca aftershocks have simple spectra with a 
well defined plateau at low frequencies, a recognizable 
corner frequency and a well defined asymptote at high 
frequencies. Therefore, stress drops calulated in the usual 
way can have meaning. A number of the events of the Rocky 
Mountain House earthquake swarm have more complicated 
spectra, even though the signals were recorded closer to the 
earthquake source than the records of the Oaxaca events and 
the structure between the Rocky Moountain House area and the 
Edmonton station is not complicated. This difference can be 
appreciated from the plot of local magnitude against seismic 
moment. The data from the Rocky Mountain House earthquake 
Swarm have more scatter than the data from the aftershocks 
from the Oaxaca event. 

This difference could result from differences in the 
state of stress in the crust and differences in the 
heterogeneity of the fault zones. The Oaxaca earthquake was 
a large simple event (Stewart et al 1981), and it has been 
found to be similar to past events in the Middle America 
Trench (Chael et al 1982). Apparently simple events can be 
associated with relatively homogeneous stress states and 
physical properties on the fault zone. Maybe the complicated 
spectra from the Rocky Mountain House earthquake swarm 
reflect an incrpient) fault or a possible reactivation of an 
ancienteraulcusvocem awh che cOoULldmbe sexpeccrecmcomnave ja more 


complex behaviour than a well established fault zone. 
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Although this study has not answered many of the 
problems of cause of the Rocky Mountain Earthquake swarm it 
appears a Supportable hypothesis that 2 kinds of seismicity 
exist. There is small magnitude shallow seismicity and 
deeper seismicity probably associated with pre-Cambrian 
basement structures, which accounts for most of the energy 
release in the area. Many unknowns remain. For example, we 
do not know fault orientation, stress orientations, or a 
precise focal depth of these events. Nevertheless, whatever 
their values clues to the nature of the process exist in the 
forces acting on the lithosphere in this area and the stress 


strain relations obeyed by the lithosphere in this area. 


5.2.1.1 Comments on Quantitative Aspects of Plate Dynamics 

There are many poorly understood problems associated 
with the dynamics of the crust and upper mantle. The driving 
mechanism of plate tectonics, energies involved in the 
crustal faulting, magnitude of deviatoric stresses (total 
stress minus hydrostatic stress), and absolute stresses, are 
only a few. 

Solomon et al (1980) summarized possible driving forces 
of plate tectonics. Among those they consider lithospheric 
cooling, latitudinal plate motion, crustal thickness 
inhomogeneities, lithospheric loading and unloading, plate 
boundary forces at ridges, trenches, transform faults, zones 
of continent-continent collision, and basal forces 


associated with viscous interaction between the lithosphere 
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and asthenosphere. 

Different magnitudes of deviatoric stresses in the 
crust can give rise to different phenomena, such as 
micro-earthquakes, large earthquakes, and topographic 
features in the crust. Hanks and Raleigh (1980) argue that 
away from plate boundaries the crust can support large 
deviatoric stresses, of the order of kilobars, in response 
to local or regional loads. For example, McGarr et al (1979) 
found evidence of high stress drops, greater than 700 bars, 
in deep mines in South Africa. Scholz (1980) argues that in 
some faults stresses of the order of a kilobar may exist. 

Zoback et al (1980) found that shear stresses 
determined in wells near the San Andreas fault increase with 
depth. The depth of those wells ranges from 200 to 1000 
meters. Apparently stresses increase linearly with depth at 
least up to 5 km. McGarr (1980) found that the maximum shear 
stress increases with depth from nearly 10-60 bars (1-6 MPa) 
in the first hundred meters, to 200 bars (20 MPa) to nearly 
5 km. He calculated a theoretical relationship of the form 
A+Bz, where A and B are constants and z is the depth, for 
the increase of stresses with depth. There is great scatter 
of data in the first kilometer or so (see also for example 
Brace and Kohlsted 1980). The constants vary with tectonic 
environment. 


McGarr (1980) comments 
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"A vertical profile of measurements from the surface 
downward might show a monotonic change of the 
horizontal principal stress directions, with the 
implication that the stress trajectories measured at 
the surface may not be indicative of directions 
throughout much of the crustal section". 

Bell and Gough (1979), and Gough and Bell (1980) found, 
from in situ meaSurements of oil-well fractures (breakouts) 
in Alberta, a remarkable aligment of the maximum principal 
stress in the northeast-southwest direction, i. e. 
perpendicular to the Rocky Mountains. However, those 
measurements were made in soft rock (siltstones, sandstones, 
and carbonate sediments) and at depths not greater than 2 
km. Therefore, it could be possible that shallow earthquakes 
are triggered by those stresses in Alberta, these shallow 
events are those events with S-P times of less than 15 


seconds and recorded at EDM, or the shallow events from the 


Rocky Mountain House swarm analysed in chapter 4. 


5.2.2 A Speculation on the Cause of the Rocky Mountain House 
Earthquake Swarm 

The suggestion that seismic activity originates in the 
pre-Cambrian formations is controversial. In order to 
explore its implications I consider, in a very superficial 
way, the tectonics of South West Alberta. Two structures 
could be relevant. They are the North American Cordillera 
and the Fond du Lac gravity trend (Walcott 1968) 

The North American Cordillera runs from the Aleutian 
arc through the western third of North America into Mexico 


and Central America. The part of the North American 
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Cordillera that extends from the 49th parallel to the 
Tintina Trench, near the border with Alaska (The Canadian 
Rocky Mountains), is a relatively simple and well studied 
orogen. Its age has been determined as Paleozoic-Mesozoic 
(Thompson 1979) and its geology and tectonics have been 
extensively studied (Bally et al 1966, Thompson 1979, Martin 
1963, Bokman 1963);=Shaw 1963, Elliott 1976, Dahlstrom 
1969,1970, Monger and Price 1979, Beaumont 1981). The 
Canadian Rocky Mountains are delimited on the east by the 
Interior Plains or Great Plains and to the west by the 
Western Omineca CrysStaline Belt ,and they are usually 
divided into the Foothills, Front Ranges, Main Ranges (East 
and West) and the Western Ranges. 

The Foothills consist of deformed Paleozoic Mesozoic 
clastic rocks. The deformation (shortening) is mainly thrust 
faulting in the southern part and folding in the north. Fold 
axes and fault strikes are parallel in most cases. Tear 
faults and normal faults are also common in the Foothills. 
The basement lies at approximately 3.0 km in the Interior 
Plains and dips gently to about 11 km below the Rocky 
Mountains (Bally et al 1966). The approximate depth below 
the Rocky mountain House area is 4.8 km based on reflection 
profiles (Bally et al 1966). 

The Front Ranges are formed from a relatively few major 
thrust layers that are bounded by faults with several tens 
of kilometers of displacement. The Main Ranges consist of 


lower Paleozoic and Proterozoic rocks thrust eastward over 
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the adjoining Front Ranges. The Rocky Mountain Trench is a 
narrow linear topographic depression extending over 1600 Km. 
In the south the trench is underlain by a westward dipping 
basement, and the location and strike of the trench is 
controlled apparently by a complex system of curved 
low-angle normal faults (Bally et al 1966). 

Although conventional wisdom suggest: that the 
Cordillera in this area is decoupled from the pre-Cambrian 
basement at the decollement zone (Bally et al 1966), 
Significant stress concentrations in the basement can result 
from the load of the Rocky Mountains and the adjacent 
foreland basin. (They might also be present as possible 
residuals from the Laramide orogeny that gave origin to the 
Rocky Mountains). Beaumont (1981) considered a model of 
lithosphere flexure for the Alberta foreland basin under 
laterally migrating loads and predicted the bend of the 
basement, and consequently the deep root of the Rocky 
Mountains. | 

Monger and Price (1979) compiling the work of several 
authors suggest a Moho depth of 50 km below the Rocky 
Mountains adjacent to Rocky Mountain House. Therefore, this 
possible flexure of the crust could create enough stresses 
to generate micro-earthquake activity in the basement of the 


Rocky Mountains. 
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5.2.2.1 Is There a Relationship to the Athabasca Axis? 

There is a gravity low extending northeast 
perdendicular to the Rocky Mountains (Douglas 1970). This 
anomaly extends from the area of Rocky Mountain House 
southwest of Edmonton Alberta to North Saskatchewan, and 
North West Territories. This feature is referred as the 
Athabasca axis in north Saskatchewan (Darnley 1981) or the 
Fond du Lac gravity low (Walcott 1968). This gravity low 
correlates with granitoids containing elevated levels of 
radioelements (Darnley 1981). 

There exists geologic speculation (Darnley 1981) that 
this rock belt is the result of intrusion into a zone of 
tensional faulting that developed before the end of the 
NewermProterozoic4 Lf this. “isin *factithetcase,.trtis 
reasonable to expect remnant zones of weakness in this 
portion of the lithosphere. Walcott (1968) studied a portion 
of this anomaly in Northeastern Alberta and Northern 
Saskatchewan. In order to explain the Fond du Lac gravity 
low, he proposed a crustal thinning, perhaps related to 
tensional stresses in the past, to 34km from 40km with a 
relatively low density below this anomaly. Such lateral 
density anomalies could lead to stresses related to 
isostatic compensation (McNutt 1980). In other regions such 
stresses may lead to seismicity, but analysis of the 
archives of the Edmonton seismic station (McGavin personal 


communication) shows no seismicity along this gravity 


anomaly. 
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The exception to this statement occurs where the Fond 
du Lac trend intersects the Canadian Cordillera. Here, in 
the region of Rocky Mountain House, the combination of 
loading stresses from the foreland and isostatic 
compensation stresses due to the crustal density gradients 
could conceivably activate zones of weakness in the 
pre-Cambrian crust. Such zones of weakness need not be large 
faults but they imply a degree of brittle yielding in the 
lithosphere at this point. 

I do not suggest that major fault zones exist in the 
Pre-Cambrian basement, for SW NE reflection profiles in 
Southwest Alberta imply that the crystalline basement has 
not been involved in the shortening of the Canadian Rockies 
(Balley et al 1966). Structure contour maps of the top of 
the Devonian and of the top of the Sub-Cretaceous show NW-SE 
and NE-SW trends in central Alberta. Those trends could be 
faults (Robinson et al 1969). 

I conclude therefore that deeper parts of the basement 
below the South West Alberta can undergo some deformation 
and consequently generate the deep seismicity I suggest is 


associated with the Rocky Mountain House Earthquake swarm. 


5.3 Summary of Conclusions 
Reasonably reliable source mechanism studies can be 
done with single digital stations at some distance from an 


earthquake. Source parameter estimates can be related 
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systematically to those derived from close-in records. In 
general seismic moments derived from body waves at distant 
Stations are smaller than those derived from near stations 

The useful band of information on 2 vastly different 
Systems, the RESMAC system, and the EDM station, turned out 
to be very much the same. This reflects a property of the 
seismic noise background and should be useful for other 
analyses of local networks. 

Localized zones of seismic activity rather than a 
random distribution (Milne and White 1958) exist in South 
West Alberta. One such zone, the Rocky Mountain House area, 
is characterized by shallow events which might be related to 
hydrocarbon recovery operations and deeper events which may 
reflect yielding of the lithosphere under a combination of 
isostatic compensation stresses and the load of the Canadian 
Cordillera. More detailed research should be done with local 
networks in order to map this activity. 

The shallow events at Rocky Mountain House are 
characterized by high stress drops, the deep events have low 
stress drops. This suggests distinctly different properties 
of the fault zones and supports the separation of this 
SCCl witvanntome durterent. types -ithe shal lowsevents are 
probably related to secondary recovery processes in nearby 
gas fields, the deep activity may reflect reactivation of 


ancient zones of weakness in the pre-Cambrian basement. 


et brane ren pete aD Serkan 
eqezeib 2a ceeew gio es) pevineh' 
sooidease sae wus hevires seadr pee 
jantg37he i these hei sonroInk Bi 2 b -xotaal 
fod Baal ¢ ee tt Ste ss eym a 
iY oven aed: abt + Sram wa Ye 

Wet 93. tee ee ods bas: bree 4 9 
caper anol 

ii chhties ie ime faa te. ‘aeaoe suit 

dh. talwe pei pit bos ey patil | 

nia URE. rohit na wd. afc, aOR Hong way atts i ‘ 
statde se! ore a gee etches veateda lai dpa sae 
vine ko Lie lesa toga shu TRe, ong tadegs, aevgeae | rite ' - 5 
jriadtihee > Set a or skigeatia Bf ea . “ 
ober: aah band gag tive ometetae nelsmen Ros, 5 
sot t dd ame te Riwege (adder, Neg hibeaae oxait i 
Pf ct 2? tag “er . oo + Sh30) wd i 
we ems ‘nded eos ae amt in nates voLiedn =o 
wol 3V68 3700 gad’ dD geet reer) fei ya aesiietg i 
asi sinyedg needa ih. giao ‘take aieepeue a ight pare: 
| Site ad inet th paged ‘ans a7 yorRgue bas aay uae 
ate aihews aXe sit ase Megat t2b S joss v4 £235, 
yoveata0! adatom’ gteveger Tapes ot Setsles vided 
+o’ Halts. + yr ome" pce eae oil doe “ooh ac3 ‘bie beat 
eng! tc relddeceg ig sade ib aendaay., Ze patos ae 


6. References 

Aki, K. Scaling law of seismic spectrum. J. Geophys. Res. 
2A 217 = 123 ty eg6 7 t 

Aki, K. and Richards, P. G. Quantitative Seismology Theory 
and Methods. Volume I and II. W. H. Freeman and Company 
San Francisco, 1980. 

Anderson, D.L., and Hart R.S., Q of the earth. J. Geophys. 
Res, 037586945662 7 ers7e 

Archambeau, C. B., General theory of elasto-dynamic source 
fields, Rev. Geophys. Space Phys. 6, 241-288, 1968. 

Archambeau, C. B., Developments in seismic source theory. 
Rev. Geophys. Space Phys. 13, 304-306, 1975. 

Bakun, W., and Bufe, C. G., Shear-wave attenuation along the 
San Andreas fault zone in central California. Bull. 
Seism. Soc. Am. 65, 439-459, 1975. 

Bath, M. The energies of seismic body waves and surface 
waves. In:,H:. Benioff; M. Ewing, B. F. Howell, Jr-.,-and 
Press (Editors), Contributions in Geophysics: In honor 
of Beno Gutenber. Pergamon, New York. pp. 1-16, 1958. 

Bath, M. and Duda, S. J. Barthquake volume, fault plane 
area, seismic energy, Strain, deformation and related 
quantities. Ann. Geof. 17, 353-368, 1964. 

Ben-Menahem, A., Earthquake similarity laws. Phys. Earth 
Planet: Sinter 1S .eP C=P 1S tore. 

Ben-Menahem, A., The role of shear mach number in earthquake 
Sourceudynamscow, -SullaesemsmarcOo. PAM. OG SITO) ati oo, 


RS:76". 


161 


Viet? Yeo Lomeiee2 aviwgdcr nang, oh Ang 2 


rho D> Dam Lema se a ie wn? Bate i 


a A 
Be tA 


8 


HOF DEiMSnvh“os sgl ge granite “leas ni ink ee 

jaer. . Bagh es 3 ad hc | ane aha. oan «aia 

etd euoyer 3 teed etie. a Lo has oe Oi arms 

270, 2hteede let dangle eae enol 

d- gite La np had He ROA 4 KPT yoke biG Wy oh “8 Bey i 
{ ent 
flow (erasott bao der debo iy ansel err tee aebtbak os 
O50) .CemeBCe eae ee ea 

snot ie Ste eo ree ens whieh ey. Re Pore oar 8 eg as 

SiG yet VER ONO: oR FY Oe ce de tees 3 «nt err i 

"Lorttost nt ‘ele Reee At eres) Such aOMOD ies Desay og auig "a 

$28). at A ae avhit, wm: nowmiggiel Jrsdnss0m conga, | 
Shaws (Sisek | amslow oxcope tee y We... \ shut Bae ie a 

5a tude? erx eit ene ath , Whee ayers  simelee er ‘ 

/ pegs. apr ie tien et cot ei aaa 

dote8 dyad .sihi ysoetans Aa Ires . «A ; 

ROE, SFSee 2) jawed tooett| ey 


Savers tee P ‘ete ier soarig - §6/ S207 a7? ..h dca haat 
eet =| a . ’ od « Dulk . soe ime 1p8 ¢hh05 g-«© ag te Ranh soim08* 
| . . ater" ay 


162 


Ben-Menahem, A., and S. J. Singh., Computation of models of 


elastic dislocation in the earth, In:B. A. Bolt 


(Editor), Methods in computational physics, Academic 


PRESS; sNews York, 12.) 299-3752 41972. 


Beaumont, C., Foreland Basins., Geophys. J. R. Astro. Soc. 
65pe29i-329. 1984. 


Berckhemer, H. Die Ausdhnung der bruchflache im Erdbebenherd 


und ihr Einfluss auf das seismische Wellenspektrum, 


Gerlands Beitr. z. Geophys., 71, 5-26, 1962. 


Berg, E., Relation between earthquake foreshocks, stress and 
mainshocks. Nature, 219, 1141-1143, 1968. 


Brace, W. F., and Kahlstedt, D. L., Limits on Lithospheric 


stress imposed by laboratory experiments. J. Geophys. 
Res. 85, 6248-6252, 1980. 


Brune, J. N. Seismic moment, seismicity, and rate of slip 


along major fault zones. J. Geophys. Res. 73, 777-784, 
1968. 


Brune,J.N., Tectonic stress and the spectra of seismic shear 


waves from earthquakes, J. Geophys. Res., 


Eo yn 9.9:7.55003 
Meio 70 


Brune, J. N., (1971) Correction, J. Geophys. Res. 76, pp. 
5002 


Brinewadedie, ALenuleta gRediepyand Haguezell 2 Se rarcfield 


S-wave spectra, corner frequencies, and pulse shapes, J. 


Geophyss (Ress 64 pn226252272 7. 1979 


Bell die Sen andGoudhoD apliepeNoreheastssouthwest compressive 


stress sin Alberta evidence from oil wells Earth 


ee 
ie eietem 26 oa a getionia, hy py 
pit mt 

tide. SA RE vdvaina pete. ck 


4, 7 va a ) r 
* rd i] pat | _ r 
simeheos .astakde Sanet sezugRes ak at RWAROPEE 


akg sPTE“ReE ‘te eee . 
eet Se ek ee eee nt 008) y (AOE SEK): a 5 


Rein ae ipaee 
MWajeoeliad wtSetaeied Ben de angitole as 
2). ADE Sg eee me Aaaton: 
iS seated wioeign yes stpopiliiised eet sien 7 
Bari > Cala ap Rib spread. wtos ia 
“ertacred 6c He es vee See .abasetiiod BAe tt es 
xno) |G) 220 02 TS. He aa 5 bapeqnd ‘2 ali, 
| POR Aatiabss ay a 
gite 20 (820% be. «ya veeiied sa) pairtietag Sgt oe 1H yy 8 
METH ee ae aN “P iis it token-oet , ‘| 


ee. gs as ces 4a ‘ai! 4) Ses time ea 7% fs ao RMS FT! 


sada Jimeiee 2s eiaiege cade tern ry) Peenctcals ed t 
‘ air 4 


ee ee «aitgoae 7 sotnnseve nord, seve a 
i AT pare 2 eeaean .\% ae bisedaed 4a Ter). vat o bye - rane 
\ : ies ame ; dopa 9) ma 
b.Saleo tes , a Ls goo pe . Dera. Ane gevetuhsiA ,. Mel? aa 
‘badgsituleniag, tytn’ ae ieine th aWkwe> SisssRe avewSs) 1 if 
| tes |, GYESHRORE. (ke ..co% eeieeeee 
ViRRergars It ewisy oe -egpessot). 4 6S doecd Bae 448 ae See, 


a ; he 


dase’ allew Lig moat enehivertesiedls ni abeaie 


G5 


Planet. Sci. Lett., 45, 475-482, 1979. 

Bokman, J. Post-Mississippian unconformity in western Canada 
basa@n: In Childs, 0. E., Editor, Backbone of ‘the 
Americas, Tulsa, Am. Assoc. Petroleum Geologists Mem. 2, 
2ONS242 4 1963% 

Burridge, R. The effect of sonic rupture velocity on the 
ratio of S to P corner frequencies. Bull. Seism. Soc. 
Ams: 65, 667567575 119755 

Burridge, R., and L. Knopoff. Body force equivalents for 
seismic dislocations. Bull. Seism. Soc. Am. 54, 
1875-1888. 7964< 

Bally A. W., Gordy P. L., and Stewart G. A. (1966) 
Structure, seismic data, and orogenic evolution of 
southern Canadian Rocky Mountains. Bull. Canadian 
Petroleum Geology Vol 14 pp. 337-381. 

Bakun, W. H. Bufe, C. G., and Stewart, R. M. Body-wave 
spectra of central California earthquakes. Bull. Seism. 
Soc. Am. 66, 363-384, 1976. 

Chinnery, M. A. Earthquake magnitude and source parameters, 
Bulls Seism.?'SocsVAmin 59 196954 96270 1960r 

Clowes R. M. and Kanasewich E. R., Seismic attenuation and 
the nature of reflecting horizons within the crust.J. 
Geophys. Res., 75, 6693-6705, 1970. 

Chaelwib © PF. andiuStewart,»G. Sv, “Recent slarge earthquakes 
along the Middle American Trench and Their implications 
for the subduction process. J. Geophys. Res. 87, 


829-836, sto 2. 


{ 


Lager an-ets ioe) 7 
Hk Aeteew Tih Saeeia cond ne fees 
si3 %o snoddeas pigs ibe “un di 


nit agatpodesd? ustines ree a BB 
ws Gin 


if 


QO Uthoolss tga assis ; 

“oe ‘toe: . Lae | Perr. ween so4 Sec 99 a sé 
| Terao. aan 
exchine 2968 Ros re oe ee: 
vA . 508) Joe opiria anol dapele it ‘ott 


; q Balt sae ie 
(Wet 4 ea eee) Eee i : 
id “Se bared polangint a) re} vedipin sia re co 
bie | 
cathe .TTUd sabeG cae ine besa ic mage | ? OR, ! 


TEES: Ree PN ee | 
ewawrehee 6 J Uo tenes Wale Ae ris ne Ht Bia 
niied . Laie. eae dae wt aaalh Sage oe io nite 
| ), eee yueentes BB: mh S | 
sv4cSactes sap iee shégiapen gdeepe t152 Tees ar ve erin tt 
Baer ,seet-eset ce ne!,ooe awed Beer ae 

aan volstrasts Tene kee 8 ,@ Mifweseted Sine oh 7) aie 

wert sat oltdew \gion Pod, pniseediie: i suudian ‘eda 
aver  abcastee, aoe . leet -edaoed! 

Seveuotdyhe e926) Thee a Ve jones 32 AS. ot teed a | 


eacicEestiqns ied? fngoltgnest seoRPegA slbbin ade pneks 
{8 .gaa .2edqae® 4 .aespege gelsicvbdas ad?) tt 
(iM Reet” BEEweeE tee 


164 


Chandra N. N., and Cumming G. L., Seismic refraction studies 
in Western Canada. Can. J. Earth Sci., 9, 1099-1109, 
FIT 

Cox, A. Plate tectonics and geomagnetics reversals. W. H. 
Freeman and Company, 1973. 

Cumming G. L. and Chandra N. N. (1975) Further studies of 
reflections from the deep crust in Southern Alberta., 
Canis (Earth Seve. i227  539-5573 

Das, S., and Aki, K., Fault plane with barriers: A versatile 
earthquake model, J. Geophys. Res. 82, 5658-5670, 1970. 

Dahlen, F.A., On the ratio of P-wave to S-wave corner 
frequencies for shallow earthquakes sources, Bull. 
Sevsme2Soc 7iAms,. 6&,.) hwS9=1. 150% x on974 

Dahlstrom C. D. A., Structural geology in the eastern margin 
of the Canadian Rocky Mountains. Bull. of Canadian 
Petroleum Geology 18, 332-406, 1970. 

Dahlstrom C. DA.) Balanced Cross sections. Can. J. Earth 
Sci. 6, 743-757, 1969. 

Douglas, B. M. and Ryall, A. Spectral characteristics and 
stress drop for microearthquakes near Fairview Peak, 
Nevada. J. Geophys. Res. 77,351-359, 1972. 

Elliott D., The motion of thrust faults. J. Geophys. Res. 
Steet 92 9-963) e576 

Enddan,. Buin .war linn Be. Aus, Gand Romney C.) FSevsmic. waves 
reflected from the earth's inner core. 

Ell isi ReNIM..,) Dragert): Hige andvOzard,- JU. Mis Sersmic 


activity in the McNaughton lake area, Canada., 


4 iy i “" 8 f} ms 7 DR: aii3k “4 y 


H.W velaaneven sakggnpsoov: oa 
‘ eres " 

ipiauga tei GELS. sa "avn ti 
Cat: perry HE) veo yaok sth mora Ines 
repece onl yoke a oT: 
vertqQ A saak teed tie suislG ahah ea iia) ne - 
aro; - @hae<Gend) Si 1 ae aoe “8 -SoBSon lite | 
boo » Lode st none te eee S. Se aitet ate \90 ets 44 
fut) nee poe elikeepelor eects LEAR 389 teins vi 
£784 | eo GENT aE Y 2B py om Rrcaita Hs | 
ean eee bh -a8 ae wpetose’s bie serie ick my ve 
pens Se EtoR sho setaiec nts 17 A, ont -am 
pine! Upaonee Bt ‘Poe la0! men Loe 
éycret ctolnedD Brel sae seca ala x yok ol ee 


wa Cat aoe, oa < 
ee | 


ae 


bis: sae setaggoe eto Taeomege AGaheye bas om a 
Abed peretist zeae soil ap omega te 702 nah 
258) | BRE~nde., eal vaqons anh, oI 
fae avdeend it au (ea rage Ad nod tom sat tie | 3 
| fe weve eae “en: aa . * 
iSvne S2Me Tae) FT 39 geoee bite am Ps tc io Oh eo 4 fribigad 
so. tanot Staak els ox bagseties | 
Dre iOk . “ i J22Ese ‘te ge pieepaxd y wit ahi 
abetted porte let Mogae 262-0 t ya tvs ee 


NGO 


Engineering Geology, 10, 227-238, 1976. 

Ellis, R. M., and Chandra, B., Seismicity in the mica 
reservoir (McNaughton lake) area: 1973-1978. Can. J. 
Earnth?Sci<) 18,4 1708-1716, 1984. 

Evernden, J. F., Study of regional seismicity and associated 
problems, Bull. Seism. Soc. Am. 60, 393-446, 1970. 

Fletcher, J. B. Spectra from high-dynamic range digital 
recordings of Oroville, California aftershocks and their 
source parameters. Bull. Seism. Soc. Am. 70, 735-755, 
1980. 

Guzman, J. E., and Zoltan de Cserna. Tectonic History of 
Mexico In Childs, O. E., Editors, Backbone of the 
Americas, Tulsa, Am. Assoc. Petroleum Geologist Mem. 2. 
pps ris i295) 96s. 

Garzaets,9Gilvd., Lomnitzec. ~~ Tubillar Ai, Estado de avance 
de la red sismologica Mexicana de apertura continental 
RESMAC. Ingenieria Vol XLVIII, Num 2 UNAM , Abril-Junio 
1978 Mexico. | 

Geller, R. J. Scaling relations for earthquake source 
Parameters and magnitudes. Bull. Seism. Soc. Am. 66, 

1 Oi O27, ea Oe. GS. 

Gibowicz, S. J. Variation of source properties: The 
Inangahua, New Zeland, aftershocks of 1968. Bull. Seism. 
Soc. GAMemoo ms 262/06; 2519 75. 

Gutenber, B. and Richter. C. F. Barthquake magnitude, 
intensity, energy, and acceleration. Bull. Seism. Soc. 


AM oy 63> 19 ise 19425 


| Tee RES Mas oe i: E 
fp rf wits ephpien t ac rare: ae’ shaw * 

% sce .2Ce Ree aibete) (adel Re ripest ae el 
ee 

tal aosRs nih, fa t+ ino dee Dheorsar Fe ba nme ’ 
Te ee ah qamale 
fabsoid epites ee va ment wa ee a 
73 Ste eSo0de3e¢ttn Aero niiad iat: | 


PAV SET Ae A eR) Re i ee eer ater 


Lewes ct sieoeset eS WORD wy veties Bap ak or 

ade He avtiioed (OReFhe Qa ct BALERS, iy 08 

A eet Sekpeiow? whelgyeeo (o928h 38 uel z =| 

oS cae 3 

niet, SB oBesas . v& bbibdicies2 a: ytd, Aa tah Qe) ‘ane ko se 

favour Hija, ehitiss Sh or Ble eatoatane he bat ie 
pauGeLiiak , MANU Rani CARO la ee Meme cates 

A ee Lob bestt, seer) 


SA5068 stuugeitige tH ened sate patios ot hue 


e 


Pe -«s soe. wakes . bias aehud kage bis. meson) 
; ot Gi chat nuhgaer | a 
: ent ‘eae: srayouy satis’ 72> petdel ‘we solvation 4 
mies shlga. ae! 26 golodtewesl4 {ensles wel sudegnan® he 7 
OVS) yet ot), 29 oe suealelt AY 

jabusinpam adden 3a 4b so. seeenik Sea oe . tedhetue, 


30% sie? -List . onde teletes) Gok ,Yousde.. yt bandana 


SG). 2' Sat SE Ty 


166 


Gutenberg, B. and Richter, C. F. Earthquake magnitude, 
intensity, energy, and acceleration (second paper). 
Bull. Seism. Soc. Am. 46, 105-145, 1956. 

Gough D. I., and Bell J. S., Stress orientations from oil 
well fractures in Alberta and Texas. Canadian Journal of 
Earth Sciences 18, 638-645, 1981. 

Ganley D. C. and Cumming G. L., A seismic reflection model 
of the crust near Edmonton, Alberta. Can. J. Earth Sci. 
PROSTOTL=1O09R' 1974. 

Geller, J. R. and Mueller, C. S. Four similar earthquakes in 
central California. Geophysical Research Letters 7, 
S21=8 247 +1980" 

Garland, G. D., and Burwash, R. A., Geophysical and 
Petrological Study of Precambrian of Central Alberta, 
Canada. Bull. Am. Association of Petroleum Geologists 
£3 4790=806; S959). 

Haskell, N. A. Radiation pattern of surface waves from point 
sources ina multilayered medium. Bull. Seism. Soc. Am. 
54 (377-394 201964: 

Hanks, T. C. and Thatcher, W. A graphical representation of 
seismic source parameters. J. Geophys. Res. 77, 
439354205), 4097 2 

Hanks, T2°C.; “and Wyss, M., The use of body-wave spectra in 
the determination of seismic~-source parameters. Bull. 
Seism. “Soc. *Ameft62, §561=5907 19725 

Hartzell Ss. H¥eand=Brune, od. *source ‘parameters for the 


January 1975 Brawley-Imperial valley earthquake swarm. 


,etwi ingen a pt 1% ft wi sei 
« (Seam, CRO Seee hetinieteace a “ | ae 
avy Nid +20: ae oe aa 
Lionee we, ta) Bhd WROIIR 9 tke + an ® r 
fame) e&4 teeae? eenen ore eared faint | 
(hae ,tearsia’ ae tegcie a3 

ae y Hoc 4Oe bes ~heneeeie 5 tue onlieedatiled oo oye | 
i+ te50k ,f cee Vee is scent) sane oe io 
vita, enF-1 ORR 
‘aaron ast ues Wer at , ag be some Mas i a 
ee tae | ADT AeA, ‘fac iaiyiqon® .ekaati tes ie thes 
O88, setts : 7 
ey tele ip sitesi Poh ORY Seana, ris She ae a 
wi kt) Paltees te. erect senpe ie yn aes 
ren heals sa Pe Believe +a Lee , Bet 
(leeey s9H- 0008 | . 


ay 


De hed Maro +e is -ueaalipiatbes o& 14) 


tite 


Tots aS 
rt “= | @e Dek, wie Ae besgyetie tus & nz Bee . 


hh HIRT teh ee 
jo mmlisiassenges Taniigase A .8 0 aia Hite Ror ane ie 


oe mi EE 


we avn R ayetquen: vw 2a areh S340 % Sages: Siam tee 


(ei \ Ce ees anna-tebaldiael 
rh autosgkh eave (hod de San ee | Ree Gre 62) at vedas 


; 
tive .2*useusiee saedecep bee GQ ools+enlwiesdbcen? iy: 
R09 DRE-D5R yg Re .ak 1900 meee: 
eis wot gxueteneteg osveR.h YRRESE Ons .koo 8) heeesiee 


-WIsve Staepdiian es.iay iat requuk-geiwese SORf ypeutiate oy 


_ ? 


167 


Pure and Applied Geophysics. 115, 333-355, 1977. 

Hanks, T. C. The corner frequency shift, earthquake source 
models, and Q. Bull. Seism. Soc. Am. 71, 597-612, 1981. 

Hanks, T. C., and Raleigh, C. B., The conference on 
magnitude of deviatoric stresses in the earth's crust 
and uppermost mantle, J. Geophys. Res. 85, 6083-6085, 
4:98.03 

Horner, R. B., Stevens, A. E., and Hasegawa, H. S., The 
Bengough, Saskatchewan, earthquake of July 26, 1972. 
Gan SeJ nebartneScl 1 1025 1805-182:1.,. 197s. 

Iida, K., and Aki, K., Seismic source time function of 
propagating longitudinal shear cracks, J. Geophys. Res. 
27, 2034-2044, 1972. 

IMSL 1979. International Mathematical and Statistical 
Libraries Inc., 6th floor, GNB Building, 7500 Bellaire 
Blvdse, «Houston; iTX.77036, UsS.A. 

Johnson, L. R., Seismic source theory. Rev. Geophys. Space 
Phys. 17, 328-336, 1979. 

Johnson, L. R. and McEvilly, T. V. Near-field observations 
and source parameters of central California earthquakes. 
Bull. Seism. Soc. Am. 64, 1855-1886, 1974. 

Kanamori, H., and Anderson, D. L., (1975). Theoretical basis 
of some empirical relations in seismology. Bull. Seism. 
Soc. “of “America 752 pps 107323095 

Kanamori, H. The energy release in great earthquakes. J. 
Geophys. Res. 8274298 142987, shohs: 


Kanaswich “Ew ~RapcClowes@GR. Mi sand {McCloughan),*G¢. 0H. A 


soehe jphoeile cin ian te 
eves. Sdetyid eRe hide. eae 
ae) Brates aR oY sh wigt ee iteaas 
eo 3 a ale Big phe ot debstan t as 

sare’ eee Ch ie ae erie 132 a 
‘odaeteds (38h). cat aeegee ..t Sn | 


ie... — 


re 

SOT wie «i aint oeaeh Rita pa, Ai en ed A =n 

POAT Se Oe ae odipunitaee newsdagades2 re shh. re y 

Hay abe ea, ‘Or ae tat tite ‘3 

Longe doen Si: OS bingiae uc mba. ona 

ook aviGcan .b ,2eoe8s ar “at bo $ tone Nena: oa 

p uinbh et digt-acae 4 
favdusid ks Oi that simaatiatt ‘tapotaacasinn, LER 

siteiies® O8et , wwis the: RING: yeah Adds abn “pelt 4 


ttt Ts 


Las dbaty ab inetenait «i Ea psi 

éoet agdgods. .tet  <Qtoans etait sine tee on ar a ont 

| 58 Laet-ase ivf, “ 

ano fdevessad bisty+ tee 7 vbtivaps bee oi ide 

Seance sae ets FLD, fasoras. 26 baetades oq a ee onc | 
bu2T sORhe nGOBt wy ti oer Taste il cor ae 

reped Lei dsinantt paeee be eee ticle a 4. DAB ps pee 

mete? sEAge, “poloemior ns aie ldetae bat tique ones S950 

NEC <2 Wea gel sivas tel soee ae 


4 


Jageucat+es Jkeke Al Stee ed, giaene odT 8 |. 2tomenan 
TkEL  faeAPGGg,, Sh .aet- ayddoss 


Suet) 2) dedgonloa Geo ot se ROO? , 8) A de spaenen 


168 


buried Precambrian rift in Western Canada. 
Tectonophysics 8, 513-527, 1969. 

Keilis-Borok, V. I. An estimation of the displacement in an 
earthquake source and of source dimensions. Ann. Geofis. 
$2) 205-24 959% 

Keating, L. F., Exploration in the Canadian Rockies and 
South West Alberta, Canadian Journal of Earth Sciences 
35, th 3- 7235 e0 066% 

Knopoff, L. Energy release in earthquakes. J. Geophys. Res. 
ip? 24-5266 1958. 

Lomnitz, C., Gil, J., RESMAC, the new Mexican seismic array. 
EOS) 5 /Otec, 1 '68-09;6 | S7'6e 

Maruyama, T. On the force equivalent of dynamic elastic 
dislocations with reference to the earthquake mechanism. 
Bull. of the Earthquake Research Institute, Tokyo 
University 48, 467-486, 1963. 

Madariaga, R., Dynamics of an expanding circular fault, 
Bull. Seism. Soc. Am., 66, 639-666 , 1976 

Madariaga, R., Implications of stress-drop models of 
earthquakes for the inversion of stress drop from 
seismic observations, Pure Appl. Geopys., 115(112), 
30g 34.62, 3 Sol 

Monger, J.) W.cH.s andi Price,’R. A.; Geodynamic evolution ‘of 
the Canadian Cordillera- progress and problems, Can. J. 
Earths Scae 678 7077/94) 1a 979° 

Mi Ines) WeG?,nand) Rogers. RiP.) sRiddihough, «McMechan, “CG. 


A. ,.and Hyndman, R. D. , Seismicity of western Canada. 


aera Asi. nih cs 

{ “baer reg an be 
Ut Yasmosetqans sda ao? al teks i YY ee 
oad. @he. . wre se lbiboruce- ao fini : 


avo 


espyatt® atzaad To cast rie Po aD esata ae $e) i | 
suituaap 7%: shea a3 romp te, vanes a . 


Lota “De eee heodart emt wd’ ,pAMEES m4 ly ihe o.. 

pmtet veaebay, ds ai 

fae Drs heey (oO) Ape si A gard ied nG 

orn pret te ebsites | Wen NA ote pet * y ; 

| Sveo? . ssn dank Pogo ont hcheanaeal aa? 30°. ‘i 

| PRED) quire he oe Le wen 

Sivet ten o te rapes ie Le by sanenya ie), 4 ar 

| STE hy Baa+ees aa" re gee a2 Fee a na iy 

io @tahem gave-wbedas To . grabs tognt , eye 

meet? yee gee7s0. ta 5 satus "ta onsnigitteibal aa ' 

Aor PY; Sepa Rak, ate dato seed blmetee ’ 

: \ | 7 ht ar e-1 bee 

i Pb eeee, line eeRD aM), a eRe ee. 4 GN eG pa 
ioe) jemeticoig ta seedy Saaeghtb105 wie aa aes 
EEE PETAUNT 6h). ER ae 

@atdpeeet (dover) Sham 4.4) Seem had’ . AWG eRcees 


4 


<Sg3040 ogei few td BFiniseleset  ..G ie . nanbart S76. ek 


169 


Canwd. thacth*Scry (Sige 1195, 1976 

Milne, W. G., The Snipe lake, Alberta earthquake of March 8, 
IS7O2 Canyred. ‘Barth "Sere 7, °1564-1567 "1970. 

Milne, W. G., and White, W. R. H., A seismic investigation 
of mine "bumps" in the Crowsnest Pass coal field. The 
Canadian Mining and Metallurgical Bulletin, 51, 678-685, 
71956.. 

Molnar, Pe, Tucker, Bob .; and Brune “JIN. , Corner frequencies 
of P and S waves and models of earthquake sources, Bull. 
SETS. .750C. TAM. 7 OS. "tooo Ntou. "tors 

McGarr, A., Spottiswoode, S.M., Gay, N.C., and Ortlepp, 
W.D., Observations relevant to seismic driving stress, 
stress drop, and efficiency, J. Geophys Res., 84, 
Doi ec Ot hore 

Munguia,L., J. Brune, A. Reyes, J. Gonzalez, R. Simons, F. 
Vernon, 1979, Digital Seismic Event Recorder Records and 
Spectra for Aftershocks of the November 28, 1978, Oaxaca 
Earthquake, Geofisica Internacional, Voom os 11 oro 

McCloughan C. H. and Kanasewich E. R. (1974). Geophysical 
data link manual. University of Alberta. 

McGarr, A., Some constraints on levels of shear stress in 
the crust from observations and theory. J. Geophys. Res. 
Soy 6231-6250, 1J0U. 

McKenzie, D., and Jarvis, G., The conversion of heat into 
mechanical work by mantle convection. J. Geophys. Res. 
854.6093 609," hooU. 


Martin L. J., Tectonics of Northern Cordillera in Canada: in 


1) qa year} ~0r1? “ ‘ | 

5 Ansa, 36 we apd, @daids rear ia | ; 
ee ee ee ee 
ang iatacvnl. Shai 1 Se re wag i i 


'e .atteite® Tagegie lian ee’ 
a. ibe 


Dp 
- 
> 


rie , lity aitaed, bik ies * piu 

eharght Tee o eheten bane! iaenw & ong | 
| rey tt we ‘eat + sath ‘ate ioe 
ageLero rik cM ERP . ee Abou abana , a 
etatze@ ORiViTb nla hal oF chat ee ee 
(8S “eer eyigaed uty Hered 2th maa, mesh es : oer 
| Ais ee (aSeag ‘ 
SER A) eal I hy eat yh, anes: ee 4 yee gt 
sp abot “seb roses sagen 2 Lirdiee apie: eet — Vv" 
Seun0 , aV8? Bay seam e add 24 sitdonerd sts 10k 
hha Unt Daw pieaweie se une vet et otaopio 

Lenbayngowo tage Poa ow eS ian wii 8 9 eats eS 
titre tha te’ tie vo etl Lita dort pa Hs 


. 
ms ~~ 


Sl nieise tome ae diwwes se ainbemreinos stn? ye 80M 
mK, i 
ah 


Ever 


tom cegdaeed ye waive Ses eiodthwteace werd +auTS ? ‘ 
O88! cat Sih 
édli sadid to nolstevnes ae a vbbviet. fetal Wi 
268 .eyharc® .% .dolensuned ci ya Sade Lesh 


waver BEOR=208 
rh ameaiiee of ‘aieliioses raat ‘apm tie 


ce 


170 


Childs, 0. E., Editor, Backbone of the Americas, Tulsa, 
Am. Assoc. Petroleum Geologists Mem. 2, 243-251, 1963. 

Mogi, K., Study of the elastic shocks caused by the fracture 
of heterogeneous materials and its relation to 
earthquake phenomena, Bull. Earth Res. Inst. 40, 

P25 M7 34 1962. 

Miyatake, T., Numerical simulations of earthquake source 
process by a three-dimensional crack model. Part I. 
Rupture process..J. Phys. Earth 28, 565-598, 1980a.. 

Miyatake, T., Numerical simulations of earthquake source 
process by a three-dimensional crack model. Part II. 
Seismic waves and spectrum. J. Phys. Earth. 28, 599-616, 
1980b. 

Marion, G. E. and Long, L. T. Microearthquake spectra in the 
Southeastern United States. Bull. Seis. Soc. Am. 70, 
1037-1054... 1980. 

McNutt, M. Implications of Regional Gravity for State of 
Stress in the Bart's Crust and Upper Mantle. J. Geophys. 
Res#.83,'637/-6396, 71980" 

Nuttli, O., The effect of the Earth's surface on the S wave 
particle. motion,-Buld. | Seism. Soc... Am. bol gecs 7 a2e6n § 
1961 

Nyland, E. (1981). Stress and strain in the North America 
Plate, Abstract .G. A. C./ C. G..U. Meeting May 1981, 
Calgary; 

Orowan, E. Mechanism of seismic faulting in rock 


deformation, A symposium. Geol. Soc. Amer. Mem. 79, 


‘erhy -Seors ene | ont Me Paodione 
fier. Tes eee yas ee Gbikieneat | 
(JRA os yo Deegea Pijpacts hsnsel HR 
ion HOE je haet 1, Bas ont vatieme 


o Ge. ‘tan 1 ie “Tidal ob Last eae 


yoe adsedilcina Sy ewe had uate ha dat ae a 

2385 Lp a, ne as s iNe? sage Repeal 24 ae 

BOBS) Sta Aud, (AR eae aOR ee Eaaaaaal rh 
SITVOE, a SUH ae te ena.crel i umbe S34 vam eae 
Laps Pa Sou Aes or Seanemtbr pod, taba eee 
ita-eee ei Pcs i Pa.) Ges mate sede a ae sa 
BLS oMeR, od 26 ohh ad a h. =) ot bas, “ni i 19 


Oo ed sa62) \etae Bey , lnenicirn Dine sotarttae. 


ns 


5 “4 
Sey 


wicse 36d eo Pvass. tenet geet. % bapieno: 


‘is ea i et YTEe 


sien & oft a anetsina2 ‘ees wag te, angites ett oe 


aoe 


(ASO PERG Ey 4 A) OR RRS Aye. snot tom Lola 
\ ; avs 


Ag i i ni ie 7 Ag 


BciCURA agsal wag At Abeoase bag BewIz2 , fr 0et a ne atv n 


epee ot 
bes 
- re 7 = 


Ve ePowi's eM sO atl Nye ot sail toeitada (oaeke a 
eae op: - aut 

toon. cl, pebeive? id 1s) me Lad oeM, 2 

| ue 

12S oR TOR OR. eran ice ha f aodinaig te 


oo mae 7 ots a yan Bia 


171 


2327545), '960. 

Ohtake, M., Matumoto, T., and Latham, G.V., Seismicity gap 
near Oaxaca, southern Mexico, as a probable precursor to 
a large earthquake, Pure. Appl. Geophys., 115, 375-385, 
TOW, 

Ohnaka, M. Earthquake-source parameters related to 
magnitude. Geophys. J. R. Astr. Soc. 55, 45-66, 1978. 

Pilant, W. L. Elastic waves in the earth. Elsevier 
Scientific Company, 1979. 

Pearson, C., The relationship between microseismicity and 
high pore pressure during hydraulic stimulation 
experiments in low permeability granitic rocks. J. 
Geophys. Res. 86, 7855-7864, 19812 

Ponce, L.,McNally, K., Sumin de portilla, V., Gonzalez, J., 
Del’ Castillo; A...) Gonzalez, L.~ Chaels, E.* French, M. 
Oaxaca, Mexico earthquake of 29 November 1978: A 
preliminary report on spatio-temporal pattern of 
preceding seismic activity in mainshock relocation. 
Georis.- Interv 7 225  TO9R1265 Fei ee 

Prince, J., H. Rodriguez, E. Z. Jawaski, G. Kilander, A 
strong Motion Radio Telemetry Network, Proceedings of 
the fifth World Conference on Earthquake Engineering, 1, 
TODS 1M Os OSs, 

Randall, M. J. The spectral theory of seismic sources. Bull. 
Seism. SOC.” AM. O57, loon lee, lass. 

Randall)” M.°U.)" Stress drop and the ratio of seismic energy 


to moment. Journal of Geophys Res. 77, 969-970, 1972. 


cae Wotoime let, Pe meee Sire yy 
yoajwoetg VIARIRiy ae i the ate fa 
Sag-ave (yet! . ae ea wind yam | 
te ' mi : . : F 


f 


os beaslay GreteMbsey po TR 

cat (sarees ae “va OA LS ’ 
“event. tees ads a fon ae ae cr 
Saat Ra mn 


lean £620 7% ret AaRe. Maal 1 Ee daeasbibaiemea 4 ee ¢ 


i 
7 
‘e 
af 

a2 
= 
% 
-*- 
, 


i gtte> ois baw Ea 
Peet y" inside ‘ad sani a 
eC oy Se bc a ort leven, 4 
M uteneth: 2 totte aA heat bam vaciiiaelt: 
vet edaneyew ey 3 apa BS ikea 6 ae 
ested) Lee yk , Reg Sia pages yrenime ye 
de JO L279 ROE GSE TFL a Sik ofhihaie eiates enibsvesd | 
| SBP RT) ad) « Oe oe . 1etat sie 
&  tagtel en Ce hapewey-we 8 VRAD sisost ‘H vot 4 | 
aeeTh beet s9 epee usu fo% pibart nay 2oM enerte “s 
ppisacai iowa waded me My songtetnod bisae ty oe baa 
ere jeonr-aaene 
ilu eecivee ataaiaes ss ‘acids tatépog2: sat GD) Mw joeninaa 
la} (OM REEPE 1,2). .24 208 Jan ten 

sotens sifcie® to abnwe eda) Sle CRePiedante Yi oe DAR 
Stir Oc @-eee 7. atl titgoat Be Seni. Samo | St 


“= 


Ie 


Reyes A., Gonzalez J., Munguia L., Nava A., Vernon F., Brune 
J. N., Locations of aftershocks of the Oaxaca earthquake 
using amoked paper and digital event recorders. 
Technical report Sept. 1979 

Reichle, M. S., A seismological study of the Gulf of 
California: Sonobuoy and teleseismic observations, and 
tectonic implications, PhD Tesis University of 
California San Diego, 1975. 

Richter, C. F., An instrumental earthquake scale, Bull. 
Seism. Soc. Am. 25, 1-32, 1935. 

Richter, C. F. Elementary Seismology W. H. Freeman and 
Company. 1958. 

Richards, T. C., and Walker, D. J., Measurements of the 
thickness of the earth's crust in the Albertan plains of 
Western Canada. Geophysics Vol. XXIV, 262-284, 1959. 

Rudnicki, J: .W: and Kanamori; H.» Effects of: fault 
interaction on moment, stress drop, and strain energy 
release. J. Geophys. Res. 86, 1785-1793, 1981. 

Rebollar, C. J., and Nyland, E., Spectra of some Oaxaca 
earthquake aftershocks from RESMAC. Geofisica 
Internacional 19, 109-127, 1980. 

Rebouwtar ec. J0+,eKanasevichy@n. iF. -andeNyiand,, B.. Focal 
depths and source parameters of the Rocky Mountain House 
earthquake swarm from digital data ar Edmonton. 
Submitted for publication to the Canadian Journal of 
Earth Sciences. 1982. 


Reboldlare,. Cine Spekanasewich) G2, R.y7 and: Nyland, E. Source 


nt ‘ af ] Ce 
ike 4 
ws et Hories 1s weit a Saal ‘ai ¥ 


supdioee soak aug te Rano s 1ad ih om 
Srebsocdd ANT had bo PB, teal “ 

gir 340 i) + depivizet 

4 Hi te ay hie i ; Pe pt 

ho Siup ete Be Nhade les Lee Tero ; 


14 i 
Ae Ces 


cae 


bog. ~Snar soy ei ci VMAS SASL SS He rie 
lo utherew toe g2ae? Gud Su if 
= aan aye a E70 Rh, 
LTRE ee es inte ih ae 
RB. Sant came pee am * 

fad. ehamie se UW eT eh pe 
ate - 
4ita 36 aandcrdastt ect 16, eer a! ues +5) *) 
zitiele, ped ies lh Ses 92, ceete Sales aint FS ae: | 
ae). Sears” Vinee See te sabes anata 
hiteee a BI Ste el twauben wat Sat MAG. ee : 
VRashke weecie Bnei Oe. see 7ee jaaences wy ae 
ihe OT Teaayy rae ek 5 emirate - 9 2 a 
PSA, SMCS TG eM peNge \ ye il Se + ah ’ 261 a i 


Get: ‘ae Gt lanot 

Laat, 5 \pogigt bah yA. aniei ger 7.0 oe 

MOK: Af ear frees ahd 1 evo Someta STI Te oo me i +a 

OCOTE 1h GA JOSS Pr ReMtOeT wes ava a 

to Larsvet we leaned wd? ag itv 2 wn Lidar 103 bets innigk’ | u 
—pebOEl ~nacnetsk, at rei? 

So7gee 8 Brass Bow 6A SA AMMMMPGA vol AD sabtodaa 


= Me ; 4 
j 


PS 


parameters from shallow events in the Rocky Mountain 
House earthquake swarm. In Press Canadian Journal of 
Earth Sciences. 1982. 

Rogers, G. C., Ellis, R. M., and Hasegawa, H. S., The 
McNaughton lake earthquake of May 14, 1978. Bull. Seism. 
Sock. Am.iui0, i727 1-7 8628 19802 

Rogers, G. C. and Ellis, R. M., The eastern British Columbia 
earthquake of February 4, 1918., Can. J. Earth Sci. 16, 
146 4=4493557979:. 

Rogers, G. C., The McNaughton lake seismicity- more evidence 
for an Anahim hotspot?., Can. J. Earth Sci. 18, 826-828, 
HOG: ies 

Savage, J.C., Relation between P-and S-wave corner 
frequencies in the seismic spectrum, Bull. Seim. Soc. 
Amt 645 k62i-Hi627 weal 972 

Saito, T. and T. Hirasawa, Body wave spectra from 
propagating shear cracks. Journal of Physics of the 
Barth: 20 41524315 °19 73: 

Sangh,» S.0K.., Astazjyl) and: Havskovisgder. Seasmuceigaps: and 
recurrence periods of large earthquakes along the 
Mexican subduction zone: A reexamination. Bull Seism. 
Soc. (AM... (ere Gev mostsy. 190.1% 

Scholz, C., The frequency-magnitude relation of 
microfracturing in rock and its relationto earthquakes, 
Buld. Seas: Soc.nAm.),58,- 399-445. qi968- 

Stevens, A. E., Some twentieth-Century earthquakes., 


Geoscience Canada, 4, 41-45, 1977. 


nistiarem vinad wk ee Bine79 ol 
iS Leewsct co Mage end ni ae 


: a 

wr, 
Wa 6 Sie 
fi 


= 


48 \eMegea One . Wht yee 
ie ‘i ce 


Nae fea | 
ties. ote) ee yaw ho ocsepatise emmE: 


* Ve 
~ 
}G 
afa 
tT | a 


\ we 
‘Sy 7 \ Le P 
4 F; a’ os, 
“ ~- : 7 : > ON - 
* es AGA sin uti pene? 
sLphrdas Sing pine ere tra padt 


nee... hus... pepaeea¢g lina tite ai euis 


eee: SbF ts | 
ARP are es disokae a (Ht hae 4 ii iit | 

CIARA xh ee: es a Be Cerin Jpas " sad 
node jeeats at | 
sso -Woin eee 4 ee oat Raw & ia . 
Ot: 23h gece PROBA Mi tert slant na | 
. ae pet 
ened ere Py) ee had pa negates +f bint 


t 
iy 


ey ° VS pre esary A" 

mgd det ome: aie dena Bhs = 
1o Eaieeas, de nore .8¥gniin gabde erkdiy 
| Be is, an 


aeeeet IA se a 


. 29d) piw- Bet 4 — 


; 
’ = . a 
, Beh eu param aeneera <<: hole ¥ 
» or i 7 rn ¥ 


174 


Stein, S., An earthquake swarm on the Chagos-Laccadive ridge 
and its tectonics implications, R. Astron. Soc., 
Geophys oo .5'S5i) 577-58 S8emes7 Ss. 

Scholz, C., The frequency-magnitude relation of 
microfracturing in rock and its relation to earthquakes, 
Bull. Seism. Soc. Am. 58, 399-415, 1968. 

solomon,” S.«C. ,'Riehardson, R. MM. andsBeroman, “2. A.; 
Tectonic stress: Models and magnitudes. J. Geophys. Res. 
85)'6086-6092, “980. 

Scholz, C. H., Shear heating and the state of state of 
Stress on faults, J. Geophys Res. 85, 6174-6184, 1980. 

Zoback?, uM. Di, Tsukanhara;oH.; and Hickman, So) Stress 
measurments at depth in the vicinity of the San Andreas 
fault: Implications for the magnitude of shear stress at 
depth. J. Geophys. Res. 85, 6157-6173, 1980. 

Zoback, M. L., and Zoback, M., State of stress in the 
conterminous United States. J. Geophys Res. 85, 
Ge 6) 568 a9 60r 

Stewart, G.S., and Chael, E. Source mechanism of the 
November 29, 1978, Oaxaca, Mexico earthquake - A large 
Simple event, Earthquakes Notes, 49, 47 , 1978 

Singh, S. K., and Havskov, J. A Moment-Magnitude scale. 
Bulll'Seism:s Soc. Ame7 7059 379-3835" 1980. 

Starrs Au TA Slapiinva cystalvand (ruprurcuarm cesolacsaue to 
shear. Proc. Camb. Phil. Soc. 24, 489-500, 1928. 

Spottiswode, S. M. and McGarr. A. Source parameters of 


tremors in a deep-level gold mine. Bull. Seism. Soc. Am. 


aie 


Ex surbetoutetiadhmll! We cones 


es: croc a pent sechtai 
om ees 


eiRURA tee oF WEES ne hi 
goer. Jey yo ha nee oa ee | | 


past ah a oa 


- 
—, 


ine oeelgee 6 .ae Ree Ga ed 


Ae aa ae 
eek Bean. aah 
| / Ge Sena seta cee 


Ber aR RNa ee aa diese) ‘he sand rie’ 
arth af We 2 Himdoee ai fs aH aT sdgduet re iy 

eed 3 Sel Oto teae4 eines Sif onty nt date 90 we ati 
yette Weds Ye ehoe a Hehe “fo% sah aon ? a | 
ey EN eerie $9 

ens te See7 7S (he eS Boe 

28) 2k wreies |b swags tone enh hich 

| 1 aes Bah a-ERE 


Yaos : 4 i L4 Br 


2 
, y's 


oix te tethaanesia odatind ie i 


$276. A - snideet ud Oca 8 
mer cin Sie wa se oi, at 


qira Bes a a *@ 
Pre: 
wa S08 a pea’ 


ey ay z ; 
Sho Gite’ s'nt eet bt 
. jie ah 


i‘ a eae sonal Day: 
amy , 


mae 

fee) ,OS~ees 29 ot ¥ oy 3 
) , 

du) Haro senesen oie a att saat 


ne) ee «348 ited an 


ee ee rie 


Li 


Cope sa i 291 975) 

Shamina, O. G., Pavlov, A. A., and Strzhkov, S. A. Model 
Studies of shear displacement along a pre-existing 
fault. Pure and Applied Geophysics 116, 900-912, 1978. 

Savage, J. C. Relation between P- and S-wave corner 
frequencies in the spectrum. Bull. Seism. Soc. Am. 64, 
MOZ1= 16.2575" “1 Si7-83 

Savage, J. C. Relation of corner frequency to fault 
dimensions. J. Geophys. Res. 77, 3788-3795, 1972. 

Sato, R. (1979). Theoretical basis on relationships between 
focal parameters and earthquake magnitude. oi Phys. 
Earth, 2% Op. 300-3725 

Steketee, J. A. Some geophysical applications of the elastic 
theory of dislocations. Canadian Journal of Physics 36, 
MV6G— PIO yes 956", 

Sahay, P. N., Seismic source theory. Master of Science 
Thesis University of Alberta Canada. 1980. 

Shaw E. W., Canadian Rockies-Orientation in time and space: 
in Childs, O. E., Editor, Backbone of the Americas, 
Tulsa, Am. Assoc. Petroleum Geologist Mem. 2, 231-242, 
$963". 

Sykes, L. R., Barthquake swarms and sea-floor spreading. J. 
Geophys. Res., 75, 6598-6611, 1970. 

Suzuka, PAPIStatngt ical study Of ties occurrences otesmal:: 
earthquakes (Fourth paper). Tohoku Daigaku, Sendai Japan 
Science Reports Fifth Series Geophysics. 11, 10-54, 


(i. 


ol 


he 


iat ’.4 18 eaagst ‘bap yok ea ; fo op 

pnlsal-e- 018 Gees snaiaoazaia kip tm | : 

ave a ie goreeEi esd bare ml Na 
ev gets Bite i ndamr sei cre a | 
me pet ata on mss age wi ict uae 
he |  PRBE a? cage a 
[vet o¢ dodegpest Yeange Bye ste 
‘ove -aae lat pet caali eee 46) Eods 
susad. ea liadytseden ie band capersadadtt: petal f 
yea .t . oho dt ipsa oer hai eta ed | 
Ses8 RR re 
viv qth to aoktieedehgen ipabagriqewy, pwe, +h “i ! 
whieh to teeteow Retbenad urcueunuteerhl Pris 
Re rysetn., Sys : 
mrad paz: 20 gese—eM »ysoed3 02 thr, aranien thy ce 
O88! pebane apene ia eritcraws aves is 


Lanwga Beasties a ats eviews iO weed Fachouass ee oa Mal 

yas aks 50, ‘te! e¢ndsond Tbe ae ae penn’: i‘ 

(PETES 2 MT Vel esters ane seaaeN soemea/iMA wantin 1 

{ fet Hae he nok 3” “pha PS Be we! osaupd i308 oa tenn” 

ote (0) Oe, @et ..a0h (2ydebed 

(jena 2c 40asttento ede we gbute) hap iin ieee A, 63 Qe 
cacG igheat .uiepie® vsade? +. (teqQug disve?) jeadeupddtes) ; 

Sg-0 \°) .apiagdqued sadeemiet’s ec soged agmeriwio 7’) 

= | er 4 


i Pa 


‘76 


Thatcher, W. and Hanks, T.C., Source parameters of Southern 
California earthquakes, J. Geophys. Res., 78, 8547-8576 
en973 

Tucker, E.B., Brune, J.N., Source mechanism and Mb-Ms 
analysis of aftershocks of the San Fernando earthquake, 
Geophys... iJ. Roy... cAstr 2§Soc #1549)0 371-426 2, B19 707. 

Tsuboi, C. Earthquake energy, earthquake volume, aftershock 
area, and strengh of the earth's crust. Journal of 
Physics of the Earth 4, 63-66, 1956. 

Thompson R. I. (1979). A structural interpretation across 
part of the northern Rocky Mountains, British Columbia, 
Canada Can £ is Earth 'Scir i16,. (pp .sel228-12451. 

Utsu, T. and Seki, A. A relation between the area of the 
aftershock region and the energy of main shock, Zisin, 
7, 233-240, (in Japanese). 

Vargas M.C SAS Minsona 4A. J fanduSinghank.. Paroiveld 
displacement from a rectangular fault with two rupture 
velocities, Pure. Appl. Geophys. 118, 991-1006, 1980. 

Vinogradov, S. D. Experimental observations of elastic wave 
radiation characteristics from tensile cracks and 
pre-existing shear faults. Pure and Applied Geophysics 
PINS ye MoyelsFaneheloh ge Whi sie 

Weertman, J. Continuum distributions of dislocations on 
faults with finite friction. Bull. Seism. Soc. Am. 54, 
(0s5 056 7arises:. 

Whitham K. Milne W. G. and Smith E. T. (1970). The new 


Seimic zoning map for Canada 1970 edition. The Canadian 


a a ex 


izengtoe Fo: cis temaisg: epauo? “<2 
Vaan sed”, 65 , eae abd. | 


Raa pes 


aI OMG ae Se ue) 


i be : 
shacpa tien =u Spee ae ee: Dima! 


PR BGP VE Ee VOR 
Mie@Istis “Sm Liev sAbioeN sss do punts ih 
So Zawpieal . Jevngt eae sid 
SOE RTS Gar GR. Eby 
S262 8: ACO GES. 20° ‘oid Sera auth 
sido (ad, Helo csi eae Pek ah 
bo 1 ~ OS) aes ie eam ah) fanm 
ade 20) adj, dfs nepe tee ode: te 8) al cane, Bae 
ti h (aS Aa Te ee ee mee obese trad 
= SETS beste a 
eevee oer vGpnody Baia. aie towne va a 
wig els Adie 2nd Setepde ree ea nomend Le 
Get. QONT= 82 A ee staat Pr astatots 
oven Stinele to eee eens: | ageekaagee .a (2 <vobex 
One etoene eflanss smal awtapeeesvaseds' aotteibexy 
éSlevigosgd ta: ’Agn Gp. gee ~esieee Tees pn idetheneag 
¥ \ RSS Ts RRSe Bee sai, a) 
00\.@at i400 [96 Foca eer muuctisoud. <b. .iteisaee 
8 gh 908 oo lod cK GmebwSEE estaba make ative, py Ved 
.200t B@0t-BE0t 
wenad? O°") 4? aBieine bas. .W onlin «A mangiae 
asibeades adf .wiidchy OTS) ebatad) gel qan oalage wikiee 


ie ie 


Underwriter, June 15th, 1970. 

Whitham, K., The estimation of seismic risk in Canada. 
Geoscience Canada, 2, 133-140, 1975. 

Whitham, K., Protection of the public from earthquake 
hazards in Canada. Geophysics. 12, 1-5, 1972. 

Walcott, R. I., The gravity field of Northern Saskatchewan 
and Northeastern Alberta. Department of Energy, Mines 
and Resources Observatories Branch, Canada 1-21, 1968. 

Wyss, M., Towards a physical understanding of the frequency 
distribution. Geophys. J. R. Astro. Soc. 31, 341-359, 
1973", 

Wyss, M., and Brune, J. N. Seismic moment, stress, and 
source dimensions for earthquakes in the 
California-Nevada region. J. Geophys. Res. 73, 
4681-4694, 1968. 

Wyss, M. Stress estimates for South American shallow and 
deep earthquakes J. Geophys. Res. 75, 1529-1544, 1970. 

Wetmiller, R. J. Microseismicity in the Rocky Mountain House 
seismogenic zone, Western Canada. Abstract G.A.C./C.G.U. 
Meeting May 1981, Calgary Energy, Mines and Resources, 


Ottawa. 


HAG OTD Tipe skin end Te ele 
SKOD BP, = cho (eyedqoad: yi ry | 

asd sadeen nvodaton I biets yiivetg a 
seis Oem aie dua te 2 ne engin ested hal ’ 
sa0 'S>? ebete. Heed aeltoseee Ha 


7 


‘sreupet edd de gr Laas ect ado fwaet 
‘PRE=POE 18 OR) LeaRRAE W: a | 
ne eaohd a’ Somme alita ~ i it es bau, i 
wis ik motenpid ice tot een 
OX) ten aye > b aa vens: abavsii~d LGD 
nt G81, - 18 4 

meee ae teen dsuioe sods wozerizes apne 
e8r ee -O98) RN wae sandgoas Boneisundozas 
19H “tetaeow Wesel, meet he iio imdiaagsouM ter re 
aD 2 Aw Jee Ta Lauieticn. sedan 1208, ohne aie 


yun ‘gmt “ee rr ee oe riage eae’ ean paiea | ay 
a 


a 


Me > ia ae 
¥) ; j 
. 
‘op 
at i 4 ‘ z= 
i 


7. Appendix 1 
Plots of the first 52 hours of aftershocks from the Oaxaca 


Earthquake. 
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8. Appendix 2: Program listings 
LEE program calculates the spectra of the aftershocks 
of the Oaxaca earthquake of November 29, 1978, recorded in 
the RESMAC system. The data is stored on tapes 001400, 


001401, 001402, and 001403. 


LEE2 program is an iterative program to read and plot 


the data recorded on the RESMAC system. 


KADATA.S program calculate the spectra of S or P body 
waves of the events from the Rocky Mountain House earthquake 


Swarm recorded in the DR-100 Sprengnether digital system. 
RDANIEL program calculate the spectra of the body waves 


of the earthquakes detected at the digital permanent station 


of Edmonton, from the Rocky Mountain House earthquake Swarm. 
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Steereeeeeesaeexresess« LEE PROGRAM Fe e882 Ee Ete EKER 


ALL 


YOU NEED TO MOUNT THE TAPE LIKE THIS 


S$MOUNT 001401 STP #T*® SIZE=S12 MODE®800 


EVENTS 019,020,021,035 ARE IN THE FILES 26,36,62,217 


OF THE TAPE 001401 


THE DATA 1S IN THE TAPES 001400 001401 001402 AND 001403 


OUTPUT FILES 


7 = -SAM SAMPLE LENGTH 
8 z -=SPEC LOG OF THE SPECTRA CORRECTED FOR INSTRUMENT 
DISTANCE AND ATENUATION (0). 
1 = -Q SEISMOGRAM 
= -SPEC SPECTRA CORRECTE BY ATENUATION AND DISTANCE 
BUT WITHOUT THE CORRECTION OF THE INTRUMENT. 
t2. = -SsWw SPECTRA WITH CORRECTION OF INSTRUMENT DISTANCE 


AND ATTENUATION (Q). 


INTEGER ILFTL,ILFTR,ILSTL,ILSTR, IARR(SO) 
DIMENSION D(8) 

INTEGER*2 18(512),1B82(512),1E(40000) 

REAL ISAM( 40000) ,AI (4000) ,AR(4000) 

INTEGER SHFTR,SHFTL,ALL,ALL1,ALL2,VYES/’YES ‘/ 
INTEGER ILF(100),1LS(100),IL,IR, IMA, IMAN, IEXP,LNUM 
INTEGER#2 LX,1S(512) 

REAL AMP(3000),FF(6000) ,DAT( 7000) ,DATC(6000) 
REAL A( 10000) ,FFLOG(6000) ,SAM(10000) 

DOUBLE PRECISION AMPLOG(6000}) ,AMPOUT(6000) 
COMPLEX C(3000) 

INTEGER NN,MINI,NBS,NB 

LOGICAL#®1 FREE(1)/‘2'/ 

DATA YES/‘Y’/,ANO/‘N‘/ 


THIS IS THE CODE REQUIRED TO READ THE RESMAC TAPES 


REAL PI/3.141592654/ 

FORMAT (A1) 

WRITE(6,10) 

FORMAT(’ WHICH FILE DO YOU WANT 7°) 
READ(5,FREE,END=64) IWF 

CALL SKIP(IWF,0,8) 


READ FIRST BLOCK, (DECODE THE EVENT INFORMATION) 


CALL READ(IB,LXK,0,LNUM, 8) 
00 1 T=1,45 
ALL1=1B(1) 
ILFTL=SHFTL(ALL1,16) 
ILF(2*1-1)=SHFTR(ILFTL,24) 
ILFTRESHFTL(ALL1, 24) 
ILF(2*1)=SHFTR(ILFTR, 24) 
CONTINUE 
WRITE(6, 300) 

AT(’ YEAR MON DAY HOUR “SECT OKT) 
WRITE(6,200) (ILF(I),1=2,7) 
WRITE(6,400) 
FORMAT(‘ NUM OF STA THAT DETECTED THE EVENT’) 
WRITE(6,499) ILF(39) 
FORMAT(15X,14) 
WRITE(6,500) 
FORMAT(’ STATIONS THAT DETECTED THE EVENT *) 
po 2 1=40,45 ‘ 
PLFCL)SILF C1) /2 
WRITE(6,200) (ILF(I),1=40,45) 
FORMAT(1X,814) 


READ SECOND BLOCK 


CALL READ(I182,LXKX,0,LNUM,8&) 
00 3 F=1,7 

ALL2=1B2(1) 
ILSTL=SHFTL(ALL2,16) 
ILS(2*#1-1)SSHFTR(ILSTL,24) 
ILSTR=SHFTL(ALL2, 24) 

ILS (2*1)SSHFTR(ILSTR,24) 
CONTINUE 

WRITE(6,550) 


DON’T TAKE IN CONSIDERATION THE SSO FORMAT BECAUSE THIS WAS 


sso 


CHANGED 


FORMAT(’ DIS OF STA IN THE SECOND BLOCK ‘) 
WRITE(G, 200) (ILStT) (123, 13), 2) 
WRITE(6,600) 
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THIS IS IMPORTANT BECAUSE IT TELL YOU HOW MANY BLOCKS IN EACH 
STATION WERE RECORDED. FOR EXAMPLE IF 6 STATIONS RECORDED AN 


EVENT YOU CAN GET 


600 


702 


700 


701 


804 


801 


13 


12 


46 


11 


1m THIS PROGRAM WH ARE ASSUM 
t Gaam ASK FOR YHE 


ARRAY IN BLOCKS OF THE STATIONS 
3 48 33 0 138 183 228 


THIS MEAN 


CHANEL BLOCKS STATION 
° 3 : 48 CRX 
1 49 . $3 Acx 
2 94 - 138 MEX 
4 138 - 183 1G Ue g 
5 186 - 228 ric 
6 229 - ENO OF FILE ITIP 


THEREFORE IF YOU WANT STAIGN 4 (IIT) THEN YOU NEED 
TO SKIP 138 BLOCKS ANO THEN YOU ARE AT THE BIGINING 
OF IIT 


FORMAT(’ ARRAY IN BLOCKS OF THE STATIONS ’) 
BG 4 12834,%80 

ALL2e1B2(1) 

ILSTLSSHPT. (ALL2, 18) 

ILS( 291-1) SSWMFITR(ILSTL, 24) 

5 ta.2 Semmninineemmemneteenaeeeeee aetna 
ILSTRESHFTR(ILSTR, 28) 
ILS(2Zel)eSwPTLIILSTR, 8) 

TARR IT) SiLS(2st-1)o1es( 2st) 
WRITE(G,98) (TARA( I), 1834, 40) 
FORMAT(1X,718) 

WRITE(S, 702) 

PORMAT(’ DO YOU WANT ANOTHER FILE?’ ) 
REAO(S,1000) DE 

IF(08.890.Y8S) GG TO 9 


AT THIS POINT, NOW, YOU KNOW IN WHICH BLOCKS ARE THE 
THE STATIONS. THEREFORE, YOU CAN CHOGSE THE STATION 
THAT YOU WANT 


WRITE(6, 700) 

FORMAT(’ HOW MANY BLOCKS O06 YOU WANT TO SKIP 7? ‘) 
REAO(S, FREE, END=64) NBS 

CALL SKIP(0O,NBS,8) 

WRITE(6, 701) : 

FORMAT(’ HOW MANY BLOCKS DO YOU WANT TO READ 7°) 
REAO(S, FREE, ENOBG4) NB 

NB2SEENB 

CALL GETSAM(N,NB,ISAM,A) 

WRITE(6,804) N 

FORMAT(’YOU HAV®’,16,’ SAMPLES !°) 

WRITE(6, 800) 

FORMAT (‘HOW MANY SAMPLES DO YOU WANT TO PLOT 7°) 
REAO(S, FREE, ENOSE64) NN 

WRITE(10,12) NN 

WRITE(6, 808) 

PORMAT( ‘FROM WHICH SAMPLE 00 YOU WANT TO PLOT?’) 
REAO(S, FREE, ENOSE64) ISAIN 

WNMNSISAINONNM- 1 

WRITE(10,101) (A(I),1SAM(I), 18I1SAIN, NNM) 
WRITE(S, 803) 

PFORMAT( ‘DG YOU WANT TO GET A SPECTRUM ?.°’) 
RBAO(S,1000) DES 

IF(08S.290.AN0) GO TS 64 

WRITE(6 , 801) 

FORMAT (‘SAMPLING INTERVAL MINIT & MFIN ‘) 

REAO(S, FREE, ENOBE4) MINI ,MFIN 

NNGMFETNoMINI 

060 13 I81,NN 

JEMINI-1¢1 

DAT(I)SISAM( J) 


TME SIGNAL IS FILTERED 


CALL SNOPAS(0.2,7.9,27.77,0,G) 
CALL FILTER(OAT,NN,O,G,1) 
WRITE(7,12) NW 

FORMAT (14) 

oo 68 181,N 

Gal 

WRITE(7,107) B,DAT(I) 

PoRmmaTt (2612.6) 

conTinug 


3-? Time, SHEAR WAV@® VELOCITY 
VS (CQ P WAVE VELGCITY VP) ABB ATTERVATICON O. 


WRITE(S, 802) 

FoRmaT(’ S-®, vS OR VP 8B @ *) 
READO(S,FREE,ENOS64) SP,V0,9 
HOIS=8.882SP 

OTIs1./38 

ma2048 

LHeM/2¢1 
DaLFst./(FLOAT(M)SOTI) 

oo 46 I23,LH 

Jute 

FR(J)e( 1-1) sOGLF 

CALL CONDI (OAT,NNW,DATC,M) 
60 11 Is1,M 
C(1)sCmPLZ(DATC(I),0.0) 
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Bsaeseezesesezzs PROGRAM LEEZN Sa2eeszesesesesessse 


THIS PROGRAM READ ANDO PLOT THE RESMAC DATA 


nannan 


INTEGER ILFTL,ILPYR,ILSTL,ILSTR, IARR(5O) 
OIMENSION 0(8) 

INTEGERS2 13(512),182(812) 

INTEGER #4 I1$AM( 40000) 

INTEGER SHPTR,SHFTL,ALL,ALL1,ALL2 


INTEGER ILF(100),I1LS(100),I1L,IR, IMA, IMAN, IEX®, LNUM 


INTEGER#2 LX,19(812) 
INTEGERS4 A( 10000) 
LOGICALs1 FPREE(1)/‘2°/ 
DATA YES/’°Y*/,AMNG/°N‘/ 
REAL PI/3.141592884/ 
CALL PLOTS 

68 TO 1000 

CALL SKIP(0,°2,8) 


HERG THE COMPUTER TELL YOU IN WHICH FILE YOU ARE 


nnane 


CALL FINCMO(’SOISPLAY *T#° ,12) 
wRITe(s, 12) 

12 FORMAT(’ DO YOU WANT TO QUIT 7? Y/N’) 
R@a0(8,13) 08 
IF (D0€.890.Y@8) GO TO 8&3 


c 

i+} SKIP OWB FILE APTER THE FILE IN WHICH THE 
c "POINTER" IS POSITIONED; OTHERWISE STOP 
c 
1 


e060 WRITE(6,10) 
10 PORMAT(°OGO YOU WANT THE NEXT FILE ? Y/N’) 
RBA0(S,13) OF 
IF(DE.£9.AaN0) GO TO 63 
Tweet 
CALL SKIP(IWF,0,8) 


READ FIRST BLOCK 


anna 


CALL R@AO(IB,LX,0, LNUM, 8) 
po 1 181,46 
ALLISIB(I) 
ILFTLSSHFTL(ALL1,16) 
ILF (281-1) SSHPTR(ILFTL, 24) 
ILFTR@ESHFTL(ALL1, 24) 
1 ILFP (281) SSHFTR(ILPTR, 24) 
WRITE(6, 300) 
300 FORMAT(’ YEAR MONTH OCAY HOUR MIN SEC ’) 
WRITE(6,200) (ILF(I),I#2,7) 
WRITE(6, 400) 
400 FORMAT(‘ NUM OF STA THAT DETECTED THE EVENT’) 
WRITE(6, 488) ILF(39) 
488 FoRmmaT(16x,14) 
wRITE(S, 500) 
see FormaT(’ STATIONS THAT OSTECTED THE EVENT ',/, 
e° OeCRK 18ACX WeMEX 4BIIT SelIC SSIIP’) 
06 2 1840, 48 
= 


—_ —-- --=- -- a ~~. 


wRITH(G,200) (ILF(1), 1840, 48) 
2@@ FORMAT(1X, 814) 


R@AD SECOND BLBCK 


nan 


CALL READ(182,LxX,0, LNUM, 8) 
oo 3 [81,7 

ALL2e1B2(1) 
ILSTL@ESMPTL(ALL2Z, 16) 

ILS( 281-1) SSHPTR(ILSTL, 24) 
ILSTRESHFTL(ALL2, 264) 
ILS(2*1) SSHFTR(ILSTR, 24) 
WRITE(6, 800) 


Event you CAN GET 


ARRAY IN BLOCKS OF THE STATIONS 
3 #48 33 #0 138 183 228 


THIS MEAN 


CHANEL BLOCKS STATION 
° 3 = 48 CRX 
1 4s ° #3 acx 
2 94 O 138 MEX 
4 139 - 183 IIT 
Ss 184 - 228 tte 
6 223 - EMG OF FILE Ite 


TREREPERE IF YOU wast STAIOM & (IIT) THEN 
YOU NEED TO SKIP 138 BLOCKS Ano 
Tram YOU ARG AT THE BIGINING 
or 11T 


nnannnananannannanannanannnn w 
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STATION WERE RECORDED. FOR EXAMPLE IF 6 STATIONS RECORDED AN 
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esa 
172 
173 
174 
175 
176 
V77 
178 
179 
189 
181 
182 
183 
184 
185 
186 
187 
188 
1839 
190 
191 
192 
ENO OF FILE 


600 


a 


39 


oaonnda 


801 


604 


605 


3000 
2001 
2000 


nana 


oan 


V1 


13 
63 
64 


HERE 


HERE THE TAPE 
BLOCKS YOU WANT 


FORMAT(’ ARRAY IN BLOCKS OF THE STATIONS 
00 4 1234,40 

ALL2=2=I1B2(1) 
ILSTLSSHFTL(ALL2,16) 
ILS(2*1-1)SSHFTR(ILSTL,24) 
ILSTRESHFTL(ALL2, 28) 
ILSTRESHFTR(ILSTR, 24) 

ILS( 221) 2SHFTL(ILSTR,8) 
IARR(I)2ILS(221-1)41LS( 221) 
WRITE(6,39?)> (IARR(1),1234,40) 
FORMAT(1X,718) 
CRXSIARR(35)-1ARR(34)+1 
CRXS=IARR( 34) 

CRXEZIARR(3S) 


IT TELL YOU HOW MANY BLOCKS HAS EACH STATION, 


BLOCK START ANDO IN WICH BLOCK END 


WRITE(6,601) CRX,CRXS,CRKXE 


) 


FORMAT(‘# OF BLOCKS IN CRXE‘,F8.3,’ STARTS’, F8.3,° 


ACX=1ARR(36)-IARR(35)+1 

_ ACXSZIARR(35)41 
ACXEZIARR(36) 
WRITE(6,602)ACX,ACXS, ACXE 


FORMAT(‘# OF BLOCKS IN ACX=’ ,F8.3,' STARTE’ ,FB.3,’ 


RMEXZ1ARR(38)-1ARR(36)%1 
RMEXSZIARR(36)+1 

RMEXESIARR (38) 
WRITE(6,S03)RMEX, RMEXS , RMEXE 


FORMAT(’# OF BLOCKS IN MEX#’,F8.3,° START#’ ,F8. 


RIITZIARR(39)-TARR(38) +1 
RIITSSIARR(38)+1 
RIITEsIARR(39) 

WRITE(6 ,GOS4)RIIT,RIITS,RIITE 


FORMAT(‘w# OF BLOCKS IN I1T2’,F8.3,'‘ START=’,F8. 


RIICZIARR(40)-IARR(39)+1 
RIICS=IARR(39)+1 
RIICE®IARR(40) 
WRITE(6,6OS)RIIC,RIICS,RIICE 


FORMAT(’# OF BLOCKS IN IIC#’,F8&.3,‘ START2’,F8. 


RIIPzS 

RIIPS&IARR(40)41 

RIIPE2S 
WRITE(6,6O6)RIIP,RIIPS,RIIPE 


FORMAT(’# OF BLOCKS IN IIP#®’,FE.2,° STARTS’ ,FE. 


ASK FOR THE NUMBER OF BLOCKS TO BE READ 


WRITE(6,501) 


FORMAT(/,’ HOW MANY BLOCKS DO YOU WANT TO READ 


REAO(S,FREE,ENO#64) NBLCKS 


3,’ 


3,’ 


cps 


2a) 


IN WHICH 


ENO#’,F8&.3) 


END#’ ,F&.3) 


ENOs’ ,F8&.3) 


ENDs’,F8&.3) 


ENOs’ ,F8.3) 


END2’,F6.1) 


CALCULATE THE NUMBER OF SAMPLES, EACH BLOCK HAS 256 SAMPLES 
THAT IS EQUIVALENT TO 8&8 SECONDS (32 SAMPLES PER SECOND) 


NSECZ2NBLCKS#8 
NS 2S62NBLCKS 
WRITE(6,3000). NBLCKS,NSEC 


FORMAT(14,’ BLOCKS IS EQUIVALENT TO ‘,14,’ 


WRITE(6, 2000) 


FORMAT(’HOW MANY BLOCKS 00 YOU WANT TO SKIP ?°) 


REAO(S, FREE, ENDZ64)I1STA 

CALL SKIP(O,ISTA,8) 

CALL GETSAM(N,NBLCKS,ISAM,A) 
TREWe ISTA4NBLCKS 


CALL SKIP(0,-IREW, 8) 
CALL FTNCMDO( ‘SDISPLAY #Tz’ ,12) 


THE OATA IS STORED IN "ISAM" 


CALL PLTTRC(ISAM,N) 
WRITE(6,11) 


IS REWINDED SO THAT YOU CAN SKIP WHATEVER 
IN OROER TO READ OTHER STATIONS 


FORMAT(‘DO YOU WANT TO PLOT OTHER STATION ? Y/N’) 


READ(5,13)08S 

IF (DBS.EO0.ANQ) GO TO 9 
WRITE(6,500) 

WRITE(G,200) (ILF(1I),1240, 48) 
WRITE(6,600) 

WRITE(6, 30) (TARATT), 1234, 40) 
WRITE(6,6801) CRX,CRXS,CRKE 
WRITE(6,602)ACX,ACXS, ACXE 
WRITE(8, 603) RMEX ,RMEXS , RMBEXE 
WRITE(6,SO4)RIIT,RIITS,RIITE 
WRITE(6,GOS)RIIC,RIICS,RIICE 
WRITE(6,SO06)RIIP,RIIPS,RIIPE 
GO TO 2001 

FORMAT(A1) 

CALL PLOT(O.,0.,999) 

sToPp 

END 
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eae eteresaxe2eeez PROGRAM KADATA S 


WRITTEN BY CECILIO y. 


REBOLLAR ON SUMMER 1979 


SSCs Eraeerewresezrzrx 


THIS PROGRAM CALCULATE THE BODY WAVE SPECTRA OF S WAVES 
OF THE EVENTS RECORDED BY £. KANASEWHICH. THESE EVENTS 


ARE STORED 
READ THE NECESARY 


IN THE TAPE GO8O025S VOL=EROCKY! 
INFORMATION FROM THE FILE EV.S.i ¢€. 


THE PROGRAM 


NBLKS NUMBER OF BLOCKS T0 BE READ 
NN MUMBER OF SAMPLES 
io) ATENUATION 
STRT STARTING POINT OF THE SAMPLING INTERVAL 
sp S$-P TIME 
vs S-WAVE VELOCITY 
OSNAME DATA SET NAME 
OUTPUT FILES 
“2 SPECTRA OF VERTICAL COMPONENT 


-N SPECTRA NORTH-SOUTH COMPONENT OR TRANSVERSE 
COMPONENT OR RADIAL 
-Z7St SAMPLE INTERVAL OF 2 N-S 8 E-W COMPONENTS 


=e SPECTRA EAST-WEST 


COMMON /PLTLOG/WIOTH,SCP,SCT 


REAL *®4 2( 1020) ,NS( 1020) ,EW( 1020) ,7(1020) ,26(1024) ,SNC( 1024) 


»WEC( 1024) ,FF(S12),A82(512),AS(512) ,AW(S12) 


INTEGER POSN(8)/’SCON‘,’ 


LOGICAL#1 OSNAME(13) 


EQUIVALENCE(POSN(S) ,OSNAME(1)) 


ST che 


P’,’QSNE’, 42° 


COMPLEX C2C(513),CSNC(S13) , CwEC(S13) 


DIMENSION 0(8) 
INTEGER 9,STRT, IWK(9) 
OATA ANO/’N’/ 
REAL PI/3.1418$926S84/ 


au 


REAL#& 271(1024) ,NS1(1024) ,EW1(1024) ,MINNS1,MAXNS1,MINEW! 


»MAXEw! 

INTEGER UNIT/4/,BLOANK/’ A Uf 
INTEGER*2 LEN 

LOGICAL#1 FREE(1)/’‘2’/,LOGX,LOGY 


NUMBER OF 


CULATE 


Cc 
S HERE WE START 
c 
OTI=1/100. 
WRITE(6,1) 
1 FORMAT( ‘ENTER THE NUMBER OF EVENTS TO BE SKIPPED’, 
* NUMBER TO BE RUN’) 
REAO(S,FREE,ENO=64) NSKIP,NEVNTS 
IF(NSKIP.EO.0) GO TO 100 
00 2 121,NSKIP 
2 READ(12,3) NBLKS,NN,O,STRT,SP,VS,DSNAME 
3 FORMAT(12,314,2F5.2,1381) 
100 00 4@ IIT=1,NEVNTS 
c 
C READ IN THE NECESSARY INFORMATION: NUMBER OF 8LOCKS, 
c SAMPLES, STARTING POINT, S-P TIME, S VELOCITY 
¢ AND DATA SET NAME 
(3 
REAOTIZ, 3) NBLKS,_NW,O,STRT,SP,VS, OSNAME 
e 2 
C SET THE TAPE TO THE APROPIATE DATA SET NAME 
Cc 
CALL FTNCMO(POSN,29) 
CALL GETDAT(Z1,NS1,EW1,N,NBLKS)} 
c 
C FINO MAXIMUN OF NS ANO EW COMPONENTS IN ORDER TO CAL 
c LOCAL MAGNITUDE 
S 
WRITE(6, 400) 
400 FORMOT(’DO YGU WANT MAX AMPLITUDE OF NS & EW COMPONENT’) 
READ(S,201) DE , 
IF (DE.EO.ANO) GO TO 402 
CALL USMNMX(NS1,N,1,MINNS1,MAXNS1) 
CALL USMNMXK(EW1,N,1,MINEW! ,MAXEW1) 
WRITE(6,401) MAXNS1,MAXEW! 
401 FORMOT( ‘NS MAX=’,F12.4,' EW MAX2’,F12.4,/) 
c 
C HERE THE PROGRAM ASK IF YOU WANT TO CALCULATE THE SPECTRA 
c 
WRITE(6, 403) 
403 FORMAT(’DO YOU WANT TO CALCULATE THE SPECTRA 7? Y/N‘) 
READ(S,201) DE : 
IF(0E.EQ.ANO) GO TO 6&4 
402 Ju20 
NENOSSTRT+NN 
00 20 JK=STRT,NENO 
JJtduet 
Z(dsIzZ1( NK? 
NS (JJ) =NS1( JK) 
20 EWl UJ) = EWI( UK) 
00 21 KK#1,NN 
21 T(KK)=KK 
WRITE(13,101) NN 
WRITE(13,102) (T(II),2¢11),1121,NN) 
WRITE(13,101) NN 
WRITE(13,102) (T4II),NS(II1),1121,NN) 
WRITE( 13,101) NN 
WRITE(13,102) (T(I1),EWCIT),11=1,NN) 
c 
C IF YOU KNOW THE AZIMUTH YOU CAN ROTATE YOUR SIGNALS 
6 
WRITE(6,300) 
300 FORMAT( ‘DO YOU WANT TO ROTATE THE SEISMOGRAMS 7°) 
READ(5S,201) OE 
IF(0E.E0.A8N0) GO TO 304 
WRITE(6, 302) 
302 FORMAT( ‘GIVE THE AZIMUTH !|1°) 
READ(S,FREE,ENO=64) AZIM 


CALL ROTATE(EW,NS,NN,AZIM) 
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180 (s 


ENO 


hates) 


181 
182 
183 
184 
185 
186 
187 
188 
189 
190 
191 
192 
193 
194 
198 
196 
197 
198 
199 
200 
201 
202 
203 
204 
205 
206 
207 
208 
209 
OF FILE 


NOW EW = 


AZIM = AZIMUTH 
ALL THREE COMPONENTS ARE FILTERED 


CALL BNOPAS(1.0,45.0,10.0,0,6) 
CALL FILTER(2,NN,0,G,1) 

CALL FILTER(NS,NN,O,C,1! 
CALL FILTER(EW,NW,D,C,1> 


HERE WE START THE CALCULATION OF 


THE SPECTRUM ANDO WE ASSUME A VP OF 6.0 KM/SEC 


HOIS=8 22*SP 


is 
c HERE WE ARE GENERATING & VECTOR IN FREQUENCIES FF 
c 


M22S6 

MMsM/2 

LHSM/241 
OELFs1./(mMF0TI) 
Bo 10 123,LK 
J=I-1 


RADIAL COMPONENT & NS = TRANSVERSE COMPONENT 


IMLS 


HENCE 


10 FF( J)SJPOELF 
c 
c REMOVE THE TRENO, AVERAGE, ADD ZEROES TO AN EVEN 
Cc NUMBER "M" BNO SMOOTH THE SIGNAL AT THE ENDS WITH 
c A COSBELL (SEE KANASEWICH BOOK) 
c 
CALL CONGI(2,NN,ZC,M) 
CALL CONOI(SN,NN,SNC,M) 
CALL CONOI(EW,NN,WEC,M) 
c 
c SET THE DATA FOR THE USE OF THE FFT FFTRC FROM THE 
c LIBRARY 
c 
CALL FFTRC(Z2C,M,C2C, IWK, WK) 
CALL FFTRC(SNC,M,CSNC, IWK,WK) 
CALL FFTRC(WEC,M, CWEC, IWK ,WK) 
BO) Wt Letts 
AZ(I1)SREAL(CZC(1)) 
AS(1)=REAL(CSNC(1)) 
"1 Aw(1)=REAL(CWEC(1)) 
LLStHet 
po 12 IT=LkL,M 
Jel-CH! 
AZ(1)=M&1MAaG(CZC( J)? 
AS(1)2AIMAG(CSNC(J)) 
iy 3 AW(1)2AIMAG(CWEC(J)) 
LOG2N28 
CALL AMPPHZ(LOG2N,482,1,P0CZ,PFNZ) 
CALL AMPPHZ(LOG2N,4S,1,P0CS,PFNS) 
CALL AMPPHZ'LOGZN,4W,1,P OCW, PFNW) 
c 
c You CAN SMOOTH YOUR SPECTRUM IF YOu WANT 
l= (SEE KANASEWICH 300K) 
WRITE(6, 200) 
200 FORMAT( ‘00 YOU WANT TO SMOCTH THE SPECTRA’) 
READ(5,201) OES 
201 FORMAT( 41) 
IF(O0ES.EQ.ANO) GO TO 202 
wRITE(S,14) 
14 FORMAT( ‘WE NEED THE WINOOW .N FREQUENCY DOMAIN FOR OANIEL '') 
READ(S,FREE,END=64) FWINO 
CALL DANIEL (AZ,M,FWINOD, 12.0) 
CALL DANIEL(AS,M,FWINO, 190.0) 
CALL DANIEL (AW,M,FWINOD,10.0) 
cS 
THE SPRENGNETHER INTRUMENT WAS WORKING WITH A GAIN FROM 60 TC 
120 0 N 6 9068 Ss Pisic HE GAIN WAS SET UP 0 66 OB. 
C we NEED TO MULTIPLY THE OUTPUT BY 2**9 IN ORDER TO RECOVER 
c THE SIGNAL 
c 
202 oO fs 225MM 
V=AMAGN(FF(I)) 
ATENUAZEXP( (PIFHOISZFF(I))/(O8VS)) 
AZ(1)2¢(AZ(1)2ATENUA)/V)=(2%29) 
AS(1)=((AS(1)*ATENUA)/V)2(2**9) 
AWw(1)=( (AWC TL) 2 ATENUA)/V)2 (2829) 
AZ(1)2ALOG1IO(AZ(1)) 
AS(1)2O8L0G10(AS(1)) 
Aw(l)sAlLoGcio(aw(!)) 
FF(I1)s@LOGIO(FF(I)) 
13 CONTINUE 
MMM= 8S 
WRITE(9,101) MMM 
WRITE(9, 102) (FF( IT), AZ(1),126,90) 
WRITE(10,101) MMM 
WRITE( 10,102) (FF(1),AS(1),1=6,90) 
WRITE( 11,101) MMM 
wRITE( 11,102) (FECT), AW(1), 126,90) 
101 FORMAT(14) 
102 FORMAT(2G12.6) 
wRITE(6,5) OSNAME 
5 FORMOT(’I FINISHED WITH On ec tI 
4 CONTINUE 
64 sTtoPp 
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1 SR =FORTG T=6 
2 c 
3 C sBeeezeseeaeeeee PROGRAM RDANIEL PES ESE ELS SES FS 
4 c 
s C THIS PROGRAM CALCULATE THE SPECTRA OF EARTHQUAKES FROM THE 
6 C ROCKY MOUNTAIN HOUSE EARTHQUAKE SWARM RECORDED AT EDMONTON 
7 C USING A DANIEL-SPECTRAL ESTIMATE FROM 
8 C THE PERIOOROGRAM, USING A SUBROUTINE BY 0. GANLEY. 
9 C THE SEISMOGRAMS ARE IN THE TAPE GO8067 VOLSMANOE7 
10 C AND THE PROGRAM READ ALL THE INFORMATION FROM 
11 (2 THE FILE “"ROCKYF” 
12 c 
13 DIMENSION Z(6000) ,SN( 6000) ,WE(6000) , VB(6000) 
14 REALs4 FF(1024),2C(1024) ,SNC( 1024) ,WEC( 1028) ,AZ(1024) ,SVC( 1024) 
15 1,812(1024) ,AS( 1024) ,ASI (1024) ,AW( 1024) ,AWI (1024) ,$V¥( 1024) 
16 2,AMPZ(1024) , AMPE( 1024) , AMPW( 1024) ,T(1024),AZIM,WR(9) 
17 COMPLEX CZC(102S8) ,CSNC( 1025) ,CWEC(1025) ,cSVC(1025) 
18 REAL®4 Av(512) 
19 INTEGER QO,NN, IWK(9) 
20 OIMENSION 0(8) 
21 REAL P1I/3.141892684/ 
22 INTEGER RC/4/, 
23 . CONREG (27) ,OSNAME(S),BLANK/' Cf Vee. LPS CLP; 2CUIP / CLIP 
24 »POSN(7)/°POSN‘ ,68° uy) 
2s INTEGER MOUNTC(6) 
26 OATA YES/‘Y’/,ANO/‘N’/ 
27 INTEGER TAPE/’TAPE’/ 
28 LOGICAL#1 FREE(1)/’2’/ 
23 c 
30 EQUIVALENCE (OSNAME(1), POSN(3)) 
31 EQUIVALENCE (MOUNTC(1), CLIP), (MOUNTC(2), DSNAME) 
32 c 
a3 INTEGER®2 CONLEN/28/, MNTLEN/24/ 
34 c 
35 c READ IN THE NECESSARY INFORMATION: 
36 c LONG OR SHORT PERIOO DATA 7, NUMBER OF SAMPLES, ATENUATION 
37 3 STARTING POINT IN SEC, S-P, S-WAVE VELOCITY 
38 c &@ OATA SET NAME 
39 c 
40 1015 FORMAT(’ ENTER THE NUMBER OF EVENTS TO BE SKIPPED, NUMBER’, 
41 .' TO BE RUN‘) 
42 WRITE(6,1015) 
43 READ(S,FREE)NSKIP,NEVNTS 
44 DO 1010 1#1,NSKIP 
4s 1010 READ(12,97,ENO#64) LPD,NN,9,STRT,SP,VS,OSNAME 
46 DO 1011 JJJJ2z1,NEVYNTS 
47 READ(12,97,END=54) LPO,NN,O,STRT,SP,VS,DSNAME 
48 97 FORMAT(A4,214,2F5.0,F5.2,6A4) 
49 Cc 
50 C IS THE DATA LONG PERIOD DATA 7? IF NOT, THERE ARE 18 POINTS PER 
$1 c SEC. OTHERWISE, THERE ARE 3 PER SEC 
$2 (s . 
$3 DTIz1./18. 
54 IF(LPO.EO.LP) OTI1#1./3. 
55 IST#2STRT/OTI 
56 c 
57 c SET THE TAPE TO THE APPROPIATE OATA SET NAME. 
S58 (3 
‘3 CALL CNTRL(POSN, CONLEN, 4,CONREG, 8990,8995 ,8999) 
ett ERD TW YHE_GHTE_ACORSTNET3 WHETER WE WANT LONG a 
PERIOD OR SHORT PERIOD 
4 2 
64 IF(LPD.NE.LP) CALL GETDAT(Z,SN,WE,VB8,NN,4,IST) 
65 IF(LPD.£EO.LP) CALL GTDTLP(2,SN,WE,VB,NN,4,1ST) 
66 p00 9 I=1,NN 
67 9 T(1)#! 
68 WRITE(8,101) NN 
69 WRITE(8,102) (T(I1),2(1),121,NN) 
70 c : 
71 c WE ORE FILTERING BELOW F120.5 HZ ANO ABOVE 
72 Cc F2=27.0 HZ 
73 Cc 
74 WRITE(6,200) 
75 200 FORMAT( ‘DO YOU WANT TO ROTATE THE SEISMOGRAMS ?°) 
76 READ(5,201) DE 
77 201 FORMAT(A1) 
78 IF(DE.£0.AN0) GO TO 204 
79 WRITE(6,202) 
80 202 FORMAT( ‘GIVE THE AZIMUTH !!°) 
81 READ(S,FREE,ENO=64) AZIM 
82 CALL ROTATE(WE,SN,NN,AZIM) 
83 ¢ 
84 C NOW WE = RADIAL COMPONENT & SN = TRANSVERSE COMPONENT 
a5 Cc AZIM = AZIMUTH, SV(KK) IS PURE SV WAVE 
86 c 
87 po 1S KK=1,NN 
88 15 SV(KK) #SORT(2Z( KK) #*2¢WE (KK) #22) 
89 204 CALL BNOPAS(O.5,9.0,55.55,0,G) 
90 CALL FILTER(Z,NN,0,G,1) 
91 CALL FILTER(SN,NN,D,G,1) 
92 CALL FILTER(WE,NN,D,G,1) 
93 c CALL FILTER(SV,NN,0,G,1) 
94 (3 
95 c HERE WE START THE CALCULATION OF 
96 c THE SPECTRUM AND WE ASSUME A VP OF 6.2 KM/SEC 
97 c 
98 c HDIS28.37*SP 
99 HOIS=325. 
100 c 
101 c HERE WE ARE GENERATING A VECTOR IN FREQUENCIES FF 
102 c 
103 Me£1024 
104 MM=eM/2 
105 LHEM/2+1 


OELFe1./(M*OTI) 
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107 00 10 123,LH 
108 Jet-) 

109 10 FF(J)S(1-+1)*08LF 

110 (= 

111 c REMOVE THE TRENO, AVERAGE, ADD ZEROES TO AN EVEN 
112 c NUMBER “M*“ AND SMOOTH THE SIGNAL AT THE ENOS WITH 
113 c A COSBELL 

114 c 

118 CALL CONDI(2,NN,2C,M) 

116 CALL CONDI(SN,NN,SNC,M) 

117 CALL CONOI(WE,NN,WEC,M) 

118 c CALL CONDI(SV,NN,SVC,M) 

119 ¢c 

120 c SET YHE OATA FOR THE USE OF THE FFT FFTRC FROM THE IMLS 
121 ie CITBRARY 

122 c 

123 CALL FFTRC(2ZC,M,C2ZC, IWK, WK) 

124 CALL FFTRC(SNC,M,CSNC, IWK,WK) 

128- CALL FFTRC(WEC,M, CWEC, IWK, WK) 

126 ¢ CALL FFTRC(SVC,M,CSVC, IWK, WK) 

127 DO 11 I#1,LH 

128 AZ(1)sREAL(C2CE(1)) 

129 AS(1)#RBAL(CSNC(I)) 

130 11 aW( 1) RBAL (CWEC(I)) 

131 Ge AV(T)sRBAL(CSVC(I)) 

132 LL3LHe1 

133 DO 12 IsLL,M 

134 JSIeLHet 

138 AZ(1)SAIMAG(CZC(U)) 

136 AS(1)SA4IMAG(CSNC(uU)) 

137 12 AW(1)sM&IMAG(CWEC( UJ) ) 

138 c AV( 1) SAIMAG(CSVC( uJ) ) 

139 LOG2N210 

140 CALL AMPPHZ(LOG2N,A2,1,P0CZ, PENZ) 
141 CALL AMPPHZ(LOG2N,AS,1,P0CS,PFNS) 
142 CALL AMPPHZ(LOG2N, AW, 1!,POCW, PFNW) 
143 c CALL AMPPHZ(LOG2N,AV,1,POCV,PFENV) 
144 WRITE(S,14) 

145 14 FORMAT(’WE NEED THE WINDOW IN FREQUENCY DOMAIN FOR DANIEL !°) 
146 READ(S,FREE,ENOs64) FWIND 

147 CALL DANIEL(AZ,M,FWINOD,SS.55) 

148 CALL OANIEL(AS,M,FWINOD,55.55) 

149 CALL DANIEL (AW,M,FWINO,5S.55) 

150 c CALL DANIEL(AV,M,FWIND,55.55) 

151 DO 13 I1=2,MM 

152 VS AMAGN(FF(1)) 

183 ATENUAZEXP( (PIZHOIS*FF(I))/(O2VS)) 
184 AZ(1)#(AZ(1)*ATENUA)/V 

155 AS(1)=(AS(1)*ATENUA)/V 

186 AW( I) 2(AW(1)#ATENUA)/V 

187 c AV(1)2(AV(I)#ATENUA)/V¥ 

158 AZ(1)#AL0G10(A2(1)) 

158 AS(1)24L0610(AS(1)) 

160 AW(1)2AL0G10(aW(I)) 

161 c AV(1)SALOG10(AV(I)) 

162 FF(I)SALOGIO(FF(I)) 

163 13. CONTINUE : 

164 MP2157 

165 MWILL=200 

186 WRITE(9,101) MWILL 

167 WRITE(9,102) (FF(I),A2(1),1=120,320) 
168 WRITE(10,101) MP 

169 WRITE(10,102) (FF(I),AS(1I),1219,178) 
170 WRITE( 11,101) MP 

171 1011 WRITE(11,102) (FF(I),AW(I),1219,175) 
172 c WRITE( 13,101) MP 

173 c1011 WRITE(13,102) (FF(I),avil),1219,178) 
174 101 FORMAT(14) 

175 102 FORMAT(2612.6) 

176 GOGO TO 64 

177 999 RCZRC+6 

178 995 RCZRC+S 

179 990 CONTINUE 

180 WRITE(6,991) RC, CONREG arate 
181 ae MA *Oxet RETURN CODE FROM [EE ee ae TA 
182 + ‘ INFORMATION: * ,27A4) 

183 64 STOP 

184 END 
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BezeezrezexGet leng period data from Charlie tapestrtszexrerezaerssessrse 
SUBROUTINE GTOTLP(2,SN,WE,VB,N,IU,IST) 
This picks up long period data from CHARLIE tapes 


Also picks up the signal into VB (refer to GETDAT for 
explanation of where time signal is located). 


OIMENSION Z(N), SN(N), WE(N), VB(N) 
INTEGER*2 I1£(2700), Lx 
Kz((N*IST)/3900) +41 
Mz-IST 
00 10 I31,Kk 
CALL READ(IE,LX,0O,LNR, IU) 
DO 10 Jvet,2700,3 
Meme! 
IF(M.LE.0)GO TO 10 
IF(M.GT.N)GO TO 11 
ree SPE CS) cm CELE CUI a) ee 6) aes /) 2 
VB(M) = ITABS(ITEMP) 
Z(M)SIE(¥) 
SN(M)DETE( U1) 
WE (M)SIE(u+2) 
CONTINUE 
GO TO 12 
11 Meme) 
12 CONTINUE 
101 FORMAT(’ ‘,18,/(3G618.5)) 
RETURN 
DEBUG SUBTRACE 
ENO 


SESUSESSSSTT REESE cE sR aa THAR esa eEEEEexsaReseeeseezeueees 


= z 
= Get short period data = 
z= = 


SSSESSSSSR SSAA LASTS RSLS eee e ees eseeeeeeeereseeeserceesarss: 


SUBROUTINE GETOAT(Z,SN,WE,VB,N,IU,IST,=) 


This picks up short period date from CHARLIE tapes 


Array VB contains the time signal. This time signal is found 
im bit 14 ef the data, where bit o is the leftmost pit, and 
bit 15 is tne rightmost bit of a 16-bit integer {integer*2). 


DIMENSION Z(N), SN(N), WE(N), VB(N) 
INTEGER#®2 I1E(2700), Lx 
K2((N*IST)/900)+1 
Mz-IST 
00 10 I81,K : 
CALL READ(IE,LX,0,LNR,IU,£12,812,8120) 
00 160 J#1,2700,18 
DO 10 KK=1,6 
Meme) 
IF(M.LE.0)GO TO 10 
IF(M.GT.N)GO TO 11 
IPNT = J * KK - 1 
ITEMP = ( IEC IPNT) - ((TE(IPNT)/4) © 4) ) Ui} 
ve(M) = IABS( I TEMP) 
Z(M) & IE(IPNT) 
SN(M) = IE(J*KK45) 
WE(M) = ITE(U+KK+11) 
10 CONTINUE 
GO TO 12 
11 Meme) 
12 CONTINUE 
RETURN 
120 RETURN 1 
ENO 
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“SUBROUTINE BNOPAS(F1,F2,0ELT, DG) =e 


SUBROUTINE BY DAVE GANLEY OWN MARCH S, 1977. 


THE PURPOSE OF THIS SUBROUTINE IS TO DESIGN AND APPLY A 
RECURSIVE BUTTERWORTH BANO PASS FILTER (KANOSEWICH, TIME SERIES 


ANALYSIS IN GEOPHYSICS, UNIVERSITY OF ALBERTA PRESS, 1975; SHANKS, 
JOHN £, RECURSION FILTERS FOR DIGITAL PROCESSING, GEOPHYSICS, V32, 
PP esi SS athe SIGS IN ORDER TO DESIGN THE FILTER & CALL MUST BE 
MADE TO BNOPAS ANDO THEN THE FILTER MAY BE APPLIED BY CALLS TO 
FILTER. THE FILTER WILL HAVE 8 POLES IN THE S PLANE AND ILS 
OPPLIED IN FORWARD AND REVERSE DIRECTIONS SO AS TO HAVE ZERO 

PHASE SHIFT. THE GAIN AT THE TWO FREQUENCIES SPECIFIED AS 


CUTOFF FREQUENCIES WILL 8E -608 ANDO THE ROLLOFF WILL BE ABOUT 
96 DB PER OCTAVE. A BILINEAR Z TRANSFORM IS USED IN DESIGNING 
THE FILTER TO PREVENT ALIASING PROBLEMS. 


COMPLEX P(4),S(8),21,22 
OIMENSION 0(8),X(1),xXO(3),xXO0(3), XE(3) 
DATA ISW/0/, TwoP!/6.2831883/ 


THIS SECTION CALCULATES THE FILTER ANDO MUST BE CALLED BEFORE 
FILTER IS CALLED 


Fil = LOW FREQUENCY CUTOFF (6 OB DOWN) 

F2 = HIGH FREQUENCY CUTOFF (6 OB DOWN) 

OELT = SAMPLE INTERVAL IN MILLISECONOS 

0 = WILL CONTAIN 8 Z DOMAIN COEFICIENTS GF RECURSIVE FILTE 
G = WILL CONTAIN THE GAIN OF THE FILTER 


WRITE (6,1) €1,F°2, DELT 

FORMAT (‘1 BANOPASS FILTER DESIGN FOR A BAND FROM ’,F8.3,° TO ',FB 
12,” WERTZ." ,//7° SAMPLE INTERVAL 1S ‘FS .2,° MILLISECONDS.” ) 
OTSC0ELT/1000.0 

TOT=2.0/0T 

FOTS4.0/0T 


lSwet 

P(1)=CMPLX(-.3826834, 9238795) 
P(2)SCMPLX(- 3826834,-.9238795) 
P(3)=CMPLX(-.923879S, 3826834) 
P(S2)=SCMPLX(- 9238795 ,-.3826834) 


WISTWOPI=F! 
W2=TWOP!2F2 
Wi=TOTZTAN(WI/TOT) 
W2zTOT*TAN(W2/TOT) 
HwIDs(w2-w1)/2.0 
WWEeW 1 £Ww2 

po 19 T21,4 
Z1=P(1)sHWwId 
Z72=221=21-ww 
Z2Z=CSORT(722) 
SAT ie 2a ez 2: 
$(1*#4@)221-22 

WRITE (632) S$ 

FORMAT “OSS UPLEANE POLES ARE GAT r/o ON) Ele Oi.) eal oe, eat. Oa) 
G= S/HWI1O 

G2G2G 

G=Gr7G 

DO VS)" hs 75.2 
B=-2.07REAL(S(1)) 


re LES a) PS) Sates) MES oe Das aS a 
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23 


31 


34 


35 


36 


C=REAL (721) 

AzTOT+B+C/TOT 

G=GzA 

O(T)S(C#OT-FOT)/a 

OC l+1)2(8-2.028)/A4 

G=G7G 

WRITE (6,3) 

FORMOT (’-FILTER If (1°-22272)2824 / B812B22B83F84G' ) 
WRITE (6,4) OD 

FORMAT 2 Cae/72) BO Ne ste ee E26, one) eee ee Ee Gai eis Siar hi) 
WRITE (6,5) G . 

FiGRM AT. Mia FolietceR GARIN tS.“ Est2..°6 2) 

RETURN 

ENTRY FILTER¢(xX,N,0,GC,1G) 


OATA VECTOR OF LENGTH N CONTAINING DATA TO BE FILTERED 
FILTER COEFFICIENTS CALCULATED BY BNOPAS 

FILTER GAIN 

G = 1 MEANS TO REMOVE THE FILTER GAIN SO TH4&T THE GAIN IS UNITY 


- Ox 
ee 


KE GUSW IEG. 12.660 10) 31 

WRITE (6,6) 

FORMAT (‘’18NOPAS MUST BE CALLED BEFORE FILTER’) 
CALL EXIT 


APPLY FILTER IN FORWARO DIRECTION 


XM2EXK( 1) 

XM1=2X(2) 

XMex(3) 

XC(1)=xXM2 

KMO(2)2eM1-0(1)2xKC(1) 
MO(3)SKM-XM2-0(1)2KC0(2)-0(2)2xKC(1) 
KO(1)2KC(1) 

KO(2)2xKCO(2)°+0(3)2KO(1) 
KMO(C3)SKC(3)-°+KC(1)-0(3)2KO0(2)-0(4)2K0(1) 
XE(1)2xKO(1) 

XE(2)2K0(2)-0(5)2KE(1) 
KE(3)2X0(3)*KO0(1)-O0(8)2KE(2)-0(6)2xXE(1) 
K(1)SKE( 1) 

M2 ee BC 2h= Oi 7 SKC WD 
XCBISKECZI-NE C1) SOCTIFRN(2)-°-0( 8) EK I1) 
oo 393 =4,N 

XM2=XM1 

XM1=KM 

XM=X( 1) 

Keeble (ies) 3) 23 

GO) 0) (347, 38 73605. 

Mel 

M1i=3 

M2=2 

GO TO 37 

M=2 

Misi 

M2=3 

GO TO 37 

Me3 

Mi=2 

M2=1° 
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“v2t 3? MOCIMD=XM-KMZ-OC1)DEXCIMT) -Di2) = KC (M2) 
122 XD(M)SKE(M)-KXC(M2)-0(3)2KD(M1)-0(4)2XD(M2> 
rare KE(M)=XO(M)-KO(M2)-0(5)2KE(M1)°0(6)®xXE( M2) 
124 39 KCL) EKE(M)-KE(M2)-0(7)8K(1-°5)-0(8) EX 1-2) 
125 c 
126 Cc FILTER IN REVERSE OIRECTION 
127 c 
128 XM2=X(N) 
129 XM1SxX(N-1) 
130 XMEX(N-2) 
ETS! XC(1)=2xXM2 
132 XC(2)2xXM1-0(1)2KE(1) 
133 XC (3) 2KXM-KM2-°0(1)2xKO(2)-0(2)2xC(1) 
134 XO(C1)SxC11) 
135 XO0(2)2xKC'2)-0(3)2xD(1) 
136 KO(3)2KXO(3)-KC(1)-0(3)8xX0(2)-014)2xX0(1) 
137 XE( 1) =SKO0(1) 
138 XE(2)2xX0(2)-°0(5)2xXE(1) 
139 XKE(3)2XO0(3)°xX0(1)-0(5)2XE(2)-0(6)2xXE(1) 
140 KON) SKXE( 1) 
141 MCN? 1) SKE(2)°D°7)8K(1) 
142 M(N*2)2KEB(3)-KE(1)-0(7)8K(2)-0(8)EK(1) 
143 DO 49 I24,N 
144 XM2BKM1 
148 XM1= XM 
146 JEN 1+) 
147 MMExX (ds) 
148 K2t-((3-1)/3)23 
149 GO TO (44,45,46),K 
180 44 Ma) 
1$1 Mis3 
182 M222 
153 GO TO 47 
154 4S Mz2 
185 Mir=zt 
1S6 M223 
187 Go TO 47 
188 46 mMe3 
1$9 Mi=2 
160 M221 
161 87 XC(M)SNM-KXM2-0(1)8xXC(M1)-0(2)2xXC (M2) 
162 XD(M)=XC(M)-xXOC(M2)-0(3)3KXO(M1)-0(4)2KD(M2) 
163 KXE(M)2=xXO(M)-XO(M2)-0(5)2=KE(M1)-0(6)2xXE (M2) 
164 89 K( J) SKXE(M) = KE(M2)-0(7)2K( 91) -068) EKX( U2) 
165 IF (1G.NE.1) RETURN 
166 oo SS 121,N 
167 $8 X(1)=2xX(1)/G 
168 RETURN 
169 ENO 
ENO OF FILE 
ate SUBROUTINE PLTTRCO(IV,N) 
2 c 
3 Cc WRITTEN BY CECILIO J. REBOLLAR SUMMER 1979 
4 (3 
5 C THIS SUBROUTINE PLOT A SINGLE TRACE i. @. A SISMOGRAM 
6 (s 
7 INTEGER#4 IV(N) 
8 REAL*4 FL(39000) 
r) DIMENSION ITITL(7) 
10 DATA ITITL/72° OY) 
1 WRITE(6, 100) 
12 100 FORMAT( ‘ENTER TITLE (28 CHAROS OR LESS)’) 
13 READ(S,101) ITITL 
14 101 FORMAT (7A4) 
1s CALLY PLO T2).,. 0%. :-'3)) 
16 CALL SYMBOL(-1.,0.,0.25,ITITL,930.0,28) 
1G) oo 2 181,N 
18 2 FL(I)sIV(1) 
19 CALL USMNMX(FL,N,1,FMIN, FMAX) 
20 Yzo. ; 
21 XBIV(1)/FMAX 
22 CALL PLOT(X,Y,3) 
23 DO 1 182,N 
24 YaV+.0027 
28 xBIV( 1) /FMAX 
28 1 CALL PLOT(X,Y,2) 
27 RETURN 
28 ENO 
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39 
100 
101 
102 
ENO OF FILE 


N 


( 


1 


zass 
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c CA 


100 


105 
110 


17s 
120 


129 


ana 


200 


209 


ana 


300 


409g 


IMAGINARY PARTS OF THE COEFFICIENTS ARE 
N/2*2 TO N ONO ARE NOT STORED FOR OC OR NYOQUIST FREQUENCY 
(WHERE THEY ARE ZERO). THE OUTPUT (INPUT) ARRAY HAS 
AMPLITUDES IN THE FIRST N/2*1 POSITIONS ANO PHASES 


FOR DC AND NYSUIST FREQUENCIES ARE RETURNED 
IN POC ANO PFN 


SUBROUTINE AMPPHZ(LOG2N,X,1C,POC,PFN) 


WRITTEN BY DAVE GANLEY ON MAY 18, 1977 


THIS SUBROUTINE CALCULATES AMPLITUDE AND/OR PHASE 


SPECTR& FROM THE REAL ANO IMAGINARY PARTS OF A 

FREQUENCY SPECTRUM OR VICE VERSA. THE INPUT (OUTPUT) 
ARRAY xX 1S OF LENGTH 2%*LOG2N WITH THE REAL PARTS OF 
THE FOURIER COEFFICIENTS FOR FREQUENCIES ZERO TO 


YOUIST IN THE FIRST N/2%1 POSTTIONS OF X THE 


IN DEGREES) IN THE LAST N/2-1 POSITIONS. PHASES 


c=1 INPUT IS FOURTER COEFFICIENTS. 
OUTPUT IS AMPLITUDE ANDO PHASE. 

=2 INPUT IS FOURIER COEFFICIENTS. 
OUTPUT IS AMPLITUDE ONLY. 

=3 INPUT IS FOURIER COEFFICIENTS. 
OUTPUT IS PHASE ONLY. 

=-1 INPUT AMPLITUDE AND PHASE. 
OUTPUT IS FOURIER COEFFICIENTS. 


RBESTSSSTCTSSS SESE SER ERTS ASSERTS AESTASTTTSESSETARESEsE 


DIMENSION X(1) 

M21 

IF(LOG2N.LT.1) RETURN 
00 3 122,LOG2N 

MSZ2=M 

NZ=25M 

MP 1=Me1 

PRLCRE LE Tat ie (GOR snG: 4 O0 
GO TO (100,200,300),1C 


LCULATE AMPLITUDE ANO PHASE SPECTRUM 


VEO GE. 107.0))) “CSF TO 10S 
M(1)E*eKE1) 
Poc=180.0 
Go TO 110 


Poc=0.0 
IF(X(MP1).GE.0.0) GO TO 115 
X(MP1)5°-X% (MP1) 


PFNZ180.0 

Go TO 120 

PFNZO.0 

IF(LOG2N.EO.1) RETURN 

oo 129 122,M 

Azx(1) 

Bex(1+m) 
X(1)=SORT(APA+B2B) 
XC1L+M)257. 298782 ATAN2(B,A) 
RETURN 


CALCULATE AMPLITUDE SPECTRUM ONLY 


LE GCL SIT OO IT eC eT 7 
IF(X¢(MP1).LT. 0.0) X(MP1)2-K (MP1) 
IF(LOG2N.E0.1) RETURN 

060 209 1=2,mM 

B=xX( 1) 

B=x( 1+M) 

X(1)=2SORT(AtA+B2B) 

RETURN 


CALCULATE PHASE SPECTRUM ONLY 


LECXG tito. OF 1660 16) S65 
Poc=0.0 

GO rT Om3 10: 

POC=180.0 

PEGCKCMP 1 GL TOO} GO TO.315 
PFN=0.0 

GO TO 320 

PFN=180.0 

IF(LOG2N.EO0.1) RETURN 

DO 329 l2z2,.M 

Bex( 1) 

Bex( lem) 
X(1*M)2S7.298782ATAN2(B,A) 
RETURN 


CALCULATE FOURIER COEFFICIENTS FROM AMPLITUDE AND 


X(1)2"(1)8COS(POC/57.29578) 
X(MP1)=2X(MP1)28®COS(PFN/S7.29578) 
IF(LOG2N.E0.1) RETURN 

DO 409 1=2,M 

Az=x(I) 

Bex( 14M)/57.29578 

X(1)2&2COS(B) 

X( 14M) 2a2S1N(8) 

RETURN 

Eno 


IN POSITIONS 


BECTRTT ST STEAL RTKASRRSARSESEKRTESEESES SST ATTESERRETTESAesseETE 


(SPECIFIEDO) 


PHASE 
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110 
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~ SUBROUTINE OONIEL(X,N,FWINO, DT) 


10 
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SUBROUTINE BY DAVE GANLEY DECEMBER 12, 1977. 


THIS SUBROUTINE CALCULATES A DANIEL-LIKE SPECTRAL 


STIMATE FROM THE PERIOGOGRAM THIS IS DONE 8Y 


AVERAGING ALL 
OSBOUT THE DESIRED FREQUENCY 


SINCE REAL ANO 


FREQUENCIES WITHIN A WINDOW CENTERED 


IMAGINARY COEFFICIENTS (OR AMPLITUDE AND PHASE VALUES) 


ARE AVERAGED THIS 
STIMATE. 


IS NOT A DANIELL POWER SPECTRAL 
ALSO NOTE THAT AT FREQUENCIES WITHIN HALF 


OF THE WINDOW WIOTH OF OC OR NYOUIST FREQUENCIES THE 
ESTIMATES ARE CALCULATED BY AVERAGING OVER FEWER 
VALUES 


INPUTS: 
X = INPUT FOURTER TRANSFORM AS OUTPUT BY FFTRE OR 


AMPPHZ. ON OUTPUT X CONTAINS THE SMOOTHED 


SPECTRAL ESTIMATE 
N = LENGTH OF THE X (N 
FWINO = WINDOW WIDTH I 
OT = SAMPLE INTERVAL, 


DIMENSION X(N) 
OIMENSION Y(S12) 

LPN S252.) Gu Ou toO 
WRITE(6,1) 

FORMAT(’ WN CANNOT EXCEE 
STOP 64 

N22N/2 


N222N2%2 

NP2ZENe2 

OF =S00.0/(WN20T) 
M22FWINO/(2.020F)+.5 
M21=M2¢1 

M=22=M21 

XM=eM 

XM1=M-1 

FWINO=2 .O2M220F 
WRITE(6,2) M, FWINO 
FORMAT (i= 4) 38% POUNTS 
‘WHICH IS & WIOTH OF ° 
IF (M2.E90.0) RETURN 

00 189 ITz21,N 

Y(T)SxKeT) 

KE=0 

00 49 [21,M2 

1I=N22-1 

L=I-N2 

LL=1L1eNn2 

JESM21-1 

1212241-1 

11212N22-121 
L21zI1214N2 
LL2Z1=NP2-121 

MPoMM-JE 

XFYEKF*1 0 

LEC €O. 1 (GS) To Je 
KEsi-1 

00 29 K=1,KE 

KG Ee MeV TG Le i ese (iT: Ae 
KRCTITISKCTI YETI OK) ove 
M(LIEXCLI*VOL eK) 
M(LLEISRCLEIFYCLE-K) 
LEK EO. KE) SO tG (23 

DS Ot RE > Hf a ag 8, I SC 

MCLE PEXCLLI*YICLL +K) 
CONTINUE 

00 38 J=!1,JE 
M(T)SM(1T)4Y¢C 12145) 
ROE GT Te ¥ ot Fez tice Sa, 
TECL:. EO. ty Go To 39 
MCE VSR EL YS V2 Pe) 

iQ eee Ley kk 2 ie shy 
CONTINUE 
XC TISEKCF 
XC TIT) =X ( 
LPC Tre Or 
KCLISKCLI/KEY 
MCCE MAK OCE )/ XE 1 
CONTINUE 

IES=N22°-M21 
IF(M21.GT.I1E) RETURN 
LEM2Z214N2 

00 S93 121,M2 
M(M21)SX(M21)4Y (M2141) 
IF(M.GT.N2.4N0.1.EO0.M2 
RCs Moe ya ero 

DF “CT £O.M2) GO T0538 
CH se NCE eco") 
CONTINUE 
X(M21)50(M21)/%M 
M(L)SKC1)/XM1 

IF (M.GT.NZ) X(L)EX(L) 
IF(IE.EO.M21) RETURN 
LElEe*n2 

00 69 1=1,M2 
X(LEDVSXCIE)+V(LE+l ) ov 
XE aCe Lay Miler iy) 

IF (IT. NE.M2) X(L)EKCLI+ 
CONTINUE 
X(TE)2X(1E)/X™M 
X(LISK(L)/XMI 

VES he<1 

IsS=M21+1 

IF(IS.GT.1E) RETURN 

pO 63 r=tSs7 Le 

L=l+N2 

oo 785 J=1,M2 

MC TIEKXCTDI+V( Tey ev ied 
XCLIEXCLIFV (LES OV( Lod 
CONTINUE 

M(T)SX(1)/XM 
M(L)BKOL)/XM 

CONTINUE 

RETURN 

ENO 


IN THE SAME FORMAT. 
1S 222M WHERE M IS AN 


INTEGER) 


N FREQUENCY DOMAIN (HERTZ) 


OF ORIGINAL TIME SERIES 


(MSEC) 


0 S12 IN SUBROUTINE DANIEL’) 


USED IN DANIEL WINDOW 
Rey Custer ERitce: 71): 


\ 


Ge is Bt 
+V(M21-1) 

») GO TO s9 
=zXM1/(XM1-1.0) 
Le = I) 

VCS) 


) 


} 
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REAO THE DATA RECORDECS 


Wr 


- SUBROUTINE GETDATI(2,NS, EWN, NBLKS) 


THIS SUBROUTINE WAS WRITTEN BY CHARLIE MCCLOUGHAN 


REAL 2(1020),NS (1020) , EW/ 1020) 
INTEGER LNR 

INTEGER®2 10(510),LN,18(3,1020) 
NBLKSS=NBLKS2170 


M20 


00 2 [=1,NBLKS 

CALL READ(ID,LN,O,LNR, 4) 
OO 0:2 Kis 5 S5:10.3 

N=O 

Meme] 

KK=K+2 

00 3 L2K,KK 

NSN] 

IS(N,M)IZIO(L) 
CONTINUE 

NzO 

00 4 I21,NBLKSS 
1FLz0 

oo § K=1,3 
PESTS Ly? 
ITSI1ABS( IT) 
VECUT. CT. 2683) GO To § 
IFLZIFL +1 

CONTINUE 
TECLEL NE. 3) GO! +6. 6 
GO TO @ 

NEN) 

SCN RONS Car Ty 
NS(N)2FIS(2,1) 

Ew (nd sIS(3,1) 
CONTINUE 

RETURN 

ENO 


‘SUBROUTINE CETSAM(N, NOCS,ISaM.a) 


IN THE SPREGNTHER DR-100 


THIS SUBROUTINE DECOD THE MAGNETICS TAPES OF RESMAC 


INTEGER SHFTR,SHFTL,ALL,IL,IR, IMA, IMAN, LEXP, LNUM 


INTEGER=2 


1S(512},I1E&( 40000) 


REAL ISAM‘ 40000) ,A( 10000) 
INTEGER#2 LX 

INTEGER ISUB,N 

ISUB=0 

00 3 12=1,N8CS 

CALL READ(IS,LX,0,LNUM, 8) 
OO 2 Jv=1,256 

ISUB=SISUB41 
TE(ISUB)=sISiau) 

CONTINUE 

CONTINUE 

oo 1 =1,N 

ALL2=IE(I) 
ILS=SHFTL( ALL, 16) 
ILS=SHFTR(IL, 24) 
IRSSHFTL(ALL, 24) 
IR=SHFTR(IR,24) 
IRSSHFTL(IR, 8) 

IMA=IL+IR 

IMANSIMA/4 
ITEXP=SHFTL( IMA, 30) 
IEXPZSHFTR(IEXP, 30) 
TSAM(1)=( IMAN-2048) 242218 XP 
A(l1)=!I 

CONTINUE 

RETURN 

ENO 
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(SUBROUTINE CONOT CTS NTS stSiIND 


THIS SUB REMOVE THE TRENO(LINEAR), THE AVERAGE, ADO ZEROS 


AN EVEN NUMBER ANDO SMOOTH THE SIGNAL AT THE ENOS WITH A 
COSBELL SEE KANASEWICH’S BOOK PAGE 38, 
IT CALL SUB KOSBELL 


REA SSCS 20 Sr. Sile): 
TRENOZ(TS(NTS)-TS!/1))/(NTS-1) 
BO? Sst INNS 
TS(1)2TS(1)-TRENOZ( 1-1) 
sSuM=0.0 

OOS) ferns 

SUM=SSUM+TS(1) 

SUM=SUM/NTS 

00 6 J=1,NTS 

TS(1)2TS(1)-SUM 

CALL KOSBEL(TS,NTS,NTS/SO) 
NZEROZ(N-NTS)/2 

OO Ss rst 

TS1(1)20.0 

INZEROSNZEROSNTS 
TF(1.GT.NZERO.ANO 1] .LT.INZERO) TS1(1)2TS‘1-NZERO) 
CONTINUE 

RETURN 

ENO 


~ SUBROUTINE ROTATE (WE,SN NN, AZIM! 
WRITTEN BY CECILIO J. REBOLLAR SUMMER 1979 


THIS SUBROUTINE ROTATE EW AND NS COMPONENTS IN ORDER 
TO GET THE RADIAL SND TRANSVERSE COMPONENTS OF THE 
SEISMOGRAM 
WEZRADIAL COMPONENT SNSTRANSVERSE COMPONENT 
AS THE OUTPUT 


REAL *4 WE(6000) ,SN(6000) ,A2 
REAL SAZ,CAZ 
CONV=57.295778 
SA22S5.(AZ2IM/CONV) 
SAZ=COS(AZIM/ENNV) 

Ooo 1 I=1,NN 
WE(T)=WE(1)*SQ2+SN( 1) -CA2 
SN(T)SSN(1T)8SA2-wWE(1)2Cac 
CONTINUE 

RETURN 

ENO 


SUBROUTINE KOSBEL(TS,NI,M) 
WRITTEN BY CECILIO J. REBOLLAR SUMMER 1979 


SMOOTH THE SIGNAL WITH A COSBELL (SEE KANASEWICHK B00K) 


REAL TS(512),P1/3.16159/ 

LPM2N1/2 

IF(2*tPM.NE.NI) GO TO 2 

Go -710' 3 

NISNI4+) 

TS(NI)=20,.0 

LPMENI/2 

L=LPM-m™ 

OO eet eye NYE 

CBEL=1 0 

NNWEZ1-LPM 

Tie Melee (Eos My CBEL=(1.+COS(PI®(NN#L)/M))/2 
IF(T.GT.2"L+m) CBELE(1.+COS(PI=(NN-L)/M))/2 
TS(IT)STS( 1) SCBEL 

RETURN 

ENO 


TO 
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199 
200 
201 
202 
203 
204 
205 
206 
207 
208 
2093 
210 
211 
212 
213 
214 
218 
216 
aie 
218 
2193 
220 
221 
222 
223 
224 
22s 
226 
227 
228 
229 
ENO OF FILE 


ee itQunerun— 


10 


ENO OF FILE 


OO dM KMeUn — 


21 
ENO OF FILE 


c 
c 


anan 


anaan 


FAST FOURIER TRANSFORM SEE CLEARBOUT BOOK PAGE 12 


301 


401 


501 


100 
64 


THIS SUBROUTINE OBCOD THE MAGNETICS TAPES OF RESMAC 


CALL FORK(M,C,1.) 
MMeEM/2 
C(mM41)20(1) 

dO 301 I=2,MM 
AR( 1) SRBEAL (C(I) 
AI(1)SAIMAG(C(I 
AMPOUT(I)sSORT ( 
CONTINUE 

DO 401 I22,MM 
VZAMAGN(FF(1)) 
ATENUAZSEXP( (PISHOIS#FF(I))/(O2#VD) ) 

AMP (I) SAMPOUT( I) Zs ATENUA 

AMPOUT(1)2( AMPOUT(I)*ATENUA)/V 

CONTINUE 

60 S01 I122,MM 

AMPLOG(1I)20LO0G10(AmMPOUT(!)) 
FFLOG(1)S2ALOGIO(FF(I)) 

NUMM2SMM- 700 

NUMM2NUMM- 20 

NUMMM® NUMM- 10 

WRITE(93,12) NUMMM 

WRITE(9,101) (FFLOG(I), AMPLOG(1I),1210,NUMM) 
WRITE(11,12) NUMMM 


) 
») 
AR(1)SAR(IT)*AI( 1) #Al(1)) 


WRITE(11,167) (FF(I), AMP(I),1210,NUMM) 
WRITE(12,12) NUMMM 

WRITE(12,101) (FF(I), AMPOUT(I),1210, NUMM) 
FORMAT(1X,10110) 

sTop 

END 


“SUBROUTINE GETSAM(N,NBCS, ISAM, A) 


WRITEN BY CECILIGO J. REBOLLAR 


INTEGER SHFTR,SHFTL,ALL,IL,IR, IMA, IMAN, [IEXP, LNUM 


INTEGER#2 18(512),1E( 40000) 
REAL ISAM( 40000) ,A( 10000) 
INTEGER*®2 LX 

INTEGER ISUB,N 

ISUBZ0 

00 3 I121,N8CS 

CALL READ(IS,LX,0,L NUM, 8) 
00 2 J#1,256 

ISUBSISUB+! 
IE(ISUB)z=IS(J) 

CONTINUE 

CONTINUE 

00 1 I21,N 

ALLZIE(1) 
ILSSHFTL(ALL,16) 
IL2SHFTR(IL,24) 
IRSSHPTL(ALL,24) 
IRSSHFTR(IR,24) 
IRESHPTL(IR,8) 

IMAZILSIR 

IMANS IMA/4 
TEXPSSHFTL( IMA, 30) 
TEXPSSHFTR(IEXP, 30) 
ISAM(1)2( IMAN-2648) 2422 1EXP 
A(I)s!I 

CONTINUE 

RETURN 

ENO 


FUNCTION AMAGN( FRE) 


THIS SUBROUTINE CALCULATE THE INTRUMENT RESPONCE OF 


CERRILLO STATION THE AUTPUT IS IN COUNTS/CM 


WRITEN BY CECILIO Jv REBOLLAR 


REAL PI/3.141882684/ 

AMP82.52E+413 

We2=ePIl FRE 

WF22"P1*10.0 

HFsO.707 

HO#0.7 

WOzZePIz1. 

TEMs ( (WFKee2-wee2)ee2-( 28H eWR ew) 222) 222 
TEMS TEM? (4. 2 (WKex2-Wee2)es22 (2H EWF ew) 22) 
TEME (W223) SORT (TEM) 

TEMDZSORT ( (WO 2-W222)822+( 28HOPWOEW) 222) 


TEMD=TEMD= ( ( (WFEe2-We22) 222+ ( 2eHF EWR ew) 222)222) 


AMAGNe ( AMP®TEM) /TEMD 
RETURN 
END 
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